








First Dissociation Constant of Phosphoric Acid From 
O° to 60° C; Limitations of the Electromotive Force 
Method for Moderately Strong Acids 


Roger G. Bates 


Ar ac irate determinatior lbw the electromotive-force method of the dissociation 

stant of an acid as strong as the first acidic group of phosphoric acid requires that the 
hvdrogen-ion concentration be established with greater accuracy than is usually attainable 
Phe uncertain quantity is the activity coefficient of hydrochloric acid in the particular 
iixture nder study It is shown that the calculation of this coefficient by the Debye 
Hiickel equation vhich affords an adequate estimate of the ma hvdrogen-ion corre 
in buffer solutions composed of weak acids, fails to vield the desired accuracy for mixture 
in which the acid is dissociated to an appreciable extent In order to determine the commor 
warithm of the first dissociation constant of phosphoric acid with an accuracy believed to 
bee 0.005 unit, the hydrogen-ion concentration of some of the buffer mixtures studied wa 
kept at low levels by use of a large excess of potassium dihydrogen phosphat Ont , 

itions contained formic acid instead of hydrochloric acid, thereby lowering the hydrogen-ion 
onecentration still further In addition particular attention was given to a choice of reasor 
able values for the activity coefficient used to estimate the hvdrogen-ion correction 

By these means the first dissociation constant, A,, of phosphoric acid was calculated 
from emf measurements of 71 phosphate buffer solution together wit! ipplementary 
data from the literature for 18 additional solutior Phe re t from O° to 60° C is given b 

709 3) 
log A 1.5535-+- 0.013486 7 
] , 

where 7 the temperature in kK The changes of heat content entropy ind heat capacity 
for the dissociation of | mole of aqueous phosphoric acid in the standard state were found 
to be 7.650 66.8 j des and 154 5 deg respectively, at 25° ¢ 

l. Introduction peared unlikely that the dissociation constant of an 


acid as strong as the first acidic yroup of phosphoric 
acid could be fixed with an accuracy of 0.01 in the 
common logarithm by this emf method, and the 
study was accordingly abandoned The cause of 
the difficulty was readily apparent; the activity co 
eficient of hydrochloric acid in the buffer solutions 
must be known accurately in order to obtain the 
hydrogen ion concentration and consequently the 
molalities of phosphoric acid and primary phosphate 
ion from which the buffer ratio is computed. This 
activity coefficient cannot readily be determined 
lions of the apeuti and molecular Spots becomes -~ exactly, and use Is usually made ol estimated values, 
arge that an exact determination of the dissociation coleaiated perhaps by some form of the Debve- 
onstant is no longer possible Hiackel ’ 


experience has proved that the electromotive- 
rece method is well suited to the accurate determi 
tion of the dissociation constants of many weak 
ids and bases in aqueous solution. However, it is 
not cenerally recognized that the accuracy obtainable 
by this method diminishes with increasing strength 
of the acid or base. When approximately 10 percent 
more of a weak acid is dissociated in solution 
the uncertainty in establishing the true concentra- 


equation This device Is completely ade 
Some years ago the author undertook to measure 


. . quate for acids so weak that the total contribution 
the electromotive force of a series of cells of the ty pe 


of the hydrogen ion is small but is unsatisfactory, 
in general, when the pH of the buffer solution lies 


Pi H, A l atm), KH PO, ri), HC] mi AoC |: Ag below 2 5 to 30 
Cell 1) It was realized later that two possibilities of re 
ducing the hvdrogen ion concentration might be 
vhere m is molality, in an attempt to extend the in- | feasible, and the investigation was resumed. The 


vestigation of Nims [1]' to other temperatures and | first and simpler approach is to increase the value 
buffer ratios It was found to be impossible to | of 7 as greatly as is consistent with adequate buffer 
obtain an accurate value for the first dissociation | capacity and to estimate the activity coefficient 
constant of phosphoric acid by analysis of the data | from independent data. The second is to substitute 


for mixtures in which 1 r <3. At the time it ap- | a weak acid of known strength for the hydrochloric 
acid which, together with potassium dihydrogen 
Figures in brackets indicate the literature refereaces at the end of this paper. | Phosphate, forms the buffer solutions, at the same 
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time adding neutral chloride to furnish the desired 
concentration of that ion. This latter procedure 
complicates the analysis of the data, for the hydro- 
gen-ion concentration is then fixed by two dissocia- 


tive equilibria instead of one. However, it is pos- 
sible to evaluate the constant for one of these 
equilibria when that for the other is known [2, 3, 4] 


The results of both approaches are presented 
Data for 46 mixtures of potassium dihydrogen 
phosphate and hydrochloric acid, in which zs varied 
from 1.2 to 6, are reported. Twenty-five mixtures 
of potassium dihydrogen phosphate and formic acid 
and 8:1) with added potassium chloride 
were also studied. The combined data, as well as 
those for 18 solutions studied by Nims [1], were used 
to evaluate A, for phosphoric acid from 0° to 60° C 
with an estimated accuracy of +0.005 unit and to 
demonstrate the manner in which the accuracy of the 
electromotive-force method depends upon the 
strength of the acid studied. 


(ratios 5:1 


2. Experimental Procedures and Results 


Two different preparations of potassium dihydro- 
The first was prepared 


gen phosphate were used 
two recrystallizations 


from reagent-grade salt’ by 
from water. The salt was dried at 110° to 130 A 
determination of the weight suffered by 
weighed samples on ignition to metaphosphate indi- 
cated that the dried salt retained than 0.02 
percent of water. This portion of salt was used to 
prepare the 32 solutions grouped in series 1 through 
4 (see table 1). The second preparation, from which 
the solutions of series 5 through 7 were made, was 
NBS Standard Sample 186—Ib. 

The potassium chloride was a fused sample free 
of bromide [5]. Reagent-grade 90-percent formic 
acid was distilled twice in an all-glass still. Only 
the middle fraction was retained for redistillation or 
for use in the preparation of the solutions. Hydro- 
chloric acid was diluted to a concentration of 6 molar 
and distilled in a similar fashion. An approximately 
0.1-molar solution of the distilled acid was stand- 
ardized gravimetrically by weighing silver chloride. 
The specific conductance of the water used in pre- 
paring the buffer solutions did not exceed 1.0x10~ 
ohm™' em In view of the low pH of all of the 
buffer solutions, no special precautions were taken to 
exclude carbon dioxide, although dissolved air was 
removed by passage of nitrogen through the solutions 
before the cells were filled. In other respects, the 
experimental procedure was essentially the same as 
that followed in the determination of the second dis- 
sociation constant of phosphoric acid [6]. 

The solutions prepared from potassium dihydrogen 
phosphate and hydrochloric acid were grouped into 
five series according to the ratio of the molalities of 
the two components, that is, Mxn,~o!%ac1. This 
ratio was 1,2414:1 for series 1; 1.2683:1 for series 2; 
1.5581:1 for series 3; 2.184:1 for series 4; and 6:1 for 
series 5. The cells containing the aqueous mixtures 
of potassium dihydrogen phosphate and formic acid 
can be represented as follows: 


loss in 


less 


rm), HCOOH (m) 


Pt; H, (g, 1 atm), KH,PO, ’ 
KCI (me,=0.005), AgCl; Ag. (Cell ID) 
The molality of potassium chloride w as 0.005 in each 
solution, but the ratio of the molality of potassium 
dihydrogen phosphate to that of formic acid was 
5:1 in series 6 and 8:1 in series 7 
The measured electromotive force of each cell was 
corrected to a hydrogen partial pressure of 1 atm in 
the usual way. The pwH value of each of the 7) 
buffer solutions at each temperature of measurement 
was computed by the defining equation * 


K—E°)F 


. oe | log i ‘ 
2.3026 R7 Emer, (1 


log liu i Ny 


pwH 


in which f represents a molal activity coefficient; E 
and £° are the electromotive force and standard 
potential, respectively, of the cell at the appropriate 
absolute temperature, 7; and FP and F are the gas 
constant and the faraday Values of E° and 
2.3026RT7/F in absolute volts are summarized in an 
earlier publication [7]. The pwH values are listed in 
table 1. 


3. Calculation of the Dissociation Constant: 
Series 1 to 5. 


3.1. The Conventional Electromotive-Force Method 


Inasmuch as the hydrogen-ion concentration in 
the buffer solutions of series 1 to 5 can be expressed 
in terms of the first dissociation equilibrium of phos- 
phoric acid, pwH can be related directly through 
the mass law with A,, the first dissociation constant, 
and the molalities and activity coefficients of phos- 
phoric acid and primary phosphate ion. Making 
the substitution in eq 1 and rearranging gives the 
formal expression for log Ay: 


MHPO, In,po,! , ! » 


/ 0 
J HyPO, 


- log kK, pwH T log 


m > 
HPO, 


This substitution is the basis of the method of de- 
termining the dissociation constant, as outlined by 
Roberts [8] and by Harned and Ehlers [9]. The 
solutions of series 1 to 5 (cell 1) were composed of 
hydrochloric acid, m moles /1,000 g of water, and 
potassium dihydrogen phosphate, zm moles/1,000 g 
of water. Because of the appreciable dissociation of 
the phosphoric acid formed, the reaction between 
these two components is incomplete, and at equili- 
brium 

Mu,Po, m-—Ma, ») 
and 
+- My. (4) 


Mu,vor = M(x —] 

rn : : 
The hydrogen-ion concentration needed to com- 
pute the concentrations of the two phosphate species 
2 For convenience, the charges carried by hydrogen (H*) and chloride (CI-) ions 


will not be indicated. Likewise, mu will represent the total concentration of 
hydrogen ion, whether hydrated or not. 
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on the right-hand side of eq 
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lhe 5 solutions for which pwH from 0° to 50° is listed make up series 5a; the remainder are in series 5b, 


5a; A and 
temperature [10] 





th 





Tasie | 


pwHl from 0° to 50° for aqueo mixtures of potass 


rlures of potlass ” fil 


KH.PO, He 


KHLPO 


and a* is the ion-size parameter. The ionic strength, 


u, of these mixtures ts given by 


uu rm muy 6 


The Debye Hike kel equation Ssuyvests that the 
last term of eq 2 will be zero, or nearly so, in the 
region of low tonic strengths Terms of this sort are 
indeed found to be small at low concentrations and 
to vary linearly with ionic strength in the region of 


moderate concentrations. The evaluation of log K 


is accordingly made by extrapolation to infinite 
dilution of the apparent log A,, or log At, derived 
from the following equation 
; ” m Y 
log A pwH + log i 
mia l Min 


3.2. Limitations of the Conventional Method as 
Applied to Moderately Strong Acids 


The procedure outlined in the foregoing section 
may be termed “the conventional emf method.” 
The choice of a* in 5a is somewhat arbitrary, 
being guided only by our knowledge of what consti- 
tutes a “reasonable”’ ion size, for the magnitude of 
this parameter is often of the same order as the ionic | 
diametey in angstrom units. If —log A for the acid | 


eq 


exceeds 3, the particular choice of a*, within the 
gos Im , where m, and ire, respect ty and charge of the 
ion of species i 


mwmic acid 


hydrogen phosphate and hydrochloric acid and for a 


and potassium chloride Continued 


oOoH K« . 


HCOOH «s Ke 


2 tos, for - 
is small compared to m, and relatively large errors 
IT Wy alter the last term of eq 7 but slightly How 
ever, When the acid is as strong as the first dissocia 
tion step of phosphoric acid, for which log A, is 
2.12 [1], my is quite large and the choice of 
a* becomes an important consideration 

Figure | demonstrates the failure of the extrapo 
lation lines computed with various values of a’ 
an ionic strength of zero. The oper 
circles are the data at 25° for series 1, the mos 
acidic solutions included in the study. The values 
of a* are indicated at the right of each ling 


reasonable range is of little concern 


to 
converge at 


eq ja 
When a* 0, the last term of eq 5a reduces to the 
Debye-Hiickel limiting law. The dots in figure | 
were computed from the data of series 5a at 25 


The ratio of phosphate to hydrochloric acid is con 

siderably higher than in series 1. Accordingly m 

is smaller, and the choice of a* has less effect upon 
the result. It is unlikely that tne difficulty could b 
resolved and the lines made to converge by obtamineg 
data at lower ionic strengths, for the degree of dis- 
sociation of the phosphoric acid would be inevitably 
increased. Furthermore, still -higher ratios of pri 

mary phosphate to hydrochloric acid might intro- 
duce uncertainties in the extrapolation, and it seems 
advisable to avoid them. 


3.3. Estimation of the Activity Coefficient 


Various possible means of choosing the correct ex 
trapolation line were next considered. Up to this 
point the determination of the dissociation constant 


130 











rr Mig 


TrOors 
How 
“>a lil 
K, is 


ce of 


rapo 
a* to 
opel 
most 
alues 
line 
> the 
ire 

con 
ym 
upon 
lal be 
nine 
f «lis- 
tably 
. pri 
ntro 


eecims 


*t eX- 
» this 
stant 











ee 
ee 
ee 
oo 
Bo 


con pletely rigorous, and it is unfortunate that a 
rtially nonthermodyvnamiec step need be introduced 


\ eq 5 suggests, different extrapolation lines with 
llerent intercepts eNXtst because the correct value 
f., in each buffer solution is unknown. Un 


tunately, the required activity coefficients cannot 
dilv be determined Nevertheless, the most prob 
ble range of a” eq 5a) can be narrowed somewhat 
the light of the values that have been found to 
deseribe the behavior of the activity coefficient of 
drochloric acid and fy, /-, in other mixtures of sim 

r composition. Thus log fy fe; in aqueous solu 
ons of hvdrochloric acid with tonic strengths less 
n 0.2 Is represented by the last term of eq oa 

If the hydrochloric acid Is re placed 

irtially or even completely by alkali chloride with 

t altering the tonic strength, the activity coeffi 
ent term for Te decreases l 1] ' ‘te represent 
these new, lower activity coefficients, a lower value 
:* is necessary Accordingly, a* for fy fe, 1s found 
to approach a limit of 4.6 to 5.4 in pure solutions of 
mibuinn chloride and potassium chloride with bone 
The principles of tonic in 


th a*—6.0 


trengths of 0.2 or less 
teraction lead us to expect, in this instance, a large: 
pecific effect of cations than of anions. Hence, it 
is possible that the activity coefficient of hvdrochlori 
acid in a mixture of the acid and primary potassium 
phosphate does not differ greatly from that in a pure 
solution of potassium chloride. The most probable 
range of a* can be fixed at 4 to 6 and a value of 5 
chosen as a first approximation 


‘Except possibly with lithium chlor r which some measurements show a 


ia = *,2 $ e ‘ 
° Q -_ . - | 
e nee » . | 
r < a 
sts a 
» a 
] 
x } 
> | 
} 
- > + 
" ~ 
} ~ 
> 
‘ 
’ 
« 
‘ 
7 nD . 
v RENGTH 
hiaure 2 P qgK,at asafur ’ ) ength 
\ ‘ ' ] ~~ 
@ oor 
Two separate caleulations ol log A were made 
at O°, 25°, and 50° by eq 7 \ value of 4 for a* 
eq 5a) was used for the first and 6 for the second 


The data of series 1 to 5 were supplemented by those 
of Nims [1], who studied three series of six solutions 
each The ratios of the molalities of potassium di 
hydrogen phosphate to hydrochloric acid in the latter 
were respectively 2.071:1 (series N11); 3.102:1 (se 
ries N2 and ,079:1 (series N3 ‘Two or three up 
proximations were hecessary ih computing log Ay, 
as the contribution of the hvdrovgen ion to the total 
ionie strength was considerable The extrapolation 
lines for a* tare shown in figure 


’ 


4. Calculation of the Dissociation Constant: 
Series 6 and 7 


In order to reduce the hvdrogen ion concentration 
below the level in the mixtures contammimng 6 moles of 
potassium dihvdrogen phosphate to | of hydrochloric 
acid (series 5), two series of solutions were prepared 
from potassium dihvdrogen phosphate and formic 
acid Kach buffer solution was also 0.005 molal with 
respect to potassium chloride The formal expres 
sion relating the product of the dissociation constants 
for the two acids to the emf or pwHL has been derived 
in another article |4] 

Mapo ney 
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Cirel seri and dot r were Ca cperinme 
using a*=4, and tl | wer ! t k 
the shift that would r i r 1 t 


the dissociation constant of formic 


is well known from 0° to 60 IZ, 131, 
derived from the product, 


Inasmuch as Ay, 
acid (HFo 
the value of A, can be 


K Ay 


If m, represents the molality of potassium dihydro- 
gen phosphate and m, that of formic acid, 


Mupo = Mr Mr, (9) 
Mypo m my (10) 

and 
My PO, ym My My { ] l } 


Furthermore, the last term of eq 8, which becomes 
zero at infinite dilution, is a linear function of ionic 
strength at low and moderate concentrations 
Hence, log K ky can be evaluated conveniently by 
extrapolation to infinite dilution of the apparent 
quantity, log (AK, Ay)’, computed by the equation 

log (A, Ay)’ = pwH 4 


(im My my Moly Myo 


1/2 log 


an tT My Mero 


(12) 


The calculation of the formate molality rests on 
the assumed equality of fyfer and fafro /fare at 
moderate concentrations and ts hence given the prime 
mark denoting an apparent or inexact quantity: 


Myo- = may /(1 + y), (13a) 
where 
log y= pwH-+-log A, (13b) 


. oe , 
Equation 5a is used once more to compute my, the 
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Ficure 4 Effect of choice of a* on the value of oa K 


obtained f ; 


om 10 agroups of solutions 
apparent hvdrogen-ion concentration, and the ioni 
strength is given by 


Mh My My Mey m My t O.005 14 


It should be emphasized that mj, mj,, and log 
(K, Ky)’ approach the true my, Mgo-, and log KA, 
as the solution becomes more and more dilute 

The data of series 6 and 7 were used to calculate 

1/2 log (A, Ay)’ in the manner described above 
As before, two different values of a*, 4 and 6, were 
used to obtain mw by eq 5a Two successive aup- 
proximations were sufficient to establish consistent 
values ol my and . The dissociation constant of 
formic acid, Ay, was taken to be 1.638X10~* at 0°, 
1.77210-* at 25°, and 1.650X107* at 50° [13] 
The results for series 6 are shown as open circles in 
figure 3 and those for series 7 as dots. The plotted 
points have been calculated from the experimental 
data, using a*=4, and the solid lines fitted to them 
The broken lines indicate the shift that would result 
from using a*=6. The intercepts are listed in table 
2. The difference is about 0.002 to 0.003 in 1/2 
log A, Ay, or about 0.005 in log A,. The uncertainty 
in a* is expected to have a smaller effect on the 
results of series 7 (8:1 ratio of primary phosphate to 


formic acid) than on series 6 (5:1 ratio). 


Tassie 2. Summary of log K, Ky for two values of the ton 
stze parameter 
Series 6 Series 7 
t = 

a*=4 a*=6 a*=4 a*=6 

0 2 Ol4> 2 Olts 

25 2. M7. 2. 0s 2. 476 2. 499 
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5. The Dissociation Constant at 25 


lable 3 summarizes the values of the first dissocia- 
tion constant computed from all seven series and the 
data of Nims |1], a total of 81 buffer solutions The 
series are arranged in the table from top to bottom 
n the order of decreasing acidity, and it is evident, 
both in the table and in figure 4, that the uncertainty 
of the determination decreases in the same way 
For a*=4, the 10 groups of solutions furnish about 
the same value of log K The mean of these is 
»> 144. and the mean deviation from this figure is 
only 0.0014 For a*=6, the values of log A 
display a general downward trend with increasing 
ratio of phosphate salt to acid. It is not believed 
that the constancy indicates the correctness of al 
value of 4 for a*, as the parameter usually changes 
somewhat with buffer ratio. The lower end of each 


bar in figure 4 corresponds to a*=4, the upper end 
to a” A) 

The mean of the 20 values of —log A, given in 
table 3 was chosen as the best estimate of log K 


that could be made. This average figure of 2.148 

0.005 at 25° is considerably higher than 2.124 found 
by Nims [1] with @*=—0, that is, by computing 
ufc) In eq 5 from the Debye-Hiickel limiting law 
Noves and Eastman [14] have reported the conduct- 
ance of seven solutions of phosphoric acid of molalities 
ranging from 0.0002 to 0.1 From these data, Sherrill 


PARLE 3S Summary of F A, at 25° for two values of the ior 


fire pa ameter 


and Noves 15] derived a value of 2.081 for log A, 
at IS In a later recalculation, Lugg [16] obtained 
2.09. This figure is in good agreement with 2.096 
given by Nims for 18° [1]. Mason and Culvern [17] 
have recently extended the conductance Peasure- 
ments at 25 From their data for the dilute solu- 
tions, together with those of Noves and Eastman, 
thev obtained 2.126 for log Ay, at 25 However, 
there is some indication that their data, treated 
independently, would vield a value near 2.140 

Bierrum and Unmack have calculated the first 
dissociation constant of phosphoric acid trom meas- 
urements of the emf of hvdrogen-calomel cells at 
IS°, 25°, and 37° [18]. Their results are compared 
in the following table with those of Nims [1] and 
those reported here 


IS 25 7.5 
Bjerrum and Unmach 2. 120 2. 161 2. 235 
Nimes 2. 006 2. 124 2. 185 
Bates 2. 119 2. 148 2. 200 


The value given by Bjerrum and Unmack at 37° has 
been corrected to 37.5° by addition of 0.008 By 
the analysis of titration curves for phosphoric acid 
in the presence of alkali chlorides, Hentola obtained 


2.172 for log A, at 25 19] 
6. The Temperature Coefficient of log K, 


Although the choice of the correct K at 25° was 
dependent upon a partly arbitrary selection of a* 


€ Ki for (eq 5a), the temperature coefficient of log Ay was 
_ A(log & found to be substantially unaffected by this selection 
‘ a* = é ; 
: In other words, the difference between log kK, at 0 
a - mene and 50° was practically the same whether a* were 4 
14 2. 158 oi or 6, as shown in table 1, where log A log K, 1s 
4 2 16 oi . . ver 
‘ 142 2 14 O12 tabulated for seven groups of solutions lhe figures 
Nien : . on = in parentheses are the values calculated with 6 for 
. 2.141 2. 14 00 a*, and the others were obtained with a* 1 The 
r 214 2190 oo . - : 
; 4 2 148 00. values of log A, at each temperature were derived by 
: ae? ent the methods described earlier in this papel ‘| he 
Meat 2. 148 average difference at each temperature, given in the 
last column, was combined with 2.148, the value of 
Prarie 4 Change of the first dissociation constant of phosphor acid with temperature calculated with a ‘; palues for a® a 
appear in pa entheses 
h Kiel ‘ 
Ave 
h j 
, . 
t 
a] ooo ooo oro) ooo i ; 
0. 018 O18 ool 0. 016 ) ‘ 
Ww ms Ons il mz 
1 O51 om Ost mi A) ” 
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log A, at 25°, to obtain —log A, at the other 12 | be about 30 | mole On the other hand. the 


temperatures The following equation, the con- | temperature coefficient at 25° seems rather well 

stants of which were determined by the method of | established, and accordingly A//°, AS°, and A¢ 
least squares, represents the “observed” log AK, from are believed accurate within 150 } mole eS 
0° to 60° with an average deviation of 0.0005: deg mole and 5 J deg mole respectively at 
25 The uncertainty in the temperature coefficient 
, 799.3] oe : of log Ay is considerably larger near the ends of the 

log A Ri O.OL34867 » . 
] temperature range covered by the measurements than 
in the middle Hence, these three thermodynamic 
The observed and calculated values are summarized quantities are given only for temperatures from 
in the second and third columns of table 5, and K I> 15 to 45 

given in the last column The value of A/7° at 25 7,650 ) mole is to 
” be compared with 7.420 | mole derived from the 
anne 5. S ary » K, and K to 60 emf data of Nims [1, 13] and with —8,160 j mole! 
from the calorimetric data of Pitzer 21] Pitzer 

j chose 7,870 j mole”! as a “best value.”’ For th 
entropy ¢ hange on dissociation, the measurements of 
Nims lead to AS 65.7 J deg mole”! at 25° and 
those of Pitzer to 66.9. The latter figure is in 

- excellent avreement with bos |] deg mol ob 

a ae ‘ tained in this Investigation 
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Determination of Ash in GR-S Synthetic Rubbers 
and Latices 


Frederic J. Linnig, Lewis T. Milliken, and Ralph I. Cohen 


Rapid methods are described for the determination of the ash content of GROSS svn 
thetic rubbers and latices Filter paper and filter pulp are used to prevent either material 
from overrunning the edge of the crucible when introduced into the furnace at the ashing 
temperature Water is used to facilitate the burning of the carbonaceous residue obtained 
in ashing GR-S latex These procedure have been put into use in the control testing of 
GR-S production hey are as precise as older methods of ashing and far lk time-consum 


ing The techniques emploved st 


l. Introduction 


The determination of ash is an important factor 
n controlling the quality of GR-S synthetic rubbers 
ind latices and in assuring the consumer that his 
naterial meets the requirements of the specification 
Methods used in ashing natural rubber [4, 8, 9] 

matural rubber latex [7, 8, 9, 18], and other similar 
organ compounds 11, 2, 3, 5, 6] eall for distillation 
of the dry sample by gradually raising the tempera- 
ire of the specimen from room temperature to a 
point somewhat below the ignition temperature of 
the gases formed. The carbonaceous residue is then 
ashed to constant weight at temperatures varying 
from 400° to 550° C, depending on the substance 
tested. This distillation procedure is too time con- 
suming for control testing. Introduction of GR-S 
rubber into a muffle furnace at the final ashing tem- 
perature, as suggested by one of the plant laboratories 
17], vields satisfactory results with many samples 
However, when applied to certain samples of GR-S 
rubber, particularly those coagulated with alum, 
this procedure yields results that may be low by as 
much as 30 percent These low results are due to 
loss of material that froths and overruns the edge 
of the crucible during the early stages of the pyroly- 
sis. In addition, when GR-S latex is ashed accord- 
ing to the procedure used for natural rubber latex, 
the carbonaceous residue formed on distillation may 
be removed only at temperatures that cause appre- 
ciable loss of inorganic material. Accordingly, tech- 
niques were developed for safely introducing either 
undried latex or rubber into the furnace at the ashing 
temperature, and for facilitating through the use of a 
washing treatment, the removal of the carbonaceous 
material formed in ashing GR-S latex. These tech- 
niques have been adopted for use in the official 
methods of the Government synthetic rubber pro- 
gram [10, 12, 13], and should be applicable to other 
organic materials that may give rise to difficulties 
in ashing 


2. Procedure 


GR-S. The sample to be tested is mill-massed on 
a laboratory rubber mill and sheeted. A 2- to 5-g 


' Figures in brackets indicate the literature references at the end of this paper 


would be applicable to other organic materia 


specimen of the sheeted rubber is weighed to the 
nearest milligram and wrapped in a single sheet of 
l5-em “ashless”’ filter paper by spreading it along a 
chord of the circle midway between the center and 
the edge. The paper is rolled into a tight evlinder 
around the rubber, and the ends of the filter paper 
are folded over the sample. The resulting roll is 
then folded over once or twice on itself and placed 
ina porcelain or aluminum crucible that has a volume 
of from 30 to 60 ml, and that has been ignited to 
constant weight (40.2 mg) at 550° C plus or minus 
25°C. The crucible containing the wrapped speci 
men is placed directly into a muffle furnace main 
tained at 550° C plus or minus 25° C The door of 
the furnace is closed and should not be reopened 
until combustion of the carbonaceous material is 
nearly complete (30 to 60 min.). After the carbon 
appears to have burned completely, the crucible is 
removed from the furnace, cooled to room tempera 
ture ina desiccator, and weighed \s a check for 
the complete removal of carbon, the specimen may 
be reheated for 830-min. intervals until it has attained 
constant weight within + 0.2 mg 

\ correction should be made for the ash content 
of the filter paper when the effect is greater than 
0.01 percent of the weight of the specimen. 

GR-S Latex \ circle of ashless filter paper 9 cm 
in diameter is folded into the form of a cone, and 
inserted into a 60-ml porcelain crucible previously 
ignited to constant weight 0.2 mg) at 550°C plus 
or minus 25° ©. A small quantity of dry ashless 
filter pulp 0.25 to 0.50 &) 1s pressed down into the 
bottom of the cone, and the crucible and its con 
tents are weighed to the nearest milligram in a 
weighing bottle or friction top metal container Five 
milliliters of latex, previously filtered through an 
80-mesh stainless steel wire screen, are poured into 
the cone and the crucible, and its contents are re 
weighed to the nearest milligram. The crucible 
and its contents are placed directly into a muffle 
furnace maintained at 550° C plus or minus 25° C 
The door of the furnace is closed and should not be 
reopened for 45 to 60 min The crucible is then 
removed and allowed to cool until it may be touched 
by hand (about 10 min.). Two to five drops of 
distilled water are added to the ash, the actual 
number depending upon the quantity of residue 
present, Sufficient time should be allowed between 
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drops for the water to he absorbed by the residue 


An unabsorbed (X- 
it will cause 
side of the 
subsequent 


before the next drop is added 
cess of water should be avoided, sines 

of the upon the 

during ashing, thus making 
more difficult. The dampened residue 
ashed for one-half hour, and the alternate wetting 
and ashing procedures are repeated until all of the 
The crucible 
and its contents are cooled to room temperature in a 
cle siccator anda weighed \s ne he k lor the complete 
removal of carbon, the wetting tech 
niques may be repeated again until the crucible and 


spattering residue 
erie ible 
washes 


Is 


carbon appears to have been burned 


and ashing 


its contents have reached constant weight within 
+0.5 mg. This usually requires from two to five 
washes and a total time of 2 to 4 hours, depending 


upon the latex tested 


3. Discussion of Techniques 


3.1. Temperature of Ashing 


\ temperature of 550° C was chosen to allow for 
the fairly rapid burning of both uncompounded 
GR-S and GR-S latex with a minimum loss of in- 
organic material through decomposition or volatili- 


zation. Reducing the temperature of the furnace 
below this point greatly slows down the rate of 
ashing. Lower temperatures are necessary only 


when the method ts applied to certain compounded 
stocks containing appreciable quantities of materials 
that decompose or volatilize at this temperature [{14] 
The furnace door should be kept closed to prevent 
rapid burning and glowing of the specimen at 
temperatures that may be above the furnace tem- 
perature. Loss from explosion of the distilling gases 
during the early stages of the pyrolysis is prevented 
by the filter paper that appears to be intact at this 
time 


3.2. Choice of Crucible 


Both aluminum ? [16] and porcelain crucibles have 
been found satisfactory for use in the determina- 
tion of ash in GR-S. Aluminum crucibles are 
more suitable for rapid routine testing because they 
are (1) much less expensive; (2) more uniform in 
weight (approximately 5.2 g for those used in this 
work); (3) more rapidly cooled. The crucibles used 
in this work have straight sides and a relatively 
large Consequently, they are more stable 
and the wrapped specimen remains on the bottom 
of the crucible even when it tends to uncurl. 

However, only porcelain crucibles may be used in 
ashing GR-S latices. Figure 1 shows the results 
obtained when both porcelain and aluminum cruci- 
bles are It was believed that the observed 
increase in the weight of the residue in the aluminum 
crucibles on repeated wetting and ashing resulted 
from a reaction between the aluminum and _ the 
water in the presence of the alkaline residue obtained 


base 


used 


The aluminum crucibles used here were obtained on special order from The 
Aluminum Company of America in accordance with “drawing Alcoa) FR 
0747-X, capacity 6 cu. in., 2 gage Material, 14%. in. overall height,1.675 in 
diameter, beveled at bottom The price at the time of purchase was approxi 


mately $14.00 a thousand 


! 
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from ashing the soap in the latex. Tests made by 
wetting pure sodium with water and 
heating as directed in the procedure lend support 


to this view 


carbonate 


3.3. Use of Filter Paper and Pulp 


In addition to preventing loss from explosion the 
use of filter paper allows all types of GR-S polymers 
to be ashed without loss frem frothing when intro- 
duced into the furnace at the ashing temperature 
This technique has been used in ashing volatile 
metallo-organic compounds [11]. Oxalie acid 
ammonium oxalate placed over the specimen in the 
crucible help to prevent “boiling over” of the poly- 
mer, but are not as satisfactory as filter paper 
Rubber charred with nitric acid or treated with 
sulfuric acid or ammonium sulfate prior to ashing 
does not effervesce, but the resulting values are 
different from those obtained by the slow-ashing 
procedures 14, 8] In addition, the ash obtained 
when the specimen is wrapped in filter paper is quite 
porous and fluffy, and can be readily extracted or 
dissolved for further analysis. 

GR-S latex may be ashed without previously 
drying the specimen if ashless filter pulp is used to 
absorb the liquid and to prevent the spattering that 
would otherwise occur during the first few minutes 
in the furnace. Not only is a fluffy residue obtained, 
but the cone prevents the residue from adhering te 
the walls of the crucible. The washing and burning 
of the carbonaceous material are thus facilitated. 


ol 


3.4. Washiag of the Ash 


It is believed that the slow and incomplete ashing 
of GR-S latex is due to a protective layer of salt 
that covers the carbon formed when the polymeric 
portion of the latex is burned. It is probable that 
treating the residue with small portions of water 
removes at least part of this coating of salt, permit- 
ting an additional portion of the carbon to be burned. 
Qualitative tests indicate the presence of carbonate 
resulting from ashing the soap present in the original 
latex. Since the coagulated polymer contains but 
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little soap, GR-S may be ashed within | hour with 
out the washing treatment. Nitric acid, sulfuric 
acid. hydrochloric acid, and aqueous solutions of 
oxali acid have been used to Wash the residue 
Onlv one treatment with any of the mineral acids 
equired to vield a carbon-free ash, but since 
hese acids alter the composition of the ash, the 
4 ilting values though reproducible ure different 
from those obtained when no treatment is used 
co oxalic acid offered no advantage over wate 
vater Was chosen This technique of wetting the 
rbonaceous residues with water to facilitate burn 
s simpler than a number of previous methods 
t) 
| ine 1 shows the effect of using two different 
Appar 
tly, the use of one-half-hour heating periods reduces 
time required for ashing by 50 percent. Furthes 
luction of the heating periods is not helpful 


ods of heating between water washes 


4. Results 


In order to determine the over-all precision of the 
when used over an extended period, duplicate 
determinations were made daily for 5 to 10 days on 
uples of GR-S and GR-S latex. The particular 
mples of Gk Ss used were chosen because ol thei 
marked tendency to overrun the edge of the crucible 
vhen rapidly ashed without the use of filter paper 
lhe samples of GR-S latex used were chosen because 
they characterize differences in latex types that 
might effect either the final results or the ease of 
testing. Both porcelain and aluminum crucibles 
were used in testing the GR-S rubber samples so 
that the values obtained might be compared For 
further comparison the same rubber samples were 
tested according to the slow-ashing process |4, 8] 
without previous extraction 
Che results of these tests are given in tables 1 and 
An analysis of variance [15] made of the data 
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showed a significant, but small, day-to-day varia 
bility for the tests used on both the rubber and latex 
This variability was not sufficiently large to affect 
appreciably the values for the over-all precision, s, 
that is, the standard deviation of a single random 
determination given in the tables. The values of xs, 
for each of the tests are of the same order of magni 
tude as that obtamed when GR-S ts analyzed by the 
older slow-ashing process In addition, the data 
obtained for GR-S rubber indicate at least as good 
retention of material by the new process as by the 
old A similar comparison for GR-S latex cannot 
be made, since the older slow ashing methods are 
not applicable to this material. Though the ash 
is not necessarily identical with the morganic ma 
terial originally present in the sample, because of the 
decomposition of SOPs ole the residue represents 
the origmal adequately for specification purposes 
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Opening of Testing Sieves' 


Frank G. Carpenter ’ and Victor R. Deitz 


It i hown tl 
than the average 
opening by meal 
type used for high 
1.000 microns (I 
spheres, the prepa 
measurement of tl 
a statistical analy 
method were four 


at the ‘effective’? opening n 
pening A method of measuring the effective 


opening or the nominal 


of testing sieves is generally somewhat larger 


of a calibrated mixture of glass spheres is described Cla beads of the 
wav markings were a suitable source of gla spheres for size from SO to 
S. Sieves No. 170 through 18 Details are given for the selection of the 
ration of the desired mixture of spheres, subdivision of the sample the 
e diameters, and the calculation of the weight-size distribution From 


of the calibration data 
din general to be about 1 


l. Introduction 


ed both for co 


tribution, testing sieves are 

Sieve analyses obtained by different labora 
es (for example, those of 
disagreement The dis 


even the same sample are 


The sper ifications for practically every powdered 
ranular material of commerce prescribe maximum 
lor minimum /imits of the particle sizes. In addi- 
nn. a knowledge of the distribution of particle SIZES 


ntrol purposes and in 


arch. In order to test materials for conformity 
or to measure the parti le size 


used by many indus 


the buver and seller 


pancy in some cases is due to inadequate sampling 
rrocedures, but when good representative samples, 


used, differences fre 


operation consists of agitating the 


openings have passed. Tl 
carried out ina satisfactorily reproducible manner 


1] The weigh 
le errors On 


f the largest Sol 


ng tolerance in the averag 


iterial on a sieve of known opening until substan 
ily all the particles that are small enough to pass 


ith shaking process can 


» that the error from this source can be made satis 


ts of the fractions re 


ned by the various sieves can be determined with 


the other hand, the 


tluation of the size of the openings of testing sieves 


irces of error in sieve 


nalysis. Consideration of this error is the subject of 


2] for testing sieves allow a manufac 


re size of the opening 


7 percent, depending upon the size of 


These tolerance 


v if only a roug 


or if it is desired 


s are rather wide but 
h idea of the particle 
that a certain material 


finer or coarser than some stated size. If an aceu 


ile parts le size distribution is required, especially 


th material with a narrow range ol partich SIZCS, 


then these tolerances are so | 


eves almost 


useless without 


arge as to render the 
an individual calibra 


on to determine the effective opening of each sieve 


both the reproducibility and accuracy of the 


ye reent.,. 


1.1. Effective Sieve Opening 


The difficulties encountered in the evaluation of 
sieve openings arise from three sources, none of which 
appear to have received sufficient attention, First, 
all the openings of a sieve are rarely of the same size 
Rather, the SIZCS of Openinys are distributed neceord 
ing to some probability law. The over-sized particles 
can pass the larger holes and, thus, the openings that 
are effective are somewhat larger than the average 
opening. Second, the separation achieved by a sieve 
is not complete, but instead seme under-sized parti 
cles alway s remain on the top of the sieve The sieve 
opening ts thus effectively somewhat smaller than the 
average opening. Third, the opening is effective in 
three dimensions, and the plane defined by the sieve 
cloth may not coincide with the plane defined by tht 
effective opening The effective opening W ill thus be 
larger than the average opening, which is measured 
yy the projection on the plane of the sieve cloth 

The first and second of these phenomena are acting 
Ih opposite directions, and the effeets of one might 
nullify the effects of the other llowever, the man 
ner and speed of shaking affect both: therefore. the 
effective opening Is not a constant for a given sieve 
but depends also upon the way the sieve is used 
The general dependence of the effective sieve open 
ing on time of shaking ts tlustrated in figure 1 
Assume that a sieve with an opening size distribution 
given in figure | A, is to be used to separate ito 
two fractions a material with the parti le size cis 
tribution given in figure 1, B After a very short 
time of shaking, the particle size distributions in the 
fractions passing and retamed are as illustrated in 
figure 1, C The separation is very incomplete with 
large amounts of under-sized particles retamed After 
a longer shaking time the distributions given in figure 
l 1!) any be tie hieved and the ette tive ize ol the 
sieve may be defined as the pom at which the curve 
cross, Which is also the pom at which the curves are 
the steepest Figure 1, KE, illustrates the effeet of 


an infimitely long shaking time As the shaking 
progresses it Is seen that the effective opening be 
comes larger It becomes evident that the effective 
opening depends upon the time of shaking In 
addition, the type ol shaking motion would have a 
similar effect Any method of calibration for effee- 
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tive opening must necessarily take these two effects 
into account 

If all the openings of the sieve were exactly the 
sume size, the shaking time required to reach a con- 
dition similar to figure 1, D, would be considerably 
shortened. It is to be noted, however, that even in 
perfectly uniform sieves the effective opening is not 
equal to the average opening unless the shaking is 
continued for an infinite time 

The usual method of calibrating sieve openings 
|3] consists in making measurements on a magnified 
projection of the sieve cloth. Observations by this 
method are rapid, and the openings may be inspected 
for uniformity during these measurements. This 
method, however, has two serious disadvantages 
First, the average opening is measured, taking no 
account of the difference between average and effec- 
tive opening. Second, the calibration requires 
special equipment and skilled operators. Weber 
and Moran [4] suggested that the projection method 
be used to measure a large number of openings, and 
the effective opening be determined by use of an 
empirical relation between the statistical parameters 
and the effective opening 


This method requires | by the authors in a previous paper [1]. 


the same special equipment and even more skilled 
operators. a a 
Fagerholt [5] showed that the effective size of the 
sieve after a shaking period of time, f, is equal to the 
average diameter of particles that passed by con- 
tinued shaking under the same conditions during 
the interval of time ¢ to 3¢. This conclusion is based 
op assumptions that are not strictly true. For 
practical purposes this method is much too involved 
because it requires an independent particle-size 
determination for each sieve for each sieve analysis 
It offers a way, however, for determining the effective 
size of the sieve for very irregularly shaped particles 


2. New Method of Calibrating Testing Sieves 


The authors in a previous paper [1] pointed out 
that the effective opening can be determined directly 
and simply by measuring the sizes of spheres that 
will just pass. In that paper exploratory work was 
described on the use of a calibrated sample of glass 
spheres for measuring the effective opening of testing 
sieves 


2.1. Particle Shapes 


It is recognized that the sieve openings are square 
or slightly rectangular in shape and that irregularly 
shaped particles can pass through even though one 
of the dimensions of the particle, or “an average” 
of all dimensions, is considerably larger than the 
diameter of the opening. This is especially true for 
needle-like shapes. The average diameter of irreg 
ular particles that pass a sieve cannot be considered 
equal to the diameter of spheres that pass the sieve 
The glass spheres used in this work were only in 
tended to measure the opening of the sieve. The 
calculation of some “average” diameter of particles, 
which deviate from a spherical shape, is a separate 
problem that introduces factors not directly related 
to the methods for evaluating effective sieve openings 


2.2. Rectangular Openings 


If the openings are rectangular rather than square, 
then the “effective” size for irregularly shaped par- 
ticles will be increased, but the effective size for 
spheres will be the same as for a square with the 
same dimension as the minimum side. For this rea- 
son, sieves that are to be used with needle-like par- 
ticles must be examined independently for squareness 
of holes. Fortunately, the openings of most testing 
sieves made in recent years are essentially square in 
shape. 


2.3. Uniformity of Openings 


The method of calibrating testing sieves by use of 
glass spheres takes into account the effect of non- 
uniformity, but it does not measure directly the uni- 
formity of the openings. If there is reason to suspect 
that the sieve openings are not sufficiently uniform, 
they should be checked independently either by the 
usual projection method or by the method proposed 
The latter 
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based on the effect of uniformity of openings on 


sieving rate. The amount of nonuniformity that 


in be tolerated in testing sieves was considered in 


re 


pl 


Sieves that do not appeal obviously deformed 

isually sufficiently uniform so that the 
ere method of calibration will correct for the 
Hl amount of nonuniformity that does exist 


3. Selection of Glass Spheres 


lhe spherical shape was chosen primarily to elim 


ite anv doubt or question about “average” diam 


or the orientation of the particle in passing the 
ning. The spheres were not intended to be repre 
tutive of the partr les to be used on the sieves 
is their diameter to be considered equal to the 
ve diameter of irreg ileus part les that might be 
d on the same sieve 


Glass was chosen because its hardness minimizes 


wh an 


paratus 


yr 


and because glass ure readily 
ible Glass beads of the type used for projec 
Gooden and | pdike 7} is 
Class beads 
markings * are a suit 
source ol class spheres for sizes [rom SO to 1,000 


es No. 170 through 18 


of glass spheres nie 


ison spheres ol 
screens were used by 
irce [or spheres of 100-u diameter 


tive type used for highway 


Larger and smatler 
being made for othe pul 


es, and a wider range of sieve sizes could be 


red if desired 


3.1. Separation of Spherical Particles 


Vhre rla = beads for h 
not ail spherical nor are the 


ssurv to separate the 


way markings as purchased 
Sizes cCOntIMUOUS 
vas mer spherical particles 
dumbell 


The spherical shapes were 


jageed 
‘ effectively sep 
ed from the other shapes by rolling the particles 
incline The whereas 
pheres either roll in circles or not at all. The 
firvure 2 It 
disk 
} rpm and tilted slightly 
sieve fraction can be 


if elliptical tear arop and 


spheres roll true 
shown in 
entially of a smooth flat 

eter, turning at about 


used Is consists 


about 20 in. in 


iv ome separated at a time 


the larger sizes (about 1,000 y) the inclination 
the disk is only a few degrees from the _ hori- 
mtal, but for the finer sizes (100) inclinations up 
10° from the horizontal are required It is desu 


ble that the particles roll over as many 


1 
! 


times as 
ssible: therefore. a disk 30 or 40 in. in diametet 
ild be desirable for partie les larger than 1,000 u 
th larger 


1MLOW 


spheres i disk made of ordinary 
glass is an acceptabl 


No difficulty 


shape as it 


material for the sep- 
ition Is experienced in cutting it 
need hot be perfect 
id to the turntable by means 
the smallet 
in that it 
A polish 
howevel 
plate 


thr desiree 
na It ear 
nu very 


t 
, : 
icles 


be faster 
viscous stopeock grease lor 
class has a serious disadvantage 
ndly builds up static electric charges 

per plate is generally more satisfactory 


in these experiments a 


polished copper 


14) 





rade I 


one-sixteenth om thick Vil Iscal 


on particle 
than about 500 yg 
fed onto the 


The part cl were rotating disk Ih a 
thin stream For the finest es the number rate 
of feed was quite large, and omm ola sLopcor hk 
made a suitable ruts Kor the larger vA only 


i 


three or four prece ata time could be fed without 


interference \ ~ 
and mstead a cardboard tube with 


touch thr 


stopeoe K 


THT titbstiite at this 


low rats a notch 


none side was supported to just urlace 


of the plate The motion of the plate beneath the 
tube caused the particle to roll out of the noteh in 
the tube one laver high and three or four wide 
depending on the width ana heimoht of the m teh 


A gentle blast of air was directed into the stream of 
particles immediately after they 
plate that stuck 
spherical piece and started all particle 

Phe spherical particl immediately 
plate into a suitable hopper while the nonspherical 
particles either did not roll at all or rolled 
slowly or in curved paths and were carried by the 
and fell into the cdiseard hopper \ 
strong blast of air was used to dislodge the particle : 
that did not roll igure >; shows the 
particles, and figure 4 the discarded particl 


10 Ib of 74 to 149 uw (No. 100 to No. 200 


were fed onto thre 


This dislodged any behind a non 


~ rolling 
rolled off the 


more 


turning cdisi 
pheri al 
About 


seVes 


material wa separa fin 20 hours on the apparatus 
deseribed About y percent ol thi particular 
material was discarded In some lot particularly 
in the larger size 0) percent was discarded Ihe 


large 312 separated more 
Nleasurements of = the 

that on the 

projections of the spheres differed bv more than 2 or 

in 4 percent of the bead 


quick ly 
howe al 
ot the 


sphere diameter 


average thie miayo! and minor axe 


} percent only Oniv one 


Ficure 3 Spherical particles separated from mixture 


bead in 14,000 was observed whose major axis was 
more than 1.5 times the minor axis. At least 90 
percent of the beads could be considered spherical 
within the accuracy of the measurements. 


3.2. Continuous Size Distribution 


In order to obtain spheres of a continuous size 
distribution the material selected as spheres was 
carefully sieved into the closest sieve fractions 
(,2 series), and then equal weights of each fraction 
used in the final mixture. This procedure 
produces a logarithmic particle size distribution. 
This distribution was chosen because the precision 


of the measurement of the sieve openings is the 


were 


same for all sizes 

The details of the manufacturing process for these 
glass beads were not known. Various lots as re- 
ceived from the manufacturers were found to contain 
a discontinuous distribution of particle sizes. It 
is not known whether this resulted from the manu- 
facturing from a subsequent particle 
size separation. For example, one lot of spheres 
was examined in which the 350- to 500-u sieve 
fraction (No. 45 to No. 35) consisted entirely of 
spheres ranging from 350 to 400 » and from 440 to 
500 uw. There were no spheres at all in the size 
range 400 to 440 uw. The No. 40 sieve has a nominal 
opening of 420 uw, and this lot consisted of equal 
weights of No. 35 to No. 40 fraction and No. 40 to 
No. 45 fraction. The sieves indicated a smooth 
continuous particle size distribution, but closer 
measurement revealed the discontinuity. For this 
reason it was necessary to obtain spheres from 
several manufacturers in different lots and to mix 
them so as to obtain a fairly continuous particle 


process or 


Figure 4 Discard particles separated from mixture 


size distribution. The continuity was checked by 
withdrawing a sample and applying a condensed 
version of the calibration procedure to be explained 
in a later paragraph. The use of a mixture of many 
different lots of spheres produces a calibration curve 
with many minor fluctuations 


4. Preparation of Samples 


A convenient total weight of spheres in the final 
calibrated sample is about 100 ¢ The total number 
of samples prepared in lot 3 was 256 (equal to 2° 
which required about 56 lb of glass spheres having a 
continuous range of sizes in the proper proportions 
After lot 3 had been assembled and the preliminary 
check calibration indicated that the particle size 
distribution would be satisfactory, the individual 
samples were prepared. The “Boerner Sampler’ 
16] was used as the sample reduction device! It is a 
riflle-type apparatus with a funnel, gate, and chute 
attached above the compartments, which are ar- 
ranged circularly. The combination of gate and 
chute makes it impossible to feed too fast or to feed 


at the wrong angle. It has been found to give at 
least as good results as any known sample-reducing 
device or technique. As an additional precaution 


to neutralize any errors of the sampling device, care- 
ful attention was paid to the side the sampler from 
which the sample emerged. Although 256 samples 
were desired, 512 samples were first made; these were 
recombined in pairs so that the final sample rep- 
resented equal quantities of material from each 
side of the sampler. Figure 5 illustrates the process 
for only four samples. 

§ Acknowledgment is gratefully made to Lawrence Zeleny, Grain Branch, Pro 


duction and Marketing Administration, U. 8. Department of Agriculture, for 
the loan of these samplers 
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4.1. Selection of Samples for Calibration 


lt was recognized that the essential problem in 
iibration The spheres 
were divided into 256 samples by the best 
ivailable, but it must be admitted that all 256 sam 
ples may not be identical. Every twenty-fifth 
ample in the orderly arrangement of numbers given 


was adequate sampling 
means 


to the samples was chosen as a sample of the 256 
Kach of these 10 samples contain 
approximately 10 million spheres, therefore, they 
must be further sampled to obtain a practical num 
Because spheres of 


making 10 im all 


ber small enough to mensure 
ipproximately the same size do not tend to segregate 
md present an problem of measurement 
each of the 10 samples was separated into 14 sieve 


eCusicl 


fractions 
The spheres for the actual 
mounted on a gelatine-coated microscope slide In 


measurement were 
to make sure that an adequate sample was 
obtained from each sieve fraction, four different slides 
were prepared, each slide representing particles from 
i different portion of the sample container. ‘Twenty 
ve particles chosen at random from each slide were 
A total of 1,400 particles in each sample 
14.000 


ordet 


ait asured 


was measured: altogether particles were 


measured for the 10 samples 


4.2. Measurement of Diameters 


The apparatus for the measurement consisted of a 
projection microscope with its accompanying light 
source and screen. The measurements were made 
on the projection of the image of the particl s on the 
The magnifications and projection distances 
were arranged so that the projections of the particles 
between 6 and 30 em in diameter. Two 


screen 


were 
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different 


32 mm, according to the magnification desired The 
ocular was 7 


Inicroscope objectives were used, 16 and 


5 power with a seale mounted inside 
The measurements were made in air A carbon are 
and condense! lens SVsterm) were used for illumination 
The magnifications were checked at frequent inter 
vals throughout the measurements bv the use of two 
different micrometers Both of these stage, 
micrometers have been calibrated at this Bureau 
and the observations obtained are accurate to within 
1 wu 

In order to speed the Pensurement 


stage 


a seale was 
drawn on a separate plece of paper that could be 
moved about the screen The distance between the 
microscope and the image was adjusted so that the 
vraduations on the hand seale exactly coincided 
with the projected scale in the eve plece The hand 
scale could then be moved about, and several beads 
hear the center of the se reen could be measured very 
rapidly. ‘Two people working 


together measured 
25 particles in about 2 minutes. Counting the time 
required for changing slides, preparing samples, and 
relaxation to prevent evestrain, 300 or 400 particles 
When nonspherical parti 
encountered, the minimum diameter 
This was done because it was 
desired to know the size of the hole through which 
i particle Pusses rather than the ave rage size of the 
particle 

The individual diameters were measured to the 
2 to 5 percent, depending upon the part of 
\lore precise PenusuretMents could 
since sampling rather 
factor that limits the 
precision Was unnecessary and 


were measured per hour 
cles were 
Was measured 


hearest 
the scale in use 
have been made. However 
than 
Hecuracy 


measurement is the 
vreatel 
would have greatly lengthened the subsequent work 
ol computation (As it was, SS points were obtained 
in the range 60 to 1200 yu 
within the range 


averaging about five points 
between eu hy pull ol steves of the 
'y 2 series 

5. Computation 


The measurements of the diameter of each spl re 
were recorded in terms of the projected seale reading 
Since it was desired to caleulate the weight fraction 
finer than each size, the first step was to determing 
the diameter frequency distribution within the 100 
spheres representing one sieve fraction of one sample 
The actual diameter for each size was evaluated by 
multiplying the 
factor These 
distribution 


scale reading by the magnification 
computations give a number-size 
On the assumption that all the parti 
cles were of the same densitv and same shape the 
weight Is proportional to the number of particles and 
to their diameter cubed The weight fraction of 
spheres of diameter d, in the sieve fraction thus 
becomes nd? /Snd If wis the weight of the sieve 
fraction and Sw the total weight of the sample, then 
w/> wis the weight fraction of the sieve fraction in the 
total sample. The weight fraction (/ 
diameter d in the total sample ts 


ol spheres ol 


nd , Ww 


rnd lw 
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5.1. Sample Computation 
Typical data are given in table 1 for one sieve 


fraction of one sample. Each diameter recorded in 
table 1 is the mean of the size range 
measurements between 1.875 and 1.9 
as 1.90. The frequency distribution, diameters, and 
the computations leading to the weight fractions are 


for example, 
25 ure recorded 


viven in table 2 

The differential type of particle size distribution 
(expressed by weight fraction per unit size range as 
a function of size) is obtained by dividing each / by 
its corresponding size range. The integral type of 
particle size distribution (expressed by weight frac- 
tion finer as a function of size sum- 
ming the f’s. The particle size associated with the 
sum of the fractions finer is the smaller limiting 
diameter of the size range. This is in contrast to 
the mean of the size range used in the above expres- 
sion for f and in the differential type of particle size 
distribution 

These computations were repeated for each fraction 
of the ten samples, and the values were averaged to 
obtain the final calibration. The calibration is given 
graphically in figures 6 and 7 


Is obtained by 


5.2. Evaluation of Sieve Opening 


the opening with a 
the entire sample is placed 
on the sieve or sieves in question. The sieves are 
shaken with the calibrated spheres in the 
manner as will be used with the unknown material 
The spheres are then carefully brushed from each 
sieve and the weight fraction finer than each sieve 
evaluated The effective size of the sieve is then 
read directly from the calibration curve (fig. 7) 
Care must be exercised to avoid loss of the spheres 
in order that the sample can be used again. Expe- 
ience has shown that the fines are most easily lost; 
therefore, if the total weight of the sample decreases 
slightly with repeated use, it may be assumed that 
the loss is in the finest sizes. The particular oper- 
ation in which there is most chance of losing par- 
ticles is the transfer from the sieve to the weighing 
container The loss of particles in this operation 
can be minimized by inverting the sieve into a deep 
funnel and removing the contents with a stiff brush. 
The stem of the funnel should fit snugly into the 


recelver 


In order to evaluate Ssleve 


sample of the spheres, 


same 
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5.3. Application to Sieve Analysis 


One of the most disturbing featuresof sieve analyses 
is the inability to obtain the same results with the 
same sample using different sets of sieves that have 
been considered identical. The small variations in 
the sieve openings inherent in the manufacture of 
the sieves are the cause of this trouble. In order to 
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when the samples are used Third, a calibration 
error occurs in the sampling and measuring involved 
in the calibration process. The magnitude of each 


each of these errors has been determined and. also 


4 
- their effect on the over-all reproducibility. and ae 
curacy of the method 
The dividing error can be considered in two parts 
c One is due to the inaccuracies of the sample-reduces 
7 
/ and the other arises from the use of only a relatively 


limited number of partich It is possible to esti 


/ 
mate the error of the latter by the application of 
If } is the number of particles ot one 


: Statistics 
f e remaining in the final sample after three or more 
3 eparations on the sample divider, then according 
to statistical theory appendix | the tandara 
f deviation associated with this number ts approxi 
4 ; mately 4 The error in the determination of size 
> ol particle from this cause can be readily evaluated 
appendix 2) and listed in table 
Ss. ( } j 
( j 
‘ 
. Beeau ol thy ood construction of the sampling 
levice and the additional precautions that) were 
" observed. tt rensonable to assume that the error 
nl 1 frome the naccurack of the amp ner aevics 
mall Indeed, it ean be hown that this error 4 
thy results of an sieve analvsis to that 20 small that tf exert ho appre ible unfluenee on 
would be obtained tf sieves with nominal the final sum of all error appendix 3 
d openings Were used, 1 isonly necessary to piol The error due to stevin depend pon thie pra 
esults as curmulative percentage finer (or coarse! ticular stevin method emplo ad In order to obtain 
n tunetion of t he calibrated Opening ol thr weve contnve tien of le eTrotl t hist nil hit bog encountered 
ch a plot for a bone char (char 32) is shown th | one sample was sieved nine tiny imtng the particu 
i! ~ Ir) my this eurve thr cumulative pereenht he procedure modu thi lnborators Ihe tandard 
finer (or coarser) than the nominal openings can | deviations representing the variations among thes 
obtained at points indicated by arrow hemes nine sheVvil iil iso fist dl in table ait ‘ 
corrected sieve analvsis mav be caleulated Varintions are due to steving alone because the same 
this Wav sieve anal es made at different lab ample tha nin weve andl the aime calibration 
orl with diff rent ol sieves may bye reported eurve were Uu i] 
The variation amor 


the calibrations of the ten 


reference to the same sieve Openings 
different iLEnh pp! l a i re ot the th of tha 
6. Evaluation of Errors dividing and calibration erro Phe dividing error 
bas alread been detern ed and can be ubtracted 
Ih errors inher nt oun th method OF mensurin fron thre ] to Vield the calibration error alone 
effective size of the openings of testing sieves by appendix 4 Chis error is also listed in tab 
ans of calibrated iss pheres ari from three 
ources lirst 2 d error occurs When the ' 
wn ples are pre pared Second, a 8 nd error occurs 
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It is noted that each of the errors of the same 
order of magnitude, indicating that 
precautions in the sample dividing or calibration 
would have been wasted effort unless some method 
of sieving with less error could have been found 
The and reproducibility the glass 
sphere method of calibrating testing sieves are tab- 
ulated in table 4 The values listed are the maxi- 
mum variations that might reasonably be expected 


ure 


any additonal 


accuracy of 


and are taken as three times the standard deviations 
The accuracy tha of the eal 
bration and dividing errors and is an indication of 
the precision with which the size of the spheres ts 


is measured by sum 


known 

There are two types of reproducibility. The varia- 
tion to be expected when the same sieve is calibrated 
several times with the same sample of glass spheres 
is that due to the sieving error. If different samples 
of glass spheres are used, then the dividing error is 
included also 

As a general summary in regard to errors it can be 
stated that both the accuracy and reproducibility 
of a glass sphere calibration is about | percent of 
the size of the sieve opening. This is significantly 
better than the 2- to 7-percent variation allowed in 
the present specifications for testing Even 
more Important, however, the glass spheres measure 
the effective opening (for particles not too far from 


spherical) rather than 


Sle Ves 


the average opening 
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8. Appendix 1 
Acknowledgment is made to J M. Cameron and Joh; 
Mandel of the Bureau for the derivation leading to the follow 
ing resuit 
\ large number of particles i plit into two groups as ever 
as possible Kach of these groups is again divided into two 
groups as evenly as possible, and so or After / ich split 
ting operatio 2* subgroups are formed 
If there were \ particles in the original ag 
would expect each subgroup to have V2 
Sut, due to randomne in the partitioning the 





ber of particles observed will vary about this expected number 
As few as zero could be found or as many as \ 

Assuming that each particle has an equal chance of 
into either group at each partition, tl probability that 
particle be found in a particular ibgroup is 1/2*. and the 
numbers of particle inn the ibgrouy will be binomina 
distributed From the mathematic of the binomial d 
tribution, the standard deviation of the inber of particle 
found 1 tM iow subgroups 

l l l 
?’ V Vai (1 -gi day" ~ae) 
It we bn ited that for values of / reater than 5 ( 
of the factor i parenthesis practica vie henes 
tandard deviation becomes 
a vr for 4A 3 
9. Appendix 2 

It veight nit the standard deviation , np 
where p the particle densits 2 6 Cl for ub +} 
shape factor 6 for sphere rnd dois the particle a 
meter If the weight of a fractio of particles 
Sarnie mc then the mnber of particl mbtained | 

. 
pad 
The standard ak . weight eens 
\ pad \ ) 
p 
Ihe Ve frac ) obtains | divid | ota 
we nt fa he tractio 
a pad 
v3e Vz 
Sinee the error is small, to a first approximation he rat 
error in weight tractior imit to the error u 17 its is equ 
to the slope of the calibration curve at that point Ihat is 
o,/a,~dfids where df/da is the ope of integral type particle 
size distribution curve (fig. 7) and is equal to the differ 
ential particle size distribution given in figure 6. The stand 
ard deviation in size units can thus be expressed in terms ¢ 
easily measurable quantities as 
pacl® /df 
Cs 
V xw / ds 
his quantity is tabulated in table 3 for each of the sieve 


fractions 











the ample 





" Column 5 ineluce rror due to oe it i“ i! I 
10. Appendix 3. Humn 5 includes ¢ Mccurat 
livider, while column 4 does not 1} licate therefore 
that tl rror ad to tl amp divider oO sta nu to bn 
e course of the calibration procedure he 10 sample x . a ‘ pou : 
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Precise Topography of Optical Surfaces’ 


James B. Saunders 





Phe unit of ke ial ised in measuring optical surface features, by interference of 
u one-half e wavele of the monochromatic light that is used \ method is 
described in wl th the of length is a much smaller fraction of the wavelength Popo 
rraphic maps of optical surfaces are made in which the contour interval is less than one- 
thirtieth the wavelengt! rhe smaliness of the unit depends upon the quality of the surfaces 
Irregularities, too small to be detected with the Fizeau or two-beam fringes, are made to 
appear very prominent Surface markings, caused by the final polishing actions, can be 
made to stand out d relief 








The accuracy of results that may be obtained with 
any type of interferometer depends upon the accuracy 
with which one can locate a line in the fringe field, 
along which the order of interference is constant. 
The intensity of light in a direction normal to these 
lines, in a set of two-beam or Fizeau ? type of in- 
terference fringes, is a smooth mathematical function 
without abrupt changes, whereas for multiple beam or 
Fabry-Perot * type fringes it has sharp maxima and 
minima for transmitted and reflected light, respec- 
tively. The accuracy attainable with the two types 
of fringes is roughly proportional to the ratio of 
fringe separation to fringe width, the fringe width 
being defined here as the width at one-half maximum 
intensity 

Accurate contour maps of one surface relative to 
a standard surface of known shape are obtained by 
photographing the interference fringes, produced by 
these two surfaces, with monochromatic light. The 
contour intervals are one-half wavelength. As in 
geodesy, the smaller the contour interval the more 
detail is revealed in the contour map and the nar- 
rower the lines the more detail is revealed along the 
line. With broad fringes, of the two-beam type, 
irregularities of several hundredths of a wavelength 
can be completely hidden in or between the dark 
broad lines. An irregularity of several tenths of a 
wavelength can be completely invisible if it falls 
within the area between two adjacent Fabry-Perot 
fringes that appears uniformly illuminated. The 
one-half wavelength interval is, therefore, inadequate 
for representing local irrregularities on high-quality 
optical surfaces, 

This paper describes a procedure whereby very 
narrow fringes are used for producing contour maps of 
optical surfaces with contour intervals of a small 
fraction of a wavelength. These intervals can be 
made as small as the narrowness of the fringes will 
permit If m resolved fringes can be equally spaced 


This work was done as part of a research project sponsored by the Air 
Fores he paper was presented at the March 1950 meeting of the Optical 
Society of America (Abstract, J. Opt. Soc. Am. 40, 258 (Apr. 1950 Since then 
the author's attention has been called toa paper by Takemara Sakurai and Kéré 


on “Test of an 
and published in the Scientific Report of 


Shishido at the Research Institute of Scientific Measurements 


optical flat by the Fabry-Perot etalon 


the Research Institute, Téhoku University, [A], I, No. 1 (May 1949 
eau. Compt. rend. 44, 1237 (1882 
C. Fabry and A. Perot, Ann. Chim. Phy 12, 459 (1807 
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between two of the original fringes, the separation 
of which corresponds to one-half wavelength or to one 
order of interference, then the contour intervals 
will be [1/2(m+1)]A. Values of m as large as 19 
have been used, for which the intervals are 4/40 

The usual procedure in precision measurements on 
glass optical surfaces is to place the unknown close 
to and almost parallel to a standard surface of 
known shape, with a thin air wedge between, and to 
observe the fringes produced by the interference of 
monochromatic light reflected normally from these 
two surfaces. Figure 1 shows a set of two-beam 
interference fringes, observed by reflection with a 
Fizeau viewing instrument, which one might have 
under observation for the purpose of measuring the 
deviation of an unknown surface from a standard 
plane surface. The orders of interference corre- 
sponding to each fringe, although not usually known, 
are indicated by the numerals. The exact position 
of the lines that represent the loci of integral orders 
of interference, are only known to lie somewhere 
within the darker portion of the fringes. This locus 
for fringe number 3 is indicated by the narrow solid 
line. Usually, as will be shown later, these lines are 
not as smooth * as the much broader fringes appear 
to be. The assumption that this line hes at the 
center of the dark fringe is usually false and, conse- 
quently, results based on this assumption are in 
error by the amount of displacement of the line from 
the center: which can amount to over a tenth of the 
separation between adjacent fringes. Even if one 
could locate the lines of integral orders, he would 
still be unable to measure the change in order of 
interference along a straight line such as the one 
shown dotted in figure 1, for the purpose of surface 
measurement. The change in order of interference 
along this straight line is an accurate measure of the 
deviation of the unknown surface from a plane, but 
if the shape of the surface between the fringes is 
unknown one has no way of determining the frac- 
tional orders of interference at points on this straight 
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2 ire a good indication 
the location of lines of con- 


n fr 


et 


the edges of the 
In f 


4 The irregular serrat long nges 


of the shape of the lines of constant orders 


stant order (although not integral) are most accurately locatable along the 
margin of the fringes where the density gradients are a maximum This is 
the most accurate means, known to the author for locating lines of integral 


differences in order of interference for surface measurements 
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line that are not also on the lines of integral orders 
Linear interpolation, for points not on a fringe, can 
be in error by several tenths of an order. The 
measurement of deviations of a surface from a plane, 
by measuring the change in order of interference 
along this line, is in error both from the inability to 
locate the lines of integral orders of interference and 
also from interpolation based on the false assump- 
tion that the order of interference varies linearly, 
with the displacement of the point from the adjacent 
fringes. Unless a fringe passes through the center 
of a surface, the deviation of the surface at this 
point, relative to cannot be accurately 
measured. 

The accuracy of measurements on the deviations 
of the unknown surface from an arbitrarily chosen 
reference plane, along any one of the lines of integral 
orders such as the narrow solid line in fringe 3, is 
hmited only by the accuracy with which the line 
can be located. The deviation of the surface at any 
point on this line, from a plane through the indicated 
straight line (which is conincident with neither of 
the reflecting surfaces but intersecting the unknown 
surface at the two ends of fringe 3) and making an 
angle @ with the standard reference plane, is equal 
to the product of @ and the perpendicular distance 
of the point to the straight line where @ equals one- 
half the wavelength of the light used divided by the 
fringe separation The fringe separation is defined 
as the average distance between the centers of acdja- 
cent fringes and is the separation that would be 
obtained if the unknown was replaced by a true 
plane, placed in the position of the arbitrarily chosen 
reference plane. The value for @ should be constant 
for all computations, although the observed fringes 
may neither be parallel nor equally spaced. The 
use of different values for fringe separation is equiva- 


others, 


lent to measuring deviations from different planes 
at various angles to the standard 

If both the standard and the unknown are coated 
with a uniformly thin coat of highly reflecting metal 
such as silver or aluminum, the fringes in which the 
lines of integral orders of interference are located 
become narrowel! The narrowness of these tringes 
relative to their separation, increases with the reflect 
ance of the surfaces and also with decreasing angle 
of incidence. The decrease in width of these narrow 
line fringes, relative to their separation, increases thi 
locating the lines of integral order of 
interference With modern methods of metal depo 
sition and with that are almost 
parallel to each other, fringes may be obtained whose 


accuracy of 


sIpoot h surfaces 


width, at one-half the maximum intensity, are only 
one-fiftieth part of their separation With such 
narrow tringes., deviations can be measured to better 


than one-hundredth of a wavelength. These narrow- 
line or Fabry-Perot type fringes are dark by reflected 
light and bright with a dark background, with trans- 
mitted light Figure 2, A, shows a set of two-beam 
fringes produced by reflected light, and B shows a 
graph of multiple reflection fringes, formed by the 


same surfaces by transmitted light after applying 
the thin coat of aluminum Note the increased 
accentuation of detail along these narrow fringes 


over that shown by the broad ones. 

When measuring features to one-hun- 
dredth of a fringe, the mechanical errors in the micro- 
meter (assuming a micrometer eyepiece to be used) 
and personal errors of settings become large unless 
the changes in order of interference with movement 
of the cross hair, are relatively small. This requires 
that the fringes have considerable separations. At 
least two fringes must appear on the surface in order 
to determine the unit of deviation of a fringe from 
a straight line corresponding to a known unit of 
deviation of the surface from a plane. This unit is 
usually one-half the wavelength of the monochro- 
matic light being used. Figure 3, A, shows such a 
set of fringes, formed by the same surface that was 
used for figure 2, B. Further enhancement of detail 
along the fringes in figure 3, A, over those of figure 
2, B,isobvious. This is partially due to the decrease 
in width of the fringes relative to their separation, 
caused by making the wedge between the plates 
smaller and, consequently, reducing the angle of 
incidence that favors the narrowing of these fringes 


surface 


In order to further enhance irregularities along a 
fringe and still be able to measure the magnitude of 
the corresponding deviations in the surface, a mul 
tiple spectral line source, such as the mercury are 
Each spectral line pro- 

If the absolute orders 


for instance, may be used 
duces its own set of fringes 
of interference of and all of the fringes are 
known, the separation of the surfaces along the 
fringes is known and, consequently, the change in 
separation of the between two adjacent 
fringes of different known, although 
this difference in separation be only a small 


each 


surfaces 
is also 


may 


color s 


8. Tolansky, Multiple-bea nterferometry of surfaces and filr Oxfor 
the Clarendon Pre 104s 
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Fiaure 2 Equal-t! 


A, Two bea ! fl i 
Fiaure 3 Vu tiple-beam interference fringes. 
A, Singh wetral line Note increase hancement of surfac tetail B, Multiple spectral line One order of each of three lines shown on plat« W edge 
too stall t ide other line Note further enhancement of surface detail 



































on of the unit of length, one-half wavelength 
thickme ss and angle ol the nil wedge between the 
fine may that ol 
more closely spaced fringes of different color pass 


‘ now be adjusted x0 on two 


ough the center of the surface and one or more of 


| others just inside the eda The surfaces are 
us heat parallel “us is prac ticable thus favoring 

| attainment of the narrowest fringe possible 
itive to fringe separation, and measurements 


| racial distances Still greatet 


i\ made at all 
nhancement of the irregularities is thereby attained 
available for 
kk 


more fringes are quantita 
valuation of these irre B 
wes a set of fringes, produced by a mercury source 
which the order of 
ys the plate three-tenths 
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interference 
of unit 


change inh 


than 


total 


l less it 


> () 97 y 
\laximum enhancement of surface irre cvularities is 
tained by adjusting th 
sm as is possible and the s paration near to some 
number of half wavelengths If hac 
reflecting surfaces (the existence of which 
bw the these 
other 
width 


surfaces to as near paral 


evral one 
0 plane 
doubts dl 


rallel 
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author could be placed 
which the for 

If the separation is now 
the of 


produces maximum 


to each Is condition 
ol 
ed 
the resultant 
ation of intensity with change in separation, then 


aused by irregularities of 


Irniannite 


ist to some value inh neiwhborhood 


interference 


all deviations in phase, « 
few angstroms, such as may be produced by 
ndividual particles of the polishing 
visible Figure 4 
actually both surfaces contribute 


iV a 


compound, 


come quite shows the above 
tioned surtaces 
jually) in which the parallelism and separation are 


t favorable the of small surface 


( 


to revelation 


tail. Note the fiber-like appearance in the darket 
ons where maximum enhancement of detail 
ecurs. Quantitative measurements of these finer 


rularities have not vet been attempted 

The manner of determining absolute orders of in 
may deseribed by of figure 
Phe large circularly enclosed fringe svstem represents 
et that might | The 
he fringes are indicated by R for red, Y for vellow, 


rierence be Deans 5 


”« under observation colors 


(; for erTreen B for blue, and \ for violet The wiaive- 
neths of the light, forming the tringes, are identified 
usually by the color of the corresponding fringes 
The small circularly enclosed fringe svstem repre 
ents the same fringes as those shown in the large 
circle. They are straight as if formed bv parallel 


urfaces and are inserted to show the relationship 
ween the fringes in the large circle and the lines, 


which correspond to fringes, in the chart above 
! we assume the surfaces that form this’ trings 
tem to be extended in the direction of their de- 


reasing separation to their line of intersection, and 
ne could observe the fringe system in this extended 


irea, one would find and could identify the zero 
order of interference for each set (color) of fringes, 
all of which coincide Extending from the zero 


ordet fringe, one would observe equally spaced red 
lringes (assuming plane surfaces), the order of which 


could be identified by counting from the zero order 
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vertical lines 


associated with the 
for all the visibl 
pond 
this 


Thus the number 
indicates the 
thre 


POstLLOns 


Interlerence 
source, and then 
ol the tringes 
rhe ighborhood of 


ordet ol 
colors Ith 
tha 
chart, on the 
plates that several fring: 
Ing 0.12 in 
This space corresponds 1¢ 


PostlLons eorre 
One 
.o-u Separation ol the 


sces On 


fallina space correspond 
f th 


ol onl 


to Inecrense separation oO plates 


et | change Oooo mn 


the ordet of interference for this red spectral linye 
If the maximum and minimum separations of thy 
plates are adjusted to 3.2 and 3.4 yu, respectively, the 
parte ular set of fringes that corr spond to the lines 
falling in this range on the chart can be seen on the 
plates We therefore have al fringe sVstem) in which 


the inerease In separation of the surtaces from one 
fringe of a viven ¢ olor to another fringe of a different 
color is a small fraction of the conventional 
"A, that is usually used for this type of measurement 


Furthermore, the width of the bright line fringes 
relative to separation of adjacent orders, ts mad 
small r by the two surtaces bemg mac hearet 


parallel than is practl al whi nh using monochromath 


light. The narrower the fringes, relative to separa 
tion of adjacent orders, the more concentrated a 
cluster or grouping of fringes can be used, without 
overlapping. Figure 3, B, shows a constellation o1 


group ol fringes whose s« parations corre spond to an 
tnerease In separation of the surfaces of 270 A from 
fringe | to 2, and \ from 

When an optical flat is placed on an unknown 
surface with only an air film separator, the 
of the top plate forces the air film to become rela 
tiverly thin 
lute contact 
of dust will 


os 


2e0 tos 
we ioht 
{ 
but the plates will not come into abso 
at any Usually 


prevent film from becoming 


point small particles 


the an ibs 
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} od i d h a nt alt pre i re ou vr 
thin as would be attained in their absence The | index of refraction of a gas varies linearly with the 
usual range in thickness, for circular surfaces of | density. For a constant temperature of 25° C 
t-in. diameter and free from dust, is from 8 to 94. | which was maintained during this experiment, the 


For larger surfaces this range may extend to 10g 
The range, 0 to 3.64, shown in figure 5 is not prac- 
tical but was used for explanatory purposes only. 
The practical range is from 6 to 94 and is shown th 
figure 6 for several different light sources. If one 
different he ria choose a 
most suitable given film-thick- 


has to sources 


constellation 


ACCESS 
for a 
hess range 

A more elegant topographic map of an optical 
surface, relative to a standard, can be obtained by 
using monochromatic light and photographing the 
image of a given fringe in several different positions 
on a single film; the shift in the fringe system being 
accomplished by changing the optical separation of 
the surfaces There are several different ways of 
changing the optical separation, only one of which 
will in detail, The interferometer 
plates are enclosed in an airtight chamber (see fig 
7) in which the pressure of the enclosed gas may be 


be described 


controlled at will. The collimator and collector 
lenses form the windows to the chamber. A gas 


line, connected to a manometer through a two-way 
permits either compression by a com- 
evacuation by a vacuum pump. A 
needle valve in the line, not shown, facilitates con- 
trol of the pressure in the chamber 

According to the law of Gladstone and Dale, the 


stop cock. 


pressor, or 


density Is directly proportional to pressure, to a 
close approximation. The relation between refra 
tive index and P(no—1)/760 
where nis the refractive index of the gas, P? is the 
pressure in millimeters of mercury, and no is the 
refractive index of the gas at 760 mm pressure and 
at 25° C 

If An is the change in n for a given change of AP 
in ?, it may be shown that the corresponding change 
in order of interference, A.V, at a point where the 
separation of the surfaces is L, is 


pressure, is n—l 


2LAn 2L (n.—1) 


n 760 af 


AN 


From this relation one may compute the separation, 
L, that is necessary in order to get a chosen change 
in order of interterence for any desired pressure 
range. For instance, suppose the available pressure 
range to be 1 atm or 760 mm of mercury and the 
desired change in order of interference to be one (the 
maximum practical change for this work), then 
AP =760, and AN=1. Consequently, 
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It, for instance, the light to be used is the green line 
for which “A 2.7 10 
then L equal O.%0 mm 


with this se paration, pressure range 


of mercury cmand 7 Ll is 


2.8% 10 Consequently 
and light source 
a fringe can be made to take on any position between 
the positions of this fringe and that of one of its ad 
neighbor at th 
If there 1 
and one Wishes to interpolate nine additional, equally 


neen! extremes ol thy ressure 
J | 


range more than one fringe on the plate 
spaced fringes, between each of the adjacent original 
fringes thus completely filling thre photographed fe ld 
with fringes whose s« parations are one-t nth that of 
the original set, then it will be necessary to change the 
pressure in steps of 76 mm for nine times, or a total 
of nine-tenths of an atmospher If the wedge be 

tween the surfaces is such as to give a maximum 
change in order of interference of one-half, then one 
half-atm change in pressure is sufficient to move the 
center of a ring 
If, for this wedge 
fringes on the plate, the necessary change in pressure 
will be nine-twe ntieths of an atmosplh re It is, of 
essential that this pressure within 
that available, which in this parti ular case is 0 to 
760 mm 


completely across the surface 


one Wishes to photograph 10 


COUPS¢ range bye 


Figure 8 represents a topographic map with con 
tour intervals of approximately 170 A, or one-thir- 
tieth of a wavelength of cadmium green light The 
separation of the two surfaces was 0.82 mm, and the 
change In separation, in terms of interference fringes, 
was approximately one-half fringe (1,300 A 

To illustrate the application of these smaller-than- 
one-half-wavelength-unit topographic maps to quan- 





titative determination of surface features 


measure 
ments have been made on the deviations of one of the 
surfaces used to produce firure } B relative to the 
other. For this discussion one surface will be as 
sumed a true plane (called the reference surface) and 
the other (called the unknown) to have all the urre 
ularities corresponding to the irregularities shown in 
this fring pattern Lines and points have been 
inserted in figure 3, B, and reproduced as figure 9 
The deviations of the unknown from a true refe1 
for several points along the lines AKB 
plotted rn 


enee plane 
and (FD) have been computed and are 
figure 10. The chosen reference plane is an imagi 
nary plane that passes through points A and B 
consequently including all points on the straight 
line AB) and making an angle @ with the reference 
surface approximating the average angle between 
the two reflecting surfaces. The unknown surface 
and the chosen reference plane are made as neal 
parallel to each other as the irregularities in the un 
known will permit. The value of @ is obtained from 
circle ABDC, which is concentric with respect to the 
center of the plates, the straight line AB, and the 
perpendicular projections of points ¢ and J) on line 
AB Circle ABDC must intersect at least two loca 
table lines, such as AFB and CFD 
which the order of interference is constant, and be 
tween which the change in separation of the reflect 
ing surfaces is known. The line ALB is the lower 
edge of fringe 2, and line CFD) is the lower edge of 
fringe 1. The change in separation of the plates 
between these two lines is approximately equal to 
that between the centers of the corresponding fringes 


along ea h of 
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of figure 9 The lower 





nd, as Was given previously, equal to 270 A The 
agnitude of @1s taken to be 2d ) ) whet ed 
quals the change in separation of the plates between 
nges | and 2, )) equals the normal distance of the 
from the straight line AB, and ) 

point J) from AB 


arbitrary and 


wit ¢ equals 


This 


In no Way 


he normal distance of 
hoice for @ is: entirely 
ffeets the computed shape of the unknown surface 
To correlate the data shown in figure 10 with the 
opographic map (fig. 9) from which it was derived 
he designated points in figure 9 are labeled likewise 
n the graph. Thus, the ordinate of point / in 
ure 10 is the deviation of the unknown at point / 
branch of solid curve 22, 
iwure LO, represents deviations along the right-hand 
portion of line 2 figure 9, and the upper branch, 2 
having the cireled pots corresponds to the left 
hand portion of this line 

These measurements are not limited to 
long a single fringe or line, but may be 
along any fringe that 

The deviation of the unknown from 
he chosen reference plane, at a point on any other 
ne, such as pom|l F on line 1, is equal to )0 Nd 
vhere Fe is the normal distance between point I 
ind line ALB, and WN is the difference in order of inter 
ference between F and any point on line 2. Thus, 
measurements may be made at all points along each 
ind every contour line shown 

The plot shown in figure 10 does not give any 
idea of the shape of the unknown surface 


points 
made at 
nv point appears on the 
photograph 


better 
than the contour 
| hae \ alt es are, how evel 


map from which it was derived 
better defined numerically 
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The failure of all computed points to fall on a single 
smooth curve indicates that the unknown surface 
is not a fiwure of rey olution about its center The 
appearance of more than two values of a given radial 
distance and on a single line Is due to local depres 
irregular shape of the 
ree of local irregularities 


sions or elevations The 

curves also indicates the deg 
in the surface. These local irregularities represent 
deviations approximating 50 A, in some cases 
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directly applicable to re 

Brenner on copper and nicke 

density reductior at the cathoce 
densities for a given distance fror 

mined and compared with the observations 


puted values is regarded as fairly satisfactory 


l. Introduction 


Neither the cognizance of the phenomenon repre- 
senting the subject matter of the present paper, nor 
the statement of the basic causes producing it is 
original with us. Some 12 years ago, W. Blum 
Chief of the Electrodeposition Section, National 
Bureau of Standards, described to us the cathode 
film phenomenon as a thin layer of moving material 
very close to the cathode surface and as_ being 
caused by the continued reduction of the densities 
due to electrodeposition. Since the phenomenon in 
question is also partly hydrodynamical in nature, 
he suggested that we should prepare a theoretical 
study with a view to determining the density dis- 
tributions in cathode films 

Blum’s first thought was that the film might con- 
stitute a seat of discontinuity in densities, since to 
the naked eye a line of demarcation between the 
contiguous layers of liquid near the cathode surface 
is distinctly visible The line of demarcation 
should be due to the different indices of refractions 
of the liquid lavers on the two sides of the line. Thus, 
near the cathode 
surface is in one or two layers having dissimilar 
densities with sharp discontinuities. But, despite 
the appearance of the region close to the cathode 
surface, the opmion was ventured on our part that 
the exsitence of the sharp discontinuities could not 
be shown on the basis of purely hydrodynamical 
reasoning. The basis for this opinion was the belief 
that the phenomenon in question was essentially 
similar to the type of flow encountered when a 
heated vertical plate in the presence of cold air 
sets up convection currents since experience im 
this particular case does not show any discontinuities 
in temperatures, it might be expected that density 
discontinuities would not be present in the cathode 
plate problem. Unfortunately, due to other pressing 
work at that time, it was impossible to proceed 
with a mathematical solution of this interesting 
problem 

About 3 vears ago, Blum repeated his request for 
an analytical study of the problem, since he was 
planning to resume the experimental investigation of 


according to these ideas the flow 


Hydrodynamics of Cathode Films 
Garbis H. Keulegan 


sults obtained by the freezing method are effected 
baths are considered 
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\ hvdrodynamical theory of the cathode diffusion laver based on the Assumption of 
a uniform density reduction at the cathode surface is presented By means of an approxi 
mate solution for the flow and the diffusion equations, the laver dimensions, the velocity 
distributions, and the densities are determined Further derivations from the solutions 


The data of 
Using the observed values of the 


urface, the maximum width of the laver, the averag 
n the cathode, and the mobilities of the cations are deter 


The agreement between observed and com 


cathode films. The computations that appear in 
this paper were completed shortly after the receipt 
of this request. Meanwhile two other papers on 
the subject, one by Agar [1]! and the other by 
Wagner [2], were published. Since the present 
work proceeds on a mathematical basis different 
from these two papers and work was 
conceived with a view of interpreting the results 
of A. Brenner, a member of the staff of the Bureau 
on cathode film densities, it is hoped that some interest 
in the results achieved will develop 

In the present treatment it is assumed that the 
reduction of densities over the vertical cathode 
surface is constant. On the whole the treatment is 
self-contained and is given in full detail. In general, 
the results shown are approximate and_ suflicient 
it is believed, to evaluate average densities in a 
vertical plane parallel to the cathode surface. In a 
future publication it is proposed to investigate the 
problem of cathode films with the basic assumption 
that the density gradient at the cathode surface is 
constant and uniform, rather than the density 


since the 


2. Molecular Character of Electric 
Conduction in Solutions 


Certain fundamental relations of electrochemical 
import that will be considered during the course of 
the analysis are based on the ideas of Larmor [3] in 
the matter of electrolytic conduction 

Taking an idealized situation it will be supposed 
that the electrolytic solution is a binary system and 
that ionization is complete. In each portion of the 
solution the number of the positive ions and the 
negative ions is the same and the conditions are 
steady, that is, independent of time. At a point of 
the solution let the number of ions of each kind 
per unit volume be n. Measuring positive y horizon- 
tally away from the plane cathode surface toward 
the anode, let a section normal to y be selected 
Let dN, be the number of the cations and dN, the 
number of anions crossing the section per unit area 
during time dt. The motions of the ions are due mm 





Figures in brackets indicate the literature references at the end of this paper 































































to the electrical fore e dP dy being the elec - 
potential, and in part to the forces producing 
on. If l', be the specific mobility of the cations 


the coefficient of diffusion of these ions 
dN 


or the anions, in similar manne! 


dP dn 
( nl h 


dN 
dy dy 


) dt 2 


ce there is no accumulation of electricity at a 


the relation between the electrical current 

sitv. I. and the motion of the ions will be 
se(d N dN dt, 3 

ve is the charge on each ion. Applying the 


ple of continuity of the ionic flows, 


] d d \ d \ dn 
2 dy | dt dt dt 


Ou the other hand, the current intensity remaining 
tant, from eq 6 
d?N, , d?N : 
t =), (5 
dydt dydt 


this together with eq 4 gives 


ad’ \ dn ¢ 

dydt dt’ 
d?N di = 
dy dt dt " 


Differentiating eq 1 and 2 with respect to y and 
y the relations just given, 


di d =) d nv g 
dt dy\ , dy dy 
| dn d (dP dn 
n +h y 
dt c Ty | dy “dy 
M iltiply the first of these by l 2; the second by l 
adding the resulting equations, one obtains 
dn kl h [ d n 10 
dt l li, dy’ 


Now the form of this equation is exactly the form of 
diffusion equation of an ordinary nonelectrolytic 
ition. Hence, the changes in the concentrations 
in electrolyte proceed with the diffusitivity con- 


nt 


s 


D = a) (11) 
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As this relation holds good however slight the electric 
eurrent may be, it hay be supposed that the relation 
holds in the limit when the current is nil. Therefore. 
}) is the ordinary coefficient of diffusion of the solu 
tion, and thus a relation exists connecting this coefli- 
cient with the ionic quantities, that is, with the spe- 
cific mobilities and the coefficients of diffusion of the 
two 10ns 


If the relations 
k Cl’, and f Ee 12) 


exist, Where Cis a constant independent of the con 


centration of the ions, then 
2U,l 
} ( (133) 
/ U'i+l 
Now, from the theory of electrolytic diffusion by 
Nernst [4] we have 
2U,l 
D=—. RT, 
/ / 14 


where 7’ is the absolute temperature and 2 the gas 
constant Accordingly, the constant (’ in eq 1s 
may be identified with the product RT in eq 14 


3. The Equations of Motion and of Diffusion 


Taking the case of a plane cathode surface, held 
vertically, a point at the lower edge may be chosen 
as the origin of the rectangular coordinate system 
(x,y) (see fig. 1 The positive s axis is drawn 
vertically upward and the positive y axis normal to 
the surface and into the electrolyte body Denoting 
the velocity components along the + and the y axes 
by uw and v, respectively, the equations of motion are 


ou Ou 1 Op 


Oru 
i 7 ; 
Ou OV p Ou 


y 
OV 


i] 


| and 
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1 Op 
0 4 16) 
p Cy 
where p is the pressure the kinematic viscosity 


of the and g the constant of 


gravitational acceleration 


solution, p the density 
These are the equations 
of motion used in viscous boundary laver theory 

The pressure must be expressed in terms of the 
densities Denoting the densitv of the electrolytic 
bath at points far removed from the cathode surface 
by p,, the density near the surface may be expressed 
by the relation 


p=p Ap 17 


where Ap is the deficsenev of the density in the laver 
The reduction of density is at the cathode 
surface, that is at the plaare y == Let the average 
value of Ap over the entire length of the cathode sur- 


vreatest 


face be denoted by Ap The study of the densities 
in the laver is done conveniently, as will be seen 
later, by considering the ratio 
Ap 
a . IS 
Ap 


Another quantity important for the analysis is the 


effective gray ity 


q G. 19 


Now, the second boundary laver equation, eq 16, 
states that in the layer next to the cathode surface 
the pressure along a horizontal plane remains con- 
stant. Moreover, all motions bemg nil in the main 
body of the bath solution, the pressures are hydro- 
static and hence 
pidz, 20 


P pP 


where pp is the pressure at points in the horizontal 
plane passing through the origin. Thus, 
Lop, ; 
p Ou ; p 


Remembering that Ap is a small quantity and mak- 
ing use of the relations in eq 18 and 19, we have 


| Op . 
_ / -<( (ra, 22) 
p Ou ' 
and thus eq 15 reduces to 
ou , Ou Oru . 
~ 7 Vv t bra. 23) 
Ou OV OY 


This is the equation of motion that is appropriate to 
the problem at hand. To this must be added the 
conditions of continuity, 


ou , OF 
t 0. (24) 
Or OY 
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Next will be considered the equation of diffusion 
In the presence of convective currents, the diffusion 
equation takes the form 
0 v Oo i ) 


On on on 
D{ 
Ou OY 


dr Oy 
where » is the number of particles of the diffusing 
material per unit volume. This relation reduces to 
Fick’s equation when convective currents are absent 
Although the equation is well known, it may be help 
ful to give a derivation in order to clarify the 
involved Neither the name of the 
giving the first derivation nor the manner of deriva 
tion is known to us 

Consider a parallelepiped of sides dr, dy, and 

Let the number of particles crossing the face at a per 
pendicular to the ¢ axis during time dt be 


ideas 
Investigator 


V.dydt, \ p= + nu, 
, Ou 
and the face at a dr be 
ON 
\ dudt + da dud 
/ ( V Ou ) yes 


Similarly, let the number of particles crossing th 
face at y perpendicular to the y axis be 


N2, dx dt, \ nos 
Oy 
and the face at y+-dy be 
NA dydt ( AN + 0.‘ dy )dx dt, 


so the accumulation of the particles in the parallek 
piped Is 
oN 
| )da dy dt 
Ou OY 


But the accumulation is also 


on 
; da dy dt. 
ol 


Hence, 
ON 
OY 


on ON 
of Ou 


or 


on D( Orn \ =) on 


, On Ou , O 
of or? * Oy "Or ' "5, + dy) 


OY Or OY 


The last term in the right-hand number vanishes 
because of the incompressibility of the solution 
This now proves the general equation of diffusion 
A tacit assumption of the derivation 
the volume of the diffusive 


given by eq 25 
is that the ratio of 
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ticles to the volume of the solution is a small 
tion. Stated otherwise, the result is valid for 
ill concentrations 
The diffusion equation is much simplified when 
lied to the cathode film boundary laver. First, 
concentrations do not vary appreciably with 
Second, the laver being very thin and com 
atively very long, the quantity 07n/0r? is small in 
parison with Orn OY Vor ordingl\ eq 25 become s 
. on \ On D on 
Ou Oy OY 
d since nis proportional to a, also 
Oa | Oa pre 2G) 
or Oy oy 
his is the diffusion equation appropriate to the 
oblem at hand. Equations 23, 24, and 26 govern 
he movement of liquid and the distribution of the 
nsities In the cathode film subject to the boundary 
maditions 


4. Boundary Conditions 


Ihe boundary conditions to be considered are of 

o kinds. One kind relates to the velocities and the 
ther to the densities, the latter having 
hemical implications 

The boundary conditions to be assigned to velocit, 
upponents are the most obvious ones. The solution 

the liquid of the film being viscous, uv and v vanish 

the cathode surface ; that is 


electro- 


u v i). y 0. 2é 


To state the condition for the velocities at large dis- 
tances away from the cathode surface it is necessary 
to introduce first the concept of the boundary laver 
The boundary laver thickness, 

is such that the density of the solution at this 
distance from the cathode differs only imperceptibly 
from the density of the bath. Since the movement 
of liquid in general is nil in the main portion of the 
bath, a second condition for the velocities is 


thickness (see fig 1) 


uo, Y=Ys (28) 
lt is obvious that the thickness y, of the boundary 
aver at the lower edge of the cathode surface, that ts 
1=0, vanishes. This means that the liquid moving 
ipward in the film must be replenished by a later 
and horizontal movement of the bath solution into 
the film region. Hence 
v>do, oy = (29) 
The assumption of vanishing velocities in the bath 
in the vertical direction is an idealization. Under the 
actual conditions of experimentation, arbitrary move- 
ments in the bath may be present even when forced 
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agitation is not resorted to. Since the movement of 
hquid in the film next to the cathode surface Is up- 
ward and in the film next to the anode downward, 
it is obvious that a circulation is induced in the bath 
liquid. Due to accidental variations of temperature 
at places the regularity of the circulatory motion will 
be disturbed and will give rise to weak eddving mo- 
These will be most noticeable near the upper 
surface of the bath and close to the cathode Such 
accidental cist urbances can be ipagined to have only 
minor effects on the main motion in the film and thus 
will be ignored 

The boundary conditions that are to be assigned 
to the densities are not obvious, especially the one 
relating to the value of the densities at the cathode 
surface. Regarding this, it seems that two different 
statements may be made depending on whether the 
transmitted electric current is of the limiting value 
or is inferior to this limiting value. Taking the case 
of weak electric currents, we may consider the ques 
tion of the densities at the cathode surface as follows 

Assume that the ionic system is binary and 
ionization is complete. Since at the cathode surface, 
that is at y=—0, dN,/dt—0, eq 


Llons 


» 
oO IVES 


dN / 
dt «e 


and thus eq | and 2 now become 


/ ] dP dn 
N¢é 7 dy ig dy 
and 
dP dn 
° nl dy ks dy’ 


where n is the number of ions of each kind per unit 
volume. Eliminating the potential gradient between 
these two equations and making use of eq 11, 


>> nity l (30) 
ay Nf 


This is the relation between the electric current 
intensity and the gradient of n at the cathode surface 

It is desirable that the quantity n in eq 30 be 
replaced by c, the latter being defined as the gram 
equivalent of the salt per cubie centimeter of the 
In a unit cube there are n cations and n 
anions each of ionic valence s, and the system 1s 
completely ionized. Introducing NV, the Avogadro 
number, the molar concentration is n/N and hence 
the concentration in gram equivalent of the salt per 
cubie centimeter of the solution is 


solution 


ns 
Cc ; 


and eq 30 now becomes 


pt U; ; a (31) 








Since eN is the universal constant F, with the 
numerical value F= 96450, and with the dimensions 
coulomb per gram equivalent 


de 1 / = 
a be si 


The same condition may be expressed in terms of 
densities instead of concentrations. For the major- 
ity of electrolyte solutions the variation of densities 


with the concentration is linear; that is 
de a 
mn oe) 
dp 
For copper sulphate m= 1.32 10-* gram equivalent 
per gram; and for nickel chloride m=1.71™ 10 
gram equivalent per gram. Equation 32 may now 
be written as 
dp l I 
mD ’ y 0 34) 
dy U.+U.} 


This is the relation between the intensity of the 
electric current and the gradient of the density at a 
point of the cathode surface. It is the only relation 
that may be derived on the basis of simple concepts 
of an electrochemical nature. Accordingly, we are 
not certain of the state of the densities at the cathode 
surface for the present 

Two simplifying assumptions may be made. First, 
it may be supposed that the current intensity is 
uniform at the cathode surface. The resulting 
boundary condition would be that the density gra- 
dient is uniform at the cathode surface. Second, it 
may be supposed that the density is constant at the 
cathode surface. For the present analysis we shall 
adopt the second alternative simply for the reason 
that the corresponding analysis may be completed 


with ease. Thus it will be supposed that, 


Ap= Aba, y=0, 
and 
Ap=0. y= 
In view of eq 18 
a l y=0 (35 
and 
a=0, y x (36) 


The set of equations 27, 28, 35, and 36 constitute the 
boundary conditions to be considered in solving the 
set of equations 23, 24, and 26 

Since the relation implied in eq 34 is valid for all 
conditions, it may be used to consider average effects. 
Expressing p in terms of a through the use of eq 18, 
multiplying the two sides of the resulting equation 
by dr, and integrating between z= 0 and r=/,, where 


1, is the length of immersion of the cathode surface, 


“lt (da . U, Il, 
—mDApm | iy dz U.+U, F’ 


(37) 
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where now 7 is the average intensity of the electri 
current measured in amperes per square centimeter 
We shall use this relation to compute the transport 
number of copper ions taking our data from the results 
of the Brenner experiments 


5. The Dimensional Form of the 
Basic Equations 


The solution of the basic equations subject to the 
boundary conditions indicated above will be greatly 
simplified if the equations are changed into forms in 
volving dimensionless quantities only. The clue to 
the proper course of analysis may be found in the 
method of Blasius for treating the viscous boundar \ 
laver of plates in a current of air 

The limiting values of the densities and the velox 
ities being uniform at the boundaries of the eathod 
film or of the diffusion layer, it is to be expected that 
the distribution of the velocities and of the density 
reductions in different normal sections are affine to 
each other. The affinities will be established by th 
introduction of the characteristic length 6 and the 
characteristic velocity U’, both of which are functions 
of I only 

Introducing the stream function y, such that 


oy O¥ 
u ’ ’ ’ 
OY Or 
the condition of incompressibility, eq 24, is automa 
tically satisfied. Interms of the stream function the 
equation of motion, eq 23, and the equation of diffu 
sion, eq 26, transform to 


Oy OY Ody o*y : oy Ga - 
Oy OyOr Or Oy? oY 
and 
OY Oa oy Oa D Ora 10) 
Oy Or OF OY oy" 
Introducing the new variables 
n=—y/6 1] 
and 
¥ 
> $2) 
H Us 
and remembering that J’ and 6 are functions of 2 
only, we have, first 
re) On d ld 13 
Oy OY dy 6 dyn 
and 
0 Ond ndid 44 
Or Ordn bdrdn 
and hence 
1 , , dé lH ‘ 
ows vou SS, (45 
Or dz dz dy 
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Oy dil 


[ 16 
OY dy 
O’y dl’ dH tyndi @H 17 
Ou OY da dy 0 da dy : 
oy i! aH 18 
OY o dy 
oy il 19 
oY y ad 
Substituting these « xpressions in eq 30 and 40 
dl';dH ld dtl lhl , 
ad A dy ) ' 0 da ) H dn ' a dy roa ”, 
50 
d da I d a 
/ 
| oO T ’ 
oO dy H dn f dy ' : 
y D 2) 


Now eq 50 and 51 will be freed of the differential 
rators involving s if one can write 


f { { d I 
Vl : V/ 0 J G 53 
da 6 da yo 
Vi, d i 
{'6 , 5 
i) da yo ,' 


where VW, and VJ. are numerical constants In view 
of these provisional relations, the two equations 
sumplify to 


PH l CH 1 (dil 
Yaa ty, ate, ae )+a=0, 


Il da d Gy 
M.dy dy 


In fact, the provisional equations 53 and 54, 
determine the dependence of the characteristic 
locity U and the characteristic length 6 on the 


variable z. The provisional equations 53 and 54 
imply 
dl Bias 
16 (8) (57) 
Mio Maz, 5), 7 
M,6 d (U8) # (58 
' * 2 2 
nd 


yl y aa. 





It is now obvious that l' and 6 each are single 
power terms of 4 (Assuming 


l Az’ oo 
Ba Ol 


Where A and 2 are constants involving the physical 
Substi- 
tuting these expressions in eq 57, 58, and 59, com 
paring the powers of + and the coeflicients for the 
individual equations, the comparisons vield 


characteristics of the electrolyti solution 


Vir=AL(s+r), (2) 
\V/ AB H oo 
2 | i) 4 
2 (65 
and 
| y 1G 66 
Equations 64 and 65, vield 7 2 and s=1/4 


Equation 62 now vields 
2M,—3M,=0 


For simplicity subject to the last derived equations, 
we may put MM,—1/2 and M,—1/3 
Substituting the values of the constants A and B 
that are obtained from eq 63 and 66 into eq 60 and 
61, there results 
l nd I : (67) 


and 
=(G-) =", (68) 


which are the functions showing the dependence of 
the characteristic velocity ( and the characteristic 
length 6 on the independent variable 4 

Substituting the values of MM, and VV, as found 
above into eq 55 and 56, these simplify into 


CH oy,@H o (dH 
7 dy ' 3H dn 2( dy sia ™ = 
and 
d cr ° da - 
dy? t 34 G, = sa 


which are the equation of motion and the equation 
of diffusion expressed in dimensionless forms 

In terms of the new variables, the velocity eom 
ponents uv and ¢ are, using eq 38, 45, and 46, 


dil 
ul , 7 
: dy a 


and 
d dé dil 


ls / T ’ | } 
dz VU 9H dx” dy . 











and the boundary conditions are 
- 
7-0. dif 0. n=0 
dy 
y 4 73 





The two equations 69 and 70 become the basic 
equations for the problem of temperature variation 
of a heated plate placed vertically in a cool air bath, 
when 


T—T 74 
ar ] 7 ‘ 
where 7 is the temperature near the plate, 7) the 
temperature of the plate, and 7; that of the cool air, 
and the parameter y Is replaced by 
o Vv i. (75) 


where & is the thermometric conductivity 

Goldstein [5] has discussed the Pohlhausen solution 
of the heated plate problem and the equations con- 
sidered are similar to eqs 69 and 70, when allowances 
are made for the variations in some of the coeffi- 
cients. Thus, in the two basically different phe- 
nomena of the diffusion layer near a cathode plate 
and of the conduction layer near a heated plate, 
the corresponding parameters y and ¢ play identical 
roles. The latter is well known the Prandtl 
number. The former will be called the diffusion 
parameter for the present, since it is a dimensionless 
number expressing the ratio of the kinematic vis- 
In forming 


as 


cosity to the coefficient of diffusion. 
the numerical value of y, both D and »v must be 


expressed in square centimeters per second. 


6. An Approximate Solution of the Basic 
Equations 


since an approximate solution is adequate for 
the present objectives, it is hardly necessary to 
effect a rigorous solution of the two equations, 69 
and 70. To prepare the way for the approximate 
solution, it is necessary to introduce the boundary 
layer thickness y,, which is a function of z. To 
the distances y and y, correspond the parametric 
distances » and 7,, such that 
(76) 


n=y/6 and n,=y,/6. 


We introduce the ratio 


(4 4é) 





It is obvious that the quantity @ will serve as the in- 
dependent variable in the place of n Indeed. the 
basic equations 69 and 70 transform to 


CH dil 





Pil 
~ + An Dy -n o@ } 7 R) 
dé Hl de ~7 ( de ) n (} is 
and 
d a P da 
aT TO” H da 0, ret) 
with the boundary conditions 
+ 
H1—0., A—O, 
dH () i ) 
dé 
dil > S( 
0. “” ] 
dé 
a ] “7 () 
a i), 0 3 
4 


in place of the set shown by eq 73 

It is obvious, intuitively, that the space rate of 
change of the velocity components, Ou/Oy and Ov dy 
each vanish at large distances from the cathode sur 
face. These conditions are fulfilled likewise quit« 
approximately at distances equaling or exceeding y 
This observation suggests the additional boundary 
conditions 


PH a o=1 
dé ‘ ’ 
S| 
da 0 a ] 
dé ’ 
the first of these resulting from eq 48 All the 


boundary conditions of eqs 80 and 81 are satisfied 
if it is assumed that 


T 
sin a. Pe 


Tl Ale cos 7ré COs 2x0. 


The above analysis involves the quantities A 
and »,, which at present are not known. For the 
determination of their values the differential equa 
tions 78 and 79 are at our disposal. Multiplying 


and 


S83) 


| these equations by d@, integrating between the 
| limits @=0 and @=1, and making use of the boun 


| dary conditions, the result is 


. {* (di 
}- Sm, ( a 


\ 


) d0-+n3 (s4 


(CH 
| —1( a9 


*] 
adé=0, 


J0 








yin 
the 


e of 
OY 
sul 
puile 


r 


mm Y 


lary 


s 


the 
she d 


3 A 
the 
qua 
ying 
the 
oun 








7 


he 


ver problem based on the assumption that 


<> 12. ee 


dé 


resorting to the rules of integration by parts 
oducing the values of @ and // from eq 82 and 


id effecting the integrations. we find 
2°) r—2 
2 ay | ” ( —) 0 
‘ T 
T 14 1 
+ - n ‘) 
2 ) Tv 
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Y *y wT Sh 
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expressions for and ci) as given by eq SO 
S7 respectively the eXpressions 
d // as given by eq 8? and 83 respectively, con 

ite the first approximate solution of the diffusion 
the 


and for a 


iction of densities or concentrations, at 


hode surface is uniform 


Thickness of the Cathode Diffusion Layer 


Inasmuch as the diffusion parameter y= v//) is o1 
narily a very large quantity, it ts sullicient to write 


thre plac c of eq St 


troducing the value of y, from eq 76 and the value 


find first 
V7 y , 1; ; 
| 14(r—2 Gy} * 


irom eq OS. we 


next introducing the value of G from eq 10 
‘ pit 
y 3.00 I Sv 
(2) 
This is the expression that vlVves the thickness 
the diffusion laver along the cathode surface 
e hig. 2 The thickness is not uniform, but 
reases from the zero value at the lower edge of 


he cathode surface to a maximum at the upper end 


the cathode surface still in contact with the 
ctrolvte The rate of increase of the thickness, 
Wever, Is verv gradual due to the fact that the 


' 


| 
I 


«wkness if proportional to the one fourth power 
the distance from the lower end of the cathode 
rain at a chosen point s the thickness of the laver 
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the maximum reduction of the 
smce 


deere uses W hen Ap 


density at the cathode surface, is increased 


the decrease of density at the cathode surface is 
proportional to the electrical current density, when 
the latter is inferior to its limiting value, it is to 


be expected that a decrease in the intensity of 
current will cause the thickness of the diffusion laver 
to be augmented 

Having established the 
laver thickness, it interest 
the normal density 
The 


understood if it is expressed in terms of the average 


magnitude of the diffusion 
to riquire next 
the cathode surface 
best 


is of into 
gradient at 
density will be 


variation of the evradient 


density gradient, the averaging being made over the 
entire cathode surface In view of eq 1S the quan 


tity 


Oa Oa > 
(3°) (32) -2 (90 


OY 


where the numerator on the left-hand side refers to 
local values at the point 2 and the denominator to 
the 


since 


average values, is the ratio to be considered 


Va Oa OF 
OV Ov OV 


we have from eq 77 and S82 


Oa us 
( Sy ) » 


~y 
or 
( Oa ) Cy 9] 
OY J 
where 
us py, 
4)” 
6 66s YAP m4 ) . 
Accordingly the average value over the length / 
of the cathode surface ts 
Oa ' 
(5°) : 0 
OV 31 
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and the desired ratio is 


R=". ( . ) > 04 

H\7 
which gives the proportionate variation of the normal 
density gradient over the cathode surface The 
variation is shown graphically in figure 3 In the 
section of the cathode from 2«//,—0.2 to r/,;=1.0 the 


ratio is nearly constant This circumstance assures 
the utility of the present solutions even for the case 
in Which the normal concentration gradient is strictly 
constant over the entire surface of the cathode, pro- 
vided that one deals with the average values of the 
densities in the planes close and parallel to the 
cathode surface 

If the diffusion layer thickness at r=/,, that is at 
the upper end of the cathode surface just in contact 
with the bath, be denoted by y then the average 
value of the normal density gradient, expressed in 
terms of the maximum density reduction, ts 





Oa 2.1 . 
(3, 7 


— 
t 


v] 


Use will be made, subsequently, of this relation 


8. Velocitus in the Diffusion Layer 
Consider first the velocity component u, that is 


the component parallel to the surface of the cathode 


Now, 


dil “ 
u=Tl de’ il 
where 
l 2G"? (67 
Since 
dil dH dé 
dy de dn’ 
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Fiaure 3. Variation of the density gradient along the cathode 


suriace. 


164 


eq 71 may serve, in view of eq 77, 83, and 87, to 
give the distribution of the velocities in a norma! 
section of the diffusion laver. There is, however 
another method for representing the velocity dis- 
tribution Moving away from the cathode surface 
in a normal section, uv is zero at the cathode surface. 
increases to a Maximum value uw, and then de- 
creases to zero at the limit of the boundary layer 
that is at y=y Thus, the distribution of the 
velocity may be obtained also by considering the 


ratio uwiu Now 
dil 
u ( - ) O6 
In terms of the variable @ 
dH 1 (dH : 
(i, ) n ep ; oe 


the value of @ that makes uv a maximum, w 


l 
A sin rO,+ , sin 270 


_ 


If 0 be 


have from eq 83, 


” 


/ 


dil 
(u,) 


The value of @ satifies the relation 


cos rh, cos 2ré 0) 
This makes 
sin 76,+.. sin 278 1.299, 
) 
and therefore 
dil 
( ) 1.299 4 98 
dé 


Substituting into eq 96, from eq 98, 97, 88, and 67 


u,,=0 s2(“ ) Pu 


? 


gives 


or using eq 19, 


Uy 


Ap, . , 
v7 : ) : 90 


0 82( 

Pry 
Accordingly, the maximum value of uw in a normal 
section changes from one section to another and 
increases with the one-half power of the distance 
from the lower end of the cathode. 

Since, 

u/U», = (dH /dn)/(dH/dy)m. 

and introducing the relations in the equations 97 and 
98, the distribution of the velocity component u in a 
normal section in the diffusion layer becomes 


’ 0 Bi. 1 
U/Um =0.77 [sin T u +; sin 29 Ai (100) 
Ys « Ys 
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The graphical representation of the distribution ts 
ven in figure 4 


Next consider v, the velocity component normal to 


he cathode surface. The velocity can be written 


dé ; dil 
l 3H —y1 
da | i 7 | 101 


rhe magnitude of vis largest at the liquid boundary 
of the diffusion layer Let be this largest value 
This occurs at 7 n Let 77, be the value of // at 
’ or at @ l At n Ns, dil dy 0, and thus 


{ =) 
rom eq ie. US 


r 33 HU, (102) 
where 
H,=2A;/z. (103) 


In view of the values of ( and 6 from eq 67 and 68, 


dé 1 
l (4°0)'*r-'" 
dr 2 ‘ 
Substituting the above result in eq 102 and also 
introducing the value of A, from eq 87, 


iP l 013(°4@ ty Sy : (104) 


Ler ordingly, the liquid of the bath moves hori 
zontally into the cathode diffusion laver The veloe- 
ity with which the liquid approaches the cathode 
surface is variable along the length of the cathode 
surface, the maximum value occurring at the lower 
end of the cathode surface After a certain elevation 
is reached the variation in the value of v7, ceases 
to be pronounced 
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Y/Ys 
Fiaure 5 Dis / m of the ho onta , omponent 
The representation of the distribution of in a 


normal section of the laver is most conveniently 
done by considering the ratio Now 


| n dil 
H G 3 1b 
l i dil 
HA #-5 a 


This gives, in view of eq 83 and 103, 


or mm terms ot @ 


r 5 | l ri) ri 
cos wi cos 2ré , sin wl sin 27r@ 
S r 1? 


1O5) 


The graphical representation of the distribution ts 
shown in figure 5 


9. The Brenner Freezing Method and the 
Laminar Mean Densities 


In the Brenner method [6] the electrolytic solution 
hear a evlindrical cathode Is suddenly frozen, Succes 
sive lavers of the solid of 0.003- to 0.004 inch thick- 
ness are carved off on a lathe, and the compositions 
of the layers are determined. Thus, in this method 
the average concentration or the average density of 
a lamina of thickness dy,, of distance y, from the 
cathode and. of length /, are determined Accord 
ingly, if theoretical results are to be compared with 
the observed results of the Brenner experiments, tt 
becomes necessary to derive a formula giving the 
average value of densities of a lamina of distance 
Y from the cathode surface, 
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Let the average density reduction in the lamina of 
length J, and of distance y, from the cathode surface 
be denoted by Ap’ Mathematically 


2. 
Ap j Apda (106 
Let /, be the value of x for which the boundary laver 
thickness equals y, (see fig. 6 Obviously, for values 
of z l,. Ap 0, and thus, 


l 


Ap’ Apda 107 


Since, according to the present solutions, the density 


reduction at y= y, is given by the relation 


Sp =Apm( 1—sin > 108 


(from eq 82), then 


‘ wy a xy 
Ap’ = Ap, ; | i da j |, sin 9 - 


dz | 109 


It must be remembered in effecting the integration 
of the second integral of the right-hand member, that 
y is held constant Putting 


and 
(110) 


It is seen from the law of the laver thickness, eq SY, 
that 


l,=(,/ 6%, 111) 
/ l,/p" (112) 

and thus, 
di 1/,d8/8 (113) 


Accordingly, in terms of the new variable 8. 


, 4 ‘ *B T dg 
Ap Ap l re Te) sin; re] 3 114 
For convenience, we may write 
“6 rT ig 
F(8,)= Bt—4, sin ” B- 115 
J 2" 8B 
and hence, 
Ap’/Ap | kip 116 
Since B v, Yo, We have also 
Ap’/Ap,»=1—Fiy,/y 117 
This is the desired formula that will be used in 


analyzing the experimental data of Brennet It 
involves the average reduction of density over a 
lamina, the average maximum reduction of density 
at the cathode surface, the distance of the lamina 
from the cathode, and the maximum thickness of 
the diffusion layer The determination of the 
function F(8,) was carried out by replacing the term 


sin (#/2)8 appearing in the integrand by a power 
series in 8 and then effecting the integration. The 
results of this determination are given in table | 


10. An Analysis of the Results Obtained by 
Brenner with the Freezing Method 


The comparison between the results of the present 
analytical theory of the cathode diffusion layer and 
the results of the Brenner experiments on the 
reduction of the densities in the diffusion layer will 
be carried out along the following lines. First, from 
the observed values of the reduction of densit Vv al 
the cathode surface, Ap,,, the maximum thickness 
of the diffusion layer, y,.,, will be computed and 
compared with the observations. Second, using the 
same observed values of Ap,,, the ration Ap’/Ap 
will be computed for the distance y,, and these will be 
compared with the observations. Third, Ap,, will be 
computed from a consideration. of the electric 
current intensities used in the expeiments and will 
be compared with the observations 


Pane | Numerical values of the function f ¥i/¥ 
Fly I 
rT) 0. 0000 oft 0. 8026 
ol w70 7 JSS 
wre4 s O40 
us ' 17" 
‘ 7M Lo 1 ooo 
RISO 


The Brenner [6] results relate to electrolytic baths 
of copper sulfate and nickel sulfate. The baths 
actually used were not pure binary systems, but 
this fact will be ignored. The complete data from 


the Brenner tests are collected in table 2. The 


166 











the theoretical formula in eq S89 





ABLE 2 Diffusion layer densities of the Brenner test 
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entrations in the original data are in terms of 
normalities The densities as derived from the 
rmality values are added. The average densities 
taken from the table are plotted against y, in 
ire 7, and smooth curves are drawn through the 


ts of observation The extreme values p and p 


the densities are read from the curves The 
rences (p po) give Ap,,, the average maximum 
iction of densities at the cathode surfaces The 
Terence p p’ where p is the observed value 


the densities corresponding to the distance y 
es Ap’, the average reduction of the densities 


responding to the distance y All these quan 
es, together with the ratio Ap’/Ap are entered 
table ) 


From the observed values of Ap and p the maxi 
im laver thickness, y» mav be computed from 
Hlowever, the use 


" 02 
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of eq 89 for this purpose requires that the parametric 
diffusion number, y, be known. This number in- 
volves the ratio of the kinematic viscosity nu, to the 
diffusion constant J), for each of the two electrolytic 
solutions. The viscosities of the electrolytes may 
be determined from the normalities using the 
Arrhenius [7] formula 


where M, Is the specific Viscosity of a normal solution 
referred to the solvent at the same temperature 
uw is the specific viscosity of the solution having the 
normality N Wagner |S] gives u 1.358 for copper 
sulfate and yg 1.361 for nickel sulfate at the tem 
perature 25° C. The normalities of the bath solu 
tions in the Brenner tests being 2 and the densities 
about p=—1.151, the kinematic viscosities of the two 
solutions are practically equal and so gu 0.0144 
cm sec Regarding the diffusion coefficient J), the 
Thovert [9] values for copper sulfate may be sup 
posed to apply as well to nickel sulfate Kon copper 
sulfate at a temperature of 17° C, N= 1.95, the dif 
fusion coefficient is J) —2.66™ 10 cm sec Henee 
the value 


0.572“ 10' 


As regards the 
cathode immersion length /,, we take it to be 10 em 
this being the value communicated to the author 
orally by Brennet 

l sing the above-mentioned values of the pertinent 


mav be used for the two solutions 


quantities the magnitudes of the maximum lave 
thickness were computed from eq 85, and these are 
shown in table 4 Also there are given in the table 
the observed values of the same quantity as read 
from the curves in figure 7 At instances equalling 
the maximum laver thicknesses the nverave cdlensities 


equal the bath densities There is some difficulty 
on layer der oe af the I 
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Ficure 7 (ra phical representation of the data from Brenner's 
teats 


in identifying the exact magnitude of the maximum 
laver thickness from the curves of figure 7. If this 
element of uncertainty be ignored, the agreement 
between the theoretical and observed values may be 
judged to be fairly satisfactory 

The second question to be considered is the dis- 
tribution of the mean densities in the layer itself, 
that is the variation of Ap’/Ap, with y;/y.. The 
observed proportional density reductions for a given 
relative distance are given in table 3. The method 
of determinations is obvious from the contents of 
the various columns. In forming the ratio y,/y., 
the theoretical values of y,, are used. Values from 
the table are plotted in figure 8. The curve drawn 
is the theoretical distribution and is the plot of eq 
117. Again the agreement between observation and 
theory may be judged to be fairly satisfactory. 

The third question to be examined is the depend- 


Tanie 4 
Bath Cathode Reduction of Cathode 
density density density at immersion 
cathode length, 

p o de h 
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Ficure 8 Distribution of densities of the diffusion 
according to theory and observation 


ence of the average density gradient at the cathod 
surface upon the intensity of the electrical current 
This question deals with the electrochemical bound 
Since the 


ary condition represented by eq 37. 
average gradient according to eq 95 ts 
(=) 2 l 
Oy /o Ys 
eq 37 may be changed to 
I 
2.1m DAp, -—— , (118) 
: U,+U, F 
Ya 
This relation may be examined in two ways. If the 


transport number of the cations (nickel and copper) 
are known, the relation may be used in computing 
Apm, the average density reduction at the cathode 
Or, using the observed value Ap,, the transport 
number of the cations may be determined. The 
second method will be used here. 

From the data in table 3 the ratios Ap,,/y., are com- 
puted and these are plotted against current density) 
in figure 9. The distribution of the points is linear, 
and the equation of the line drawn 1s 
(119) 


Ap» (Yt) 0.5 <1. 


Theoretical and observed maximum diffusion layer thickness in the Brenner experiments 


Observed Theoretical 


Current Metallic 
maximum layer maximum layer 
thickness, thickness, entonatty —_ 
v.. + I 
em em amp/icm 
0.032 0. 302 o.04 Copper 
as2 O3&2 ol Nickel 
037 0347 02 Do 
a4 OS15 4 Do 
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tuting these values in eq 118, it is seen that 
O.356 120 


ording|y the transport of the cations, nickel, is 
0.644, a value verv close to the generally ac- 
ted value determined by other methods Thus 
his respect also there is agreement between theory 
observation 
Another element of the theory that could be ex- 
ined with considerable signficance is the question 
he vertical velocities. Unfortunately in the ex 
porupents of Brenner this matter was not touched 
lo give an idea of the order of magnitude of the 
velocities involved, We may compute the maximum 
ward velocity In a normal section of the diffusion 


er using eq QQ Consider the case of copper sul- 


fate. Corresponding to the current intensity of 
1-410? amp/cem’, the density reduction at the 
cathode surface is Ap,,/p;=— 5.31107 since 1/y 
1.74107, the formula gives for maximum velocities 
Un 
» 50X10 
vg2 


Thus at a point r=10 em, u, equals 2.48 mm/sec 

The discussion presented above must lead one to 
the following conclusion. The approximate analyti- 
cal theory of the diffusion layer based on the assump- 
tion of a uniform density reduction along the entire 
cathode surface gives fairly satisfactory results 
regarding the maximum width of the laver, the aver- 
age densities in planes parallel to the cathode surface, 
and the transport number of ions. Brenner is now 
engaged in further work on the variation of densi- 
ties In planes parallel to the cathode surface. In- 
formation of this kind may possibly point to defects 
in the present theory. In that event it might be 
hecessary to base the theoretical investigation of 
the problem on the assumption that the gradient 
of the density and not the density is uniform along 
the cathode surface 


lhe author expresses his sincere appreciation to 
W. Blum for the guidance shown in the beginning of 
the investigation, and to A. Brenner for advice given 
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Absorption Spectra of Thioindigo Dyes in 
Benzene and Chloroform’ 





Wallace R. Brode and George M. Wyman 


Phe spectral absorption curves of 10 purified thioindigo dves, containing mostly methyl 
and ethoxy! groups or halogens as substituents, in benzene and chloroform solution were 
determined over the raviolet and visible parts of the spectrum The existence of an 
equilibrium bet weet is and trans isomers was indicated by the study of the effect of light 
on the absorption spectra. Some of the results of this investigation have been published 
rhe present paper supplements the earlier paper and contains the spectral absorption curves 
that have not previously been published, along with the data on the purification of the three 
dves that had not been reported there In general, the first (long-wave) absorption band of 
the cis-form of each dye occurs at ashorter (by approximately 60 mu) wavelength than that of 
the corresponding isomer rhe spectra are displaced slightly toward longer wave 
lengths in chloroforn i than in benzene Substituents in the five and seven positions 
cause a bathochromic shift in the spectra of both isomers of each dve in each solvent without 
affecting the shape of the spectral absorption curves The introduction of substituents in 
the four or six positions results in a change in the shape of the absorption curves and a shift 
of the absorption to shorter wavelengt! 


l. Introduction 


The absorption spectra of the 10 purified thioindigo 
dyes listed in table 3 were determined as a part of 
the National Bureau of Standards research project 
on the spectra of pure dyestuffs. It was observed 
that the absorption spectra of these dves in organic 
solvents were a function of the wavelength range 
of the illumination to which the solutions had been 
exposed prior to the measurement of the spectrum 
This phototropic behavior was attributed to the 
existence of a dynamic equilibrium between the 
as shown here for thioindigo 


cis and trans forms 


} TRANS 


It was also possible to separate the two isomers in 
Some of the results 


the case of two of these dves 


of this investigation on seven of these dves have | 


been published (see footnote 2 


The object of this paper is to reproduce the com- 
plete spectral absorption curves of the dyes for which 


only the position of the absorption bands was 
reported in the earlier paper (cf. tables 2 and 3 of 
the reference cited in footnote 2); to describe the 
purification and the absorption spectra of three 


additional dyestuffs of this type; and to attempt to 
establish further correlations between the spectral 


absorption and the structure of these molecules. 
Paper X XIII in the serieson “The relation between the absorption spectra 

and the chemical constitution of dye Paper XXII referred to in footnote 2 
Gg. M. Wyman and W. R. Brode. J. Am. Chem. Soc. 73, 1487 (19051 
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2. Experimental Methods 


2.1. Purification of Dyes 


commercial 


work were 


The dves used in this 
products of known structure 
by adsorption on silica gel followed by 
elution. The solutions were then 
dryness and the residue reervstallized from benzene 
or chloroform. Table 1 lists the conditions of the 
purification of the three dyes that had not been 
reported previously 


was purified 
selective 


evaporated to 


Taste 1 Chromatographic purification of thioind 
i , 
) ae ! I iw ; 
- . 
retramethylthio- 3 CCl | 
’ J ( ( ( 
tig ve ' {Col jOHu 
1.4 7 Hexachloro {87.5 of CCk } 
cc) CHe 
t nd tl2 CH j i 
“yr rece i 
Vat Sea G r CH pO CHe 


2.2 Preparation of Solutions 


Approximately 0.01-g samples of each dye were 
weighed and dissolved in 600 to 700 ml of benzene 
or chloroform by heating at the boiling point of the 
solvent for several hours. The solutions were then 
cooled to room temperature and diluted to 1 liter 


2.3. Irradiation of Solutions 


The dye solution contained in a fused quartz ab 
sorption cell was exposed for 5 to 10 min to colored 
light, using a 100-w Spencer microfilm projector with 
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Ing vlass filters These filters are listed in 
2, along with the manufacturer's numbers and 
observed transmission cut-off (wavelength at 
h the transmission of light is less than 1° 


2.4. Measurement of Absorption Spectra 


he absorption spectra were measured at room 
perature by means of a Cary Recording Quartz 
trophotometer Model 12), using matched 
orption cells equipped with fused quartz windows 


3. Results 


Figures 1 to 16 show the experimentally deter- 
ed absorption spectra of eight of these thioindigo 
s in benzene and in chloroform solutions under 
stated conditions of irradiation. Figures 17 to 
> present the calculated probable absorption curves 
for each pair of pure cis and trans isomers of each 
which had not previously been reported, in the 
two solvents (concentration: 0.0100 g/liter, cell 
ength: 2.00 em Table 2 lists the approximate 
amounts of each isomer present at equilibrium, under 
various conditions of illumination, for the three dyes 
that had not previously been reported 
The spectral absorption curves in figures 1 to 16 
demonstrate graphically the phototropic effect dis- 
cussed in the previous paper (see footnote 2). This 
effect, attributed to the existence of a dynamic 
equilibrium between the cis and trans forms of these 


dyes, is apparently common to all thioindigo dyes 
BLE 2 ipprorimate amounts of the two isomers of each dye 
Wa In benzene I 
filte i eT tl 
. . nes . 
tra nad i u)® 
M4 it 2 
Ww xu i! ~ » 
‘ al ‘ 
‘ ~ 1 “4 ‘ 
4 Hex ! n tN TDs 
s 22 sl 1v 
x 44. a2 Is 4 6 
a 2 as ‘ aan 
J 4° Ly 45 
Vat Sea i Ths 
a” 2 ‘ 
n i 6s {2 ‘ 
y 4 “) “) 1 wv 
aA » 72 an 1 
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in benzene and in chloroform solutions, and it was 
observed that even the hemi-thioindigoid dye, Vat 
Searlet G (Colour Index No. 1228), exhibited this 
behavior. It is believed that interconversion of 
cis-trans isomers in general can be effected by selec 
tive irradiation, provided that the first: absorption 
bands (at the long-wave end of the spectrum) of the 
two isomers occur at appreciably different wave 
lengths. It can be readily seen from figures 17 to 32 
that the cis isomers of these dyes have their first 
absorption bands at a considerably higher frequency 


than do the corresponding trans forms. This 
frequency difference is: remarkably constant and 


ranges from 1,800 to 2,250 em~' in chloroform: and 
from 2.000 to 2.400 em in benzene solutions for the 
10 dyes reported here, as shown in table 3.4 

The absorption spectra of each dye in the two 
solvents used appear strikingly similar. There is a 
slight displacement of the first absorption bands 
toward lower wave numbers in chloroform. This 
displacement amounts to less than 200 em~! for the 
trans isomers and 200 to 400 em~'! for the eis forms 
It seems probable that the more polar solvent (chloro 
form) lowers the energy of the first excited states 
(which are also polar, cf. the reference cited in foot 
note 2), and this effect is more pronounced for the 
relatively more polar cis form than for the correspond- 
ing frans isomer. The ratios of the absorption bands 
are constant for each isomer and, for each dye, inde- 


pendent of the solvent (ef. table 3) 

It is that, although in’ the 
trans isomers, the most intense ultraviolet absorption 
peak is always at slightly less than double the fre- 
quency of the first absorption band (corresponding 
to the second harmonic of that fundamental vibra- 
tion), this ratio is diminished considerably (to ap- 
proximately 1.6) in the cis isomers (ef. table 3) 
Substitution in the 6,6’-positions reduces these ratios 
considerably and, in addition, 
absorption bands in the near ultraviolet 


interesting to note 


rives rise to new 

It is possible to evaluate the effect of substituents 
by comparing the calculated absorption curves for 
the pure cis and trans isomers (figs. 17 to 32). It is 
immediately apparent from such a comparison that 
that the introduction of substituents in thioindigo in 
the 5 and 7 positions causes the usual shift toward 
longer wavelengths in both isomers The effect of a 
methyl group is approximately equal to that of a 
chlorine atom (ef. figs. 17 to 23). This bathochromic 
shift is also observed when two bromine atoms are 
substituted in the 5,5’-positions in 6,6'-diethoxy- 
thioindigo (ef. figs. 27 to 29) 1,4’, 5.5’, 7,7’- 
Hexachlorothioindigo also shows this expected batho 
chromic shift in the ultraviolet region; the first 
absorption band in the visible region, however, is at 
about the same wavelength as that for tetramethyl 
thioindigo. In its ultraviolet spectrum it shows the 


absence of the second, auxiliary absorption band in 
benzene solution only, thus contributing the sole 
‘ For an interpretation of tl fre juency difference cf the references n footnote 
An explanation wa uggested for these bathochromi hift n the earlier 
paper (footnote 2 
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example in this series for a compound showing a 
significant difference in absorption spectrum in the 
two solvents used (ef. figs. 7 and 8, also 24 and 25) 

It had been pointed out by Formanek ° that the 
introduction of substituents in the 6,6’-positions in 
thioindigo causes a shift of the first absorption band 
toward higher frequencies. Although the existence 
of cis and trans isomers of thioindigo dyes was not 
known at that time, this statement has been shown 
to be applicable to the cis compounds also and an 
explanation for this, apparently anomalous, hypso- 
chromic shift was suggested in the previous paper 
(see footnote 2). In general, the introduction of 
substituents in the 6,6’-position in the thioindigo 
nucleus gives rise to the appearance of new absorp- 
tion peaks in the near ultraviolet region (ef. figs. 
26 to 29). In addition, the most intense ultraviolet 
absorption band undergoes a slight shift toward 
lower frequencies in both isomers of each of these 
dyes (except frans 6,6’-diethoxythioindigo). Thus, 
it appears likely that in the second excited state, 
the 6,6’-positions are more strongly conjugated with 
the negatively polarizable oxygen atoms, than they 
are in the first excited state. A comparison of the 
spectra of 5,5’-dichloro-7,7’-dimethylthioindigo (figs. 
17 and 18) and of 4,4’-dimethyl-6,6’-dichlorothioin- 
digo (figure 26 and figure 7 in the reference cited in 
footnote 2 shows, however, that the bathochromic 
shift of the main ultraviolet absorption band in the 
case of the latter is still far less than in the case of 


137 (1928). 


*J. Formanek, 7. Angew, Chem, 41, 
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its more “normal” isomer. It is interesting to note 
that the visible spectrum of bis-4,5-benzothioindigo 
(fig. 30 and fig. 9 in the reference cited in footnote 2 
resembles the spectra of 6,6’-diethoxythioindigo and 
of 5,5’-dibromo-6,6’-diethoxythioindigo more closely 
than that of any of the other thioindigo dyes. This 
suggests that the two phenyl groups act as electron 
donors similar to (though not as strong as) the ethoxy 
groups; contributing structures of this type to the 
resonance hybrid of the ground state 





® 
08 0° 
C — X C, . 
®@ \ | a . bf sal \ 
c=¢ c==Cc 
* a? y Y ee 
oe bo VU 
& 
~~ 
’ ( 7 
aN / P rf por / 
C==rC C=C 
af \cAy y o ve - 
oe : } 





The strong ultraviolet absorption of bis-4,5-benzo- 


thioindigo (fig. 13) is probably due to the strong 
absorption normally associated with phenyl groups 
in this region of the spectrum. 
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A Sonic-Flow Pyrometer for Measuring 
Gas Temperatures 


George T. Lalos 


Evaluation of the performance of gas turbines and jet engines requires knowledge of 


the temperature of the working medium Although conventional temperature-measuring 
instruments are adequate for some purposes, they possess a number of disadvantages that 
become more serious with increasing gas velocity and temperature The most troublesome 
sources of error are the directed motion of the Las and the loss of he at to the surroundings 
by radiation and conduction An instrument that is relatively free of these errors is described 
and evaluated It has been called a sonic pyvrometer because the measuring element is 
located in the throat of a nozzle through which the gas, whose temperature is to be measured 
is caused to flow at sonic velocity In accordance with theory and experiment, the total 
temperature of the gas is obtained by multiplying the temperature indicated by the measuring 
element by an essentially constant correction factor Tests on experimental units indicate 


that these have many advantages over conventional instruments when used in the evaluation 


and control of gas turbines and jet engines 


l. Introduction The relation between these two temperatures of a 
flowing gas can be expressed by the well-known 
While conventional instruments for measuring | ¢qUations 


the temperature of the working medium of gas 


T,.—T,=V?/2gJC,, 


irbines and jet engines are adequate for some and 

purposes, they have certain inherent disadvantages VT. 1+ (y—1) M2 - 
that become more serious with increasing gas temper- 

iture and velocity. Particularly troublesome are | in which 7 is absolute temperature (°R), V is velocity 
the effects on temperature-sensing elements of | (ft/sec), g is the acceleration of gravity (ft/sec*), J 1s 


the mechanical equivalent of heat (ft Ib/Btu), (, 1s 
heat capacity at constant pressure (Btu/lb deg R), 4 
is the specific heat ratio, Wis Mach number, and the 
subscripts s and ¢ signify static and total 


directed motion of the gas and of heat loss to the 
irroundimgs by radiation and conduction 
As will be shown, these effects become negligibly 
small in an instrument of the type described In 
such an instrument the measuring junction of a | respectively 
thermocouple is located within the throat of a small When an error-free, conventional temperature- 


nozzle, through which the gas whose temperature is | Sensing Instrument ts immersed in a flowing gas, it 


to be measured is caused to flow at sonic velocity reaches and indicates a temperature (7',) inter- 
For the latter reason. such a modification of the | mediate between 7, and 7 Neither 7, nor 7, can 
familiar suction pyvrometer 1s here called a sonic be deduced from an observed value of 7, without 


pyrometer further knowledge of the characteristics of the 


When the velocity of the gas in the free-stream | instrument, 
Experiments have shown that each actual instru- 


ment, in the absence of other errors, attains a tem 
perature that exceeds 7, by a nearly constant fraction 
of the interval 7,—7,. This fraction, known as the 
recovery factor (r), is a function of the configuration 
of the instrument and changes little with the com 
position, temperature, pressure, or velocity of the 
flowing gas. The above definition of recovery factor 


exceeds a few hundred feet per second, the kinetic 
energy of directed motion becomes an appreciable 
fraction of that due solely to random molecular 
motion For a complete energy accounting under 
these conditions it is therefore necessary to recognize 
two temperatures, pressures, and enthalpies in 
specifying the energy state of the moving gas. These 
are frequently identified by the adjectives static and 
total. As an example, the static temperature would 
be indicated by an error-free thermocouple moving 
in the same direction and at the same velocity as the 
gas, that is, by one that is static with respect to the 
gas; the total temperature would be indicated by an 
error-free instrument immersed in the gas ~fter the 
latter is brought to rest adiabatically. 


is expressed by the equation 
r=(T,—T,)(T,—T 5 
Combining eq 1 and 3 gives 


T,=T,4+(1—r) V?2/2gJC,, 





Phe material contained in this report has been accepted by the Catholic | from which it is seen that the actual temperature of 
t t f America in partial fulfillment of the requirements for the degree of . 
; om a high-velocity gas stream can be calculated from the 
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temperature indicated by a conventional instrument 
only when the velocity and heat capacity of the gas 
and the recovery factor of the instrument are known 
This requirement obviously 
limitation to the application of conventional sensing 
elements at high gas velocities 

It is to be noted that the higher the recovery factor 
of the instrument, the smaller is the difference 
between 7 and 7, and the less is the need for 
accuracy in the values of V and (,. Total-tempera- 
ture probes, which indicate 7, directly, are in com- 
mon use in the low-temperature range llowever 
they are not practical at high temperatures and under 
the extreme mechanical stresses existing in gas tur 
bines Nevertheless, an instrument designed for 
engine applications should have the highest prac- 
ticable recovery factor 

Conventional instruments have a second limitation 
that can be understood by considering the heat 
transfer to and from the surroundings. In most 
applications, the gas turbine being a common one, 
the retaining walls are colder than the gas. A measur- 
ng instrument immersed in the stream thus attains 
a steady temperature below that of the gas, because 
energy 1s lost by the sensing element to the walls by 
radiation and conduction. The temperature indicated 
by the element is determined by the relative rates of 
heat transfer from the gas by foreed convection and 
to the surroundings by radiation and conduction 
The heat exchange with the surroundings is a func- 
tion of the temperatures of the element and the 
walls, and that part that is due to radiation increases 
approximately as the fourth power of the gas tem- 
perature From the standpoint of thermometric 
error, it is obviously desirable that the rate of 
transfer to the element be high as compared with the 
rate of loss 


constitutes a serous 


For a conventional element installed in a gas 
turbine, no independent control can be exercised over 
the mass rate at which gas flows over the instrument, 
which in turn determines the rate of heat transfer to 
the unit by forced convection. Actually this rate 
is a function of the operating conditions alone, and 
these, for aircraft gas turbines, involve flow rates 
that decrease with altitude. Thus, although the 
error of a conventional instrument due to radiation 
and conduction could be reduced by increasing the 
mass flow past the sensing element, not only is it 
impossible to exercise such control in an engine, but 
also the accuracy of the measured temperatures de- 
creases with decreasing ambient pressure (i. e. with 
increasing altitude) 

From the above facts, it is clear that improved 
performance and simplified over-all instrumentation 
may be possible by modifying conventional sensing 
elements so that (a) the recovery factor is increased ; 
(b) the rate of heat loss from the sensing element by 
radiation and conduction is decreased; (c) the mass 
flow across the element is as high as practicable; 
and (d) the temperature indicated is independent 
of the free-stream velocity. 

It will be shown that a sonic pyrometer can be 
designed to have most of these advantages, and that 


it therefore deserves careful consideration in the 
instrumentation of gas turbines. Previous reports 
on the sonic py rometer appeared tts Classified docu- 
ments, and more recently a paper published by Allen 
and Hamm [1] described one instrument of this type 
The present detailed study was undertaken because 
none of the previous reports contain a careful evalua- 
tion of the several important design variables 

The present results, obtained with three experi- 
mental units, are in excellent agreement with theory, 
and indicate that a single calibration in air at a rela- 
tively low temperature, together with a knowledge 
of the y of the gas in which the instrument is applied, 
suffices for calculating reliable values of 7, from 
observed values of 7, over the range experienced in 
gas turbines 


2. Theory of the Sonic Pyrometer 


Sonic gas velocity (V.) is always maintained over 
the measuring junction of the sonic pyrometer 
This velocity may be expressed as follows in terms of 
the properties and state of the gas 


V2. = yghT,, (5) 


in which # is the gas constant in ft lb/lb deg Rankine, 
and the other quantities have been defined previ- 
ously 

Since C,—C,= R/J, eq 1 and 


of sonic flow may be combined to give 


5 under conditions 


T,=2T,/(y¥+1) (6) 


Combining eq 4, 5, and 6 gives 
T = (¥+1)7,/[2+r(y—1)]=aT,, (7) 


where a= (y+ 1)/|24+r(y—1)] 

Calibration and use of the sonic pyrometer would 
be extremely simple if the @ of a given instrument 
were completely independent of operating conditions. 
This cannot be so, because y decreases somewhat 
with increasing temperature. Application of the 
instrument remains relatively simple if r is sensibly 
constant, even though the effect of changing y must 
be taken into account. It is to be noted that both 
a and r approach unity simultaneously, and that 
the dependence of @ upon y becomes vanishingly 
small as r approaches unity. 

To illustrate these points, assume that a sonic 
pyrometer, for which @ has been determined by 
using air having a y= 1.400, is applied in high- 
velocity streams of air at elevated temperatures, 
and that the latter temperatures are calculated by 
multiplying each observed value of 7; by the same 
constant a. Under the calibration conditions, a@ is 
simply the ration 7,/7,, of two measured quantities, 
and since the y of the calibrating air is 1.400, 


a=6/(5-+r). (8) 


This a@ is not exactly correct at any other tempera- 
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Fieure 1 Te m perature errors that arise if the variation of 2 


vith temperature ts neglected, as calculated for sonic purom 


eters having various values of recovery factor, 1 


ture, and from eq 7 and 8 the fractional error, /, of 
each observed value of 7, due to deviations in 7 
from the calibration value of 1.400 is 


} V7 1)(5y 7) 


- (9) 
(r-+-5)(y + 1) 


Such errors, expressed in percentage, are plotted 
in figure 1, A, as funcions of y at assumed constant 
values of r over the range of the latter, which is to 
be expected from a thermocouple junction in a sonic 
pyrometer. Using tabulated values of the y of air 
as a function of temperature [2], the percentage 
errors can be interpreted in terms of absolute error, 
as has been done in figure 1,B. It can be seen that 
changes in # due to changes in r and in y are small, 
and that it is desirable to make ras high as is practi- 
cable 

In the preceding example it has been assumed 
that the recovery factor of a given instrument has 
the same value, r=6/a—5, during the calibration 
and in actual use. If this assumption is correct, 
then the errors discussed in the foregoing example 
ean be avoided as follows: In a calibration under 
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convenient conditions (as with air at a relatively 
low temperature), determine r as 67,/7,—5; apply 
the instrument in any gas, and calculate 7’, from eq 
7, using the actual value of y and the experimental 
value One of the purposes of the present tests 
was to determine whether such a procedure was 
permissible with the sonic pvrometer 

It is to be expected that the recovery factor of a 
bare thermocouple junction located in a sonic stream 
and obstructing only a small fraction of the flow 
area will approach the theoretical value for a flat 
plate, which, according to Pohlhausen [3], and more 
recently to Emmons and Brainerd [4], approaches 
0.87. Another purpose of the present tests was to 
determine how closely this theoretical value could 
be approached in a practical sonic pyrometer 


of 7 


3. Apparatus and Results 
Three experimental sonic pyrometers, designated 
as 1, 2, and 3, were built and evaluated. Unit 1 
was designed for low temperatures only, with the 
idea that it would yield enough information to de- 
termine whether the construction of more elaborate 
instruments seemed worth while. Unit 2 was built 
for use up to 1,500° F, but without regard to its 
suitability for use in gas turbines, particularly as to 
its size. Unit 3 was designed ov a smaller scale, and 
its performance was investigated in much greater 
detail than was that of the other two 
The calibration equipment, with the No. 1 unit 
mounted for test, is shown diagrammatically in 
figure 2. The sole purpose of the pressure taps 
was for initial evaluation of nozzle performance and 
to insure that the flow was sonic at station 2 during 
use. The flow area at station | was so large that 
the impact effect was negligible, and the thermocouple 
there served for the accurate measurement of 7\,— 7), 
From this chamber the gas passed to the nozzle of 
the sonic pyrometer and was accelerated to sonic 
velocity at station 2 by maintaining the ratio 
P,/P,, always below the required critical value. 
This could be done in the test system either by 
compressing to increase ?,,, or by pumping to reduce 
P,,. Since it would complicate a service instrument 
to install a pressure tap at the throat of the nozzle, + 
second downstream pressure tap was installed at 
station 3. If the flow area between stations 2 and 
3 is made large as compared with the area of the 
throat, the pressures at these two stations are ap- 
proximately equal, and only the tap at station 3 1s 
needed. It is further advantageous to keep the 
pressure drop between stations 2 and 3 small in 
order that sonic flow may be maintained with a 
minimum of pressure at station | during discharge 
to the atmosphere or with a minimum of pumping 
on the discharge of the unit 
In tests at room temperature the heat loss between 
stations | and 2 is entirely negligible. Since air 
behaves nearly as a perfect gas, there is no change 
in total temperature between these stations, so that 
7, is identical with the measured 7;,. The calibra- 
tion value of @ is simply the ratio of the observed 
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the unit is 6/a—5 when the »y of the calibrating air 
is 1.400 


3.1. Experimental Sonic Pyrometer 1 


The first unit, designed for preliminary studies at 
temperatures below 400 
tion in figure 2, had a convergent-divergent 
inlet nozzle with a evlindrical throat three-eighths in 
long and 0.165 in. in diameter. A bare thermo- 
couple junction of No. 28 gage Chromel and Alumel 
wire was located centrally at the downstream end 
of the throat (station 2), directly alongside a static 
pressure tap. The junction itself was approxi- 
mately spherical and had a diameter about twice 
that of either thermocouple wire. An _ identical 
thermocouple was installed at station | in a com- 
paratively large inlet chamber, for the measurement 
of T,,=7T,,=—Ty,. The static pressures in the inlet 
chamber, nozzle throat, and exit were 
measured with mercury manometers 

Air supplied to the inlet chamber could be main- 
tained by electric heaters at any desired steady 
temperature in the range from room temperature to 
100° F. Thermal electromotive forces were meas- 
ured with a semiprecision potentiometer. Sonic 
flow through the nozzle could be maintained either 
by establishing the necessary pressure in the inlet 
chamber and discharging into the atmosphere, or 
by pumping at the exit when required 


F and shown in cross sec- 
brass 


passage 
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Ficure 3 Performance of sonic pyrometer 1 

The characteristics of the instrument, as deter- 
mined with atmospheric air (74° F; y=1.400) and 
with P;,=1 atm are shown in figure 3. Figure 3, A 
in which the pressure ratio P?,,/P,, is plotted agaist 
P,, in inches of mercury gage, shows that the nozzle 
of the pyrometer behaves in accordance with theory, 
since P,,/P,, attains a constant value of 0.53. At 
all values of P,, greater than 20 in. mercury gage, 
the flow is therefore sonic in the throat, and the 
desired condition for the operation of the pyrometer 
prevails. In figure 3, B, the observed difference ip 
temperature 7,,—T7,, is plotted against upstream 
static pressure. Again in accordance with theory, 
this difference increases with the velocity in the 
throat (i. e. with increasing P,,) until this velocity 
becomes sonic at P,,=20 in. mereury gage. For all 
higher values of P,,, the total and indicated tempera- 
tures of the gas in the throat of the nozzle remain 
constant and differ by 15.6 deg F. The measured 
values of 7,,=533.7° R and 7,=—518.1° R constitute 
a calibration of the unit. Inserting these values, and 
y= 1.400, in eq 7 gives a=1.030 and r=0.825 

In all experiments of this type, the thermocouples 
at stations | and 2 were made from adjacent lengths 
of wire to avoid possible errors from inhomogeneities 
The emfs of the thermocouples were always read 
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vidually and by bucking one against the other | four runs made at initial temperatures of 74° (cali- 
differentially), and in no case did the observed | bration), 279.0°, 395°, and 400° F, and of results cal 
erence in temperature depart more than 0.1 deg | culated on the aforementioned bases. The next to 


om that obtained by subtraction the last column gives values of the static temperature 


sing eq 7, the results may be computed in the | (7,,) in the nozzle throat, caleulated from eq 2 for 
wing Ways: (a) on the assumption that a always | the condition of sonic flow, that is, when \/= | 
the same value as that determined in the ecalibra- The tabulated values show that the sonie pyrom 


, run; (b) on the assumption that r always has the | eter gives values of the total temperature of the 
value as that determined in the calibration run, | high-velocity air that differ in no case from the value 
that @ varies with y;or (c) the appropriate values | observed before the air is accelerated by more than 
l,,, Ty, and y may be used in solving for a value | 0.6 deg F, regardless of the assumptions made as to 
for each run. Table 1 summarizes the data for | @ and y, and despite increasing differences between 


TABLE 1 Results obtained with sonic pyrometer 1 
>? 
Ru i / I , ’ 1 " P , 
\ hi ‘ un 7 i um . 
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k wn e 144.8 an 
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Mi s 
{ ' " wna Ls 
and y The calculated values of recovery factor ny parts were made of IS-S staink ss steel, and the 
sensibly constant, the maximum deviation from | nozzle was surrounded by an outer tube (1's in. 
mean value of r being only 0.2 percent outside diameter), which served as a radiation shield 


Although the range of temperatures used in the | To aid in bringing both the nozzle and its shield to 


tests described is very small, the tabulated results | the gas temperature, gas was brought to stagnation 


were believed to indicate that unit No. 1 showed suf- | in the space between these parts, and small bleed 
ficient) promise to warrant the construction and | holes were provided to renew this stagnant gas slowly 
evaluation of more practical units, designed for The throat of the nozzle had a diameter of 0.21 in 
higher temperatures and was cylindrical for a length of three-eighths in. 
A pressure tap in the throat was used only in pre- 
3.2. Experimental Sonic Pyrometer 2 liminary tests made to show that sonic flow was at 
tained at all pressure ratios below 0.53, in the man- 

In this unit, shown photographically in figure 4 | ner already illustrated in figure 3 
and schematically in figure 5, the nozzle and support- A No. 22 gage Chromel-Alumel thermocouple 





junction, filed to the form of a flat disk having a 
thickness equal to and a diameter twice that of a 
single wire, Was mounted centrally at the downstream 
end of the evlindrical throat. Thermal electromotive 
forces in the tests of sonic pyrometers 2 and 3 were 
measured with a Brown Electronik potentiometer 


‘ 


After calibration in the manner already described, 
the instrument was tested in the exhaust gas stream 
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from a Jumo 004 combustor, using the equipment 
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Fieure 4 Sonic pyrometer 2 Figure 5 Diagram of sonic pyromete , 














shown diagrammatically in figure 6. A’ shielded 
thermocouple, for which radiation corrections [5] 
were known accurately, and the sonic pyrometer, 
were inserted radially to the center of the 6-in. test 
section Cras velocities here did not exceed 300 


ft/sec, so that the corrections for impact on the 
thermocouple were always less than 3 deg F 
Observations were made at gas temperatures from 
100 1,300° F, in intervals of 100 deg F. In most 
runs, the temperature of the walls was made equal 
to that of the gas, so that radiation could not affect 
either the shielded thermocouple or the sonic pyrom 
to determine the effects of radiation 
runs were made with the walls 


to 


ete! In order 
on the latter, other 
cooler than the gas by measured amounts 

All of the results with sonic pyrometer 2 are sum 
The calibration data, obtained 
are given in the 


marized in table 2 
with air at about room temperature 


first horizontal line For this unit, the calibration 
gives a=1.0205 and r=—O0.8795 Starting at the 


left, columns in the table show in order the values of 
gas temperature indicated by the shielded thermo- 
couple (7\,), the corrections applied to the shielded 
couple for radiation (A7,) and impact (A7,), the 
corrected total temperature (7 r AT,+ AT, 
of the gas in the test svstem as determined with the 
radiation-shielded thermocouple, and the tempera- 


ture indicated by the thermocouple of the sonie 
pyvrometer (7: Observed quantities not listed in 
table 2 are the pressure in the test section and the 


pyrometer discharge pressures 
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Fiaure 6 Schematic diagram of haust-ygas susten 
a(obs) the ratio 7,,/7), is given for each run. Whilk 


the maximum spread in the observed values of this 


ratio is only 0.4 percent, the ratio shows a definite 
decrease as the gas temperature increases. As is 
evident from the last twe columns of table 2 
the ratio decreases linearly with temperature, in 
accordance with the empirical equation 7,,/7 

a (smoothed) = 1.0225—38 « 10°'7,,, in which 7), is 


expressed in deg KF. 

If the value a=1.0205, determined in the calibra 
tion with air at room temperature, is applied directly 
to the other runs, the temperatures listed in the 
column headed 4 are obtained The correspond 
ing errors, 7,,— 7), include both the experimental 
errors and those introduced by neglecting changes 


i 


ae) . . pe 
The remaining columns of table 2 show the effects | in y. Even so, the maximum error is less than 7 
of various possible methods of computing 7,, from | deg F, which is probably small enough for most 
the observed values of 7, In the column marked | purposes. 
TaBLe 2 Results with soni pyrometer 2 
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i 
\ more accurate method of applying the calibra 
m data is to ealeulate r for the calibration run, 
i 


IGURE 7 bD ’ pyromete 
ind to use this value with the appropriate values of 
This has been done 


neq7 in the column headed 

and the maximum error r is found 
» be half as great as in the previous method of 
duetion Furthermore r ip does not con- 
nue to increase with temperature. This merely 


hows that a small but definite improvement is made 


> taking the changes in y into account It also 
lustrates, once again, that the y need not be known 
precisely, since a change from H 100 to 1.326 alters 
the result by vee than 6 deg F 

Most of the differences, Tl. TY, are of the same 


sign, and in a direction indicating that the recovery 
factor determined in the calibration run was higher 
than that for the other runs. This can be examined 
detail by substituting the measured values 
and 7;,, along with appropriate values of 
to calculate an individual value of r for each 
As the tabulated results show, the maximum 
spread in r is about 2.5 percent; the calibration r is 
higher than most of the others and about 0.5 percent 
higher than the average. Possible explanations for 
this are discussed later in connection with the results 
obtained with pyrometer 
If the effects of radiation and conduction on the 
erformance of the pyrometer are significant, 
hese should be apparent in the last four horizontal 
lines of table 2. In these runs, the walls were from 
78 to 559 deg colder than the gas, vet the differences 
/ T;, and 7,,—TY are no greater than they are 
t the same gas temperatures in the runs without 
Hence it is concluded that the shielding 
plus the high mass flow rate over the 
easuring junction, combine to render the heat losses 
ntirely negligible. It is obviously a great advantage 
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hi ne S Dp 
of this tvpe of instrument that this feature can be 


built in, so that no separate determination or evalua 
tion of corrections is required 


On the basis of the highly encouraging results 
obtained with pyrometers | and 2, it was considered 
desirable to design and evaluate a smaller unit, 
Which might be installed in a restricted passage 
without blocking off the flow area unduly To this 
end, unit No. 3 was built within a evlindrieal tube 


of In addition 


cdles« ribe al 


having an external diameter 
to tests like the effects of 
the following features were investigated: axial loca 
tion of the measuring junction; shape of the junction; 


inh 


those already 


the technique of calibration, with particular atten 
tion to the effect of moisture condensation if un- 
heated and undried air is used; and the rate of re 


sponse to sudden changes in temperature 


3.3. Sonic Pyrometer 3 


In this unit, shown photographically in figure 7 
and diagrammatically in figure 8, the nozzle was 
made of 18-8 stainless steel and the outer tube of 

in. Inconel. The end of the latter was closed, 
and a in. diameter hole was drilled near the closed 


tube for admission of gas 

The nozzle, having a evlindrical throat 0.157 in. in 
diameter and \ in. long, was fitted snugly within the 
Inconel tube. except where notches were cul at either 
end admit the intervening space 
The tube was peened in to anchor the nozzle at 
both ends 


end of the 


to some to 


vas 


a. Effects of Axial Location of the Junction 
During tests of this unit 
Chrome! and Alumel were 


Provision was made for locating the 


various junctions of No 
mounted axially 
junction repro- 


o> 
vave 
22 gage 


duciblv at anv one of the following positions 
Distances 
from up 
Location Ne stream end 
of evlindri 
cal throat 
Inches 
l 0 
9 
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= 
th 














It was expected that sonic velocity would prevail 
at the downstream end of the cylindrical throat, but 
exact location of this plane in the hand-made nozzle 
is impossible. Hence were made with the 
junction at locations only one thirty- 
second in, apart in this region 

The calibration results obtained with air initially 
at about 100° F. are presented in figure 9, in which 
the observed temperature differences (7, Ty) are 
plotted against the suction, in inches of Hg, applied 
at the outlet of the pyrometer, for each of the seven 
junction The for locations 1 
through 3 become horizontal after the critical 
sure is reached The plane of soni velocity appears 
to be in the region between locations and 
probably between 3 and 4, since the curve for location 
4 continues to rise slightly 
The slope of the curves at higher suctions increases 
progressively at indicating 
supersonic flow there 

The curves of figure 9 indicate that the gas velocity 
increases progressively from location | to location 7 
It is sonic somewhere in the one-sixteenth in. between 
locations 3 and 5. Only at this place do the equa- 
tions for the sonic pyrometer apply exactly, because 
only here is the square of the gas velocity equal 
to ygRT,. However, this means merely that the 
junction should not be located downstream from 
position 5, and not that locating it farther upstream is 
seriously disadvantageous. As will be shown later, 
the instrument performs equally well with the 
junction at locations 2, 3, 4, and 5, and it is believed 
that location 1 would also have been found satis- 
factory if runs had been made with that configura- 
tion. This is so because the Mach number at any 
upstream station remains constant when the critical 
pressure ratio is maintained across the nezzle, and 
wcause the effect on recovery factor of variations in 
velocity at constant WW is insignificant 

Perhaps the most important conclusion to be 
drawn from figure 9 is that the junction should be 
kept at one and only one location during calibration 
and subsequent Full appreciation of the 
importance of junction location did not develop until 
after the pyrometer 2 were completed 
Therefore, some of the small variations in the ratio 
T,,/T» noted with that unit may have been due to 
insufficient care in bringing the junction always to 
the same axial location 


tests 


measuring 


locations curves 


pres- 
3 and 5, 
as the suction is increased 


locations 5, 6, and 7 


use 


tests of 


b. Performance Tests 


Tests of sonic pyrometer 3 in exhaust gas were 
made in the same manner and with the same equip- 
ment used for unit 2. A series of runs was made with 


the measuring junction at locations 2, 3, 4, and 5 


Again, most of the runs were made with the gas and 
the walls at the same temperature, and a limited 
number were made with colder walls to see whether 
the heat loss affected the instrument significantly 
All of the results are presented in table 3, which has 
the same general form as table 2, and therefore need 
not be described in detail 
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Ficure 9. Effect of arial location of the 


on the calibration of unit 3 Ts 


In general, unit 3 showed the same characteristics 
as No. 2, namely that it was capable of measuring 
accurate values of the total temperature of streaming 
exhaust that no for radiation 
conduction were required; and that the recovery 
factor determined with air at room temperature was 
slightly higher than the average of values determined 
at elevated temperatures 

In addition, the following characteristics are evi 
dent from the table. If an error of 15 deg F can be 
tolerated, it is necessary only to determine a@ with 
air at room temperature and to consider that a@ ts 
constant. The total temperature is then the con- 
stant @ times 7, and y need not be known. If the 
recovery factor is calculated from the calibration 
and this value is held constant for all other runs, the 
error is reduced to less than 3 deg F in the range up 
to 1,500° F when the appropriate values of y are 
used. These need not be known precisely. 

The error can be reduced to less than 2 deg F if 
the instrument is calibrated over the range of tem- 
peratures in which it is to be used, and the resulting 
values of ¢ are used in conjunction with the proper 
values of y. A separate accurate instrument is 
required for this kind of calibration, which is certainly 
not required for applications in engines where large 
temperature gradients are present normally. 
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Fictre 10 Effects of the shape of the measuring junction 


c. Effects of the Shape of the Junction 


The measuring junction used in obtaining the 
results shown in table 3 was approximately spherical 
This accounts for the fact that its recovery factor is 
lower than that of the flat junction used in unit 2. 
To further check on the effects of junction shape, 
unit 3 was calibrated with a junction in the form of a 
spherical bead, with one having the bead filed to a 
disk with thickness equal to the wire diameter, and 
with a third dressed until the upstream end had the 
area of the two thermocouple wires plus the two 
small welding fillets. These junction shapes can be 
visualized from the diagrams in figure 10, which also 
shows the results obtained with each of these june- 
tions at location 4. The recovery factor of the junc- 
tions increase in the order in which they have been 
described, so that the junction with most of the bead 
filed away is best for use in a sonic pyrometer 


d. Effect of Moisture Condensation on Recovery Factor 


As stated previously, each value of recovery factor 
determined with air at room temperature was slightly 
higher than the average of values determined at 
elevated temperatures. It was thought that this 
might be caused by condensation of water vapor 
during expansion through the nozzle, which is accom- 
panied by a decrease of about 100 deg F in the static 
temperature of the gas. To check this hypothesis, 
a calibration system was devised whereby the tem- 
perature of the calibration air could be increased at 
will from ambient to about the boiling point of 
water. Depending upon the relative humidity of 
the ambient air, somewhere in this range the tem- 
perature becomes high enough that the dew point 
will not be reached in the throat of the nozzle 

The results of this calibration, made with a junc- 
tion of minimum size at location 4 in pyrometer 3, 
are shown by the open circles in figure 11, in terms 
of recovery factor as a function of upstream tem- 
perature The filled circles were obtained in a series 
of runs to be described later The observ ed values 
of r decrease as the temperature of the inlet air rises 
and seem to approach a constant minimum value 
Effects of moisture condensation are believed re- 
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tures, showing the effect of moisture condensation in the no é 


sponsible for the initial marked decrease in 7 It is 
recognized also that changes in the axial location of 
the junction may result from differential expansion, 
though these are surely small Chromel, 
Alumel, and Inconel have nearly identical coefficients 
of thermal expansion 

Actually the curve of figure 11 is of more impor- 
tance in showing how small, rather than how large, 
are the changes in r The converging dashed lines 
at the right of the figure show the variation in r, 
which would produce an error of +0.2 deg F. As 
another example, at 400° F a change in r over the 
entire range covered by the curve, that is from 0.820 
to 0.833, would introduce an error of only 1.9 deg F 


because 


e. Effect of Free-Stream Velocity 


Although the free-stream velocity was not expected 
to affect the performance of a sonic pyrometer, this 
was checked in simple tests with available equip- 
ment. Compressed air, heated in an exchanger, was 
passed through a 4-in. tube and discharged through 
a 2-in. nozzle into the atmosphere. A total-tempera- 
ture probe was installed at the center of the 4-in 
tube just upstream of the nozzle. The %¢-in. gas 
inlet to pyrometer 3 was placed in the axis of the 
jet, about 2 in. from the nozzle. As in the tests just 
described, the measuring junction had a minimum 
area and was at location 4. The total temperatue 
of the air was varied from about 330° to 420° F, 
and its velocity at the pyrometer from about 260 
to 1,260 ft/sec. The results are presented in table 4 

The average value of the recovery factor over the 
range of these tests was 0.820, as can be seen from 
the filled circles in figure 11, as well as from table 4 
In the next to the last column of this table, it can be 
seen that the gas temperatures observed with the 
total temperature probe and those calculated from 
the sonic pyrometer by eq 7 with r=0.820 are 
practically identical. This agreement at all velocities 
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wh sonic is proof that the performance of the 
pyrometer is completely independent of the 
stream velocity 


{. Rate of Response 


Vhen a temperature-sensing Instrument immersed 
an engine is to be used for 
rate of sudden 


he gas stream of 


purposes, its response to 


rol 


nges In temperature is of greater Importance 
n its absolute accuracy lt was therefore of 
rest to determine the characteristic time T 
of sonic pyvrometer 3 Using a method that has 


this pyvrometer, with a 
was 


l n dese ribed elsewhere 6] 
erical measuring junction at 
ind to have a characteristic time of 0.65 see at 


location 5 


eral free-stream mass flow rates (@ between 
mad 5 Ib sec ft | ndet the conditions of these 
moausurements, the mass flow rate in the throat 


the sonic nozzle was approximately 30 Ib/see ft 
From previous tests, the bare junction alone in the 

stream was known to have a characteristic time 
ot L.1 see rate of 3 Ib/see ft Thus 

sonic pyrometer tested has a characteristic time 
six-tenths as great as that which 
iis measuring junction would have in the 
stream. This advantage is greatest at 
low stream velocities, and decreases to zero when 
the free-stream velocity is sonic, as may be 
In this example, it is assumed that a 
given thermocouple junction has a rt of 1.0 
at a mass flow rate of 4 lb/see ft®; that r is inversely 
proportional to G"? that the volume of gas within 
the pyrometer and upstream of the junction is zero; 
and that the total temperature and static pressure 
of the free stream are held constant at 2,000° R 
and | atm, respectively , as the stream velocity Is 


at a mass flow 


proximately 
alone 
me free 
seen 
in figure 12 


sec 


varied 

The observed advantage of the sonic pvrometet 
as to the rate of response Is not as great as might 
be expected from a consideration of the mass flow 
rates alone On this basis, the r of the junction used 
during the determinations of characteristic time 
should be about 0.4 see at a mass flow rate of 30 
Ib/sec ft Since the r of the complete instrument 
was 0.65 see, the difference must be due to the time 
required for the displacement of gas already within 
the pyrometer. It was therefore of interest to 
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determine whether the rv of the instrument could be 
materially by removing the stagnation 
chamber, that is, that portion of the Inconel tube 


that projected upstream from the nozzle 


decreased 


g. Effects of Removing the Stagnation Chamber 


Actually this modification reduced the character 
istic time of the pyrometer by only about 10 percent 
In addition, the unit became sensitive to the free 
stream velocity, as can be seen in figure 13 Here 
the apparent recovery factor of the modified unit, 
plotted in terms of stream velocity, Is 
decrease to about 0.38 at sonic velocity This cannot 
be a real recovery factor, and must reflect the effect 
of velocity It is concluded that the stagnation 
chamber performs an essential function in a practical 
instrument, and cannot be dispensed with for the 
sake of a sheht decrease in characteristic time 


seen to 
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4. Conclusion 


The so-called sonic pyrometer is an instrument in 
which a thermocouple junction is located in the axis 
of a critical-flow nozzle, at or near the plane of sonic 
flow. Adequate shielding from radiation is simple 
to accomplish, and the entire instrument can be 
built within a tube \ in. in outside diameter. As 
compared with a bare thermocouple in a free stream, 
the sonic pyrometer has the following advantages: 

(a) Errors due to heat loss by radiation can be 


made negligible under normal engine operating 
conditions 
(b) It measures free-stream total temperatures, 


regardless of free-stream velocity, so that the latter 
need not be known 

(c) It is relatively simple to calibrate, and accurate 
to within about 3 deg F up to 1,500° F when approxi- 
mate values of y are known and used in calculating 
total temperatures from indicated values. 

(d) Calibration with air at room temperature, and 
the application of a constant correction factor so 
determined, gives turbine exhaust gas temperatures 
that are accurate to within about 15 deg F up to 
1,500° F, without any knowledge of or corrections 
for the y of the gas. 

(e) Because the velocity over the measuring junc- 
tion is always sonic, the mass flow rate past the 
junction is greater than that of the free stream until 
the latter also becomes sonic. Hence a given june- 
tion in a sonic pyrometer responds more rapidly to 
sudden changes in temperature than would the same 
junction in the free stream. This advantage be- 
comes greater as the free-stream velocity decreases. 

As compared with a bare thermocouple in the free 
stream, the pyrometer has the following 
disadvantages: 

(a) It is more complicated and expensive to build. 

(b) A pump is required whenever the ratio of the 
pressure in the throat of the nozzle to the free-stream 
total pressure is greater than 0.53, and in this case 
the sonic pyrometer is not convenient for flight tests. 

(c) The bare thermocouple can be made smaller. 

The following precautions in the design and use of 
sonic pyrometers are suggested by the present results: 

(a) The measuring junction should consist of a 
bead of the smallest practicable dimensions. Such 
a junction has a higher recovery factor than any 
larger beads that were tried; one which, in fact, is 
within a few percent of the theoretical value for a 
flat plate. 

(b) The throat of the nozzle should be large 
enough, and the thermocouple junction small enough, 
that the blocked area is not excessive. In the present 
unit No. 3, the smallest and the largest junctions 
blocked about 6 and 11 percent of the area of the 
throat, respectively. Both configurations of the 
device functioned in accordance with theory. 


sonic 


(c) The measuring junction should be located 9 
near the plane of sonic flow as is practicable. Loca 
tion somewhat upstream of this plane is satisfactory 
but the junction should never be located downstrean 
in a region where the flow may be supersonic 

(d) The measuring junction should be fixed axial] 
so that its location cannot change during calibrati« 
and subsequent use 

(e) When rapid response to sudden changes 
temperature is desired, the volume within th 
pyrometer and upstream from the measuring jun 
tion should be kept to a practicable minimum 

(f) The resistance to flow between the dischare: 
of the nozzle and the external port of the unit should 
be kept low, in order that sonic velocity may }y 
attained in the nozzle with a mmimum of pumping 
effort. 

(g) Use of a convergent-divergent nozzle is also 
advantageous in reducing the pumping effort 

(h) The sonic pyrometer should never be used for 
measuring gas temperature unless it is known that 
critical flow prevails in the nozzle 

(i) Where high accuracy is required, it is recom 
mended that the sonic pyrometer be calibrated with 
air initially at or above 200° F. In this way errors 
due to condensation of water vapor during the 
expansion through the nozzle are avoided 

Characteristics of the sonic pyrometer that have 
not been investigated include the following 

(a) Its ability to withstand mechanical vibration 

(b) The minimum practicable size. 

(c) The maximum permissible throat area that 
may be blocked by the measuring junction. 

(d) The effect of lateral displacement of the 
measuring junction from the axis of the nozzle. 

(e) The effect of using dry air for calibration. 

Each of these items can be treated more logically 
if and when practical needs arise that cannot be met 
by the instruments and test methods described above. 
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Convergence of Cauchy-Riemann Sums to 
Cauchy-Riemann Integrals’ 


Otto Szasz and John Todd 


Iwo general theorems giving conditions to insure the truth of the relation 
lim 3 f (nh | f (x) da 
h- aa 
are established In addition, several cases involving Bessel functions are discussed, 
1. Introduction We note that (2) covers such eases as f(a re 
where f’(a 0 for r>1: that (3) covers such cases 
We are concerned with the validity of the equation | 8% ¢(@ r*, O0Ca<l (which is continuous except 
at z=0) and that (4) covers such cases as g(a J 
° F(xr)=sin*s 


lim h S} f(vk fix\dr, 


vhere the integral is taken in the Cauchy-Riemann 
nse, that is, 


are cases when (1) is meaningless if the sum 
ind the integral are interpreted in the ordinary way, 
it it becomes meaningful and correct if suitable 
immability methods are employed. A very simple 
isc Is when 


Chere 


ha sin 7, 


is true if both the series and the integral 
We do not pro- 
We shall also 


which (1) 
ure interpreted in the (C,1) sense 
pose to consider such cases at present 
concern ourselves mainly with the when the 
ntegrand has no singularity at the origin. Detailed 
nvestigations of the case when there is a singularity 
at the origin have been carried out by A. Wintner 
8}? and A. E. Ingham [3] 

lhe problem was studied by Bromwich and Hardy 

| in 1908, and some of the results they obtained are 
found in textbooks. Their results are of the 
following type: 


Cuse 


now 


If f(x) decreases steadily and if the integral (2) 
exists, then the series converges and (1) 
holds 

If v(x) is continuous and decreases steadily to (3) 


zero then (1) is true for f(z) 
If v(x) is continuous and decreases steadily to (4) 
; 


g(r) sin x. 


zero and if g(x)dz exists, if furthermore 


‘ i 
F(x) isa bounded continuous function, then 
1) is true for f(a g(r) F(x). 


Che preparation of this paper wa by the Office of Naval 
rch 


F igures 


ponsored (in part 


in brackets indicate the literature references at the end of this paper. 
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The re 2) has been generalized he aus 
The result (2) has been generalized to the quasi- 


monotone case by O. Szisz [6] as follows 
If for some a, for ] l 
have 


and 0 y 


and if the integral exists, then (1) holds 
Krishnan [9] has pointed out that in certain 
trivial namely, 


non 
cases a stronger result than (1) is true, 
hS* f(vh f(xr)\da t)) 
_—— 
for all sufficiently small positive A; indeed 
If f(x) is a function whose Fourier transform 7 
vanishes outside a finite interval [—a, a] 
and if the Poisson formula can be applied, 
then we have (for all @ 


hs Ilva 0 / 3 f(vh firjda 
provided O0<h<2ria 
A simple case when this applies is) when 
(r)=2*sin’tv, We use this technique to give an 
alternative proof of one of our results on Bessel 


functions 

In sections 2 and 3 we give two sets of conditions 
that are suflicient to insure the truth of (1) and that 
include many of the known cases. The remainder of 
the paper, w hich is independent of these two sections, 
various methods, of special 
Bessel functions, in which 


gives discussions, by 
mainly involving 
(1) is true 


CaSCS, 


2. xf(x) =y(x), a Periodic Function 


Our first theorem asserts the truth of (1 
f(x) =r 'W(2) under certam conditions on (4 
simplest case is that when (2) =sin z. 


where 
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Theorem 1 Let 
Wa >> 4 sin nz, » + b ‘ S) 


Let the Cauchy-Riemann integral 


ae / 
v dt (4) 
Jo t 
exist Then (1) holds for f(a x7 'W(r), 7. €., we have 
- (nh . l 
lim ,S wink WV") dz 10) 
a+ —_ nh J 
Proof We begin by establishing convergence of 


the series on the left of (10 We do this by first 
establishing that it is summable-A and then noticing 
that since ¥ is bounded, its general term is O(n7! 

and we can therefore use a standard Tauberian theo- 


rem [10, p. 501]. We therefore consider 
: (nh) ss 
5) eK Sb, sin wnh 
n=} n 7 n 
(] 1) 
> a > » s m 
b, int in vi h 
This interchange of summation being permissible 


because of the absolute convergence of the double 
series 


bh & b, (r"/n) 


S*b,! S\(r"/n) 
ee — * hm 


whenever O0<r<l. We require some ele- 
mentary properties of the inner series in (11). It 


is well known [11, p. 285] that 


now 


* . rsinz 
S*— sin nz—arctan 
ne} i rcos az 
when 0<r<1. This gives 
* oa r sin vh 
S*— sin yvnh—arctan 
ani nl l1—r cos vh 


when 0<r<1. Since for r=0, the left-hand side is 
zero, we must choose that determination of the arctan 


that vanishes for r=0. Since the denominator 
l1—r cos vk never vanishes (0<r<1h—-+0, v=—1, 
2, ...) we cannot leave this branch and, in par- 
ticular, we have 

ee sin vnh arctan ———— vi <he 

ani nh 1—r cos vh\— a 


Hence, we have 


= ¥(nh) r sin vh 
> > r*= >> b, arctan 
— a — l1—r cos vh 


where the series on the right is, for fixed h>0O, w 
formly and absolutely convergent in 0<r<1, bei 


majorized by 


pie 


}xS5\b 
’ l 


Hence, as r—1, we have [10, p. 339 
4 r sin vh , sin vi 
>> 5, arctan y > b, arctan 
| r cos vh / cos p/ 
>> 6, arctan (cot }vh 


This means that >) ¥nh)/n is, for all A>O, sun 


’ 


mable—A, and therefore convergent, to the sum 


>> 5, arctan (cot byh 


We shall now show that, as h—0, 


‘ ] l . 
pa arctan ( cot 3 wh ) »— «> b, 


l - - 


The difference between these two expressions can 
be written as 


St ee Pe Wee eB 
>) > arctan( cot 5 vh ) g FOr, 3 
> > = } ) 
2, arctan ( cot wh )- sr >, b, 13) 
’ N+1 - - N+1 
oo ; l . 
We first choose N so large that Pp > b, (and there- 
2 NF 


b, arctan ( cot z vh ) is small. With 


1 - 


fore also 3 
\ 


N fixed we then can choose A so small that the first 
sum in (13) is small. This establishes (12). To 
complete the proof of Theorem 1 we have first to 
establish the existence of the integral on the right of 


im . 
(10) and then show that its value is 57>) 5}, 


Our hypothesis (9) insures convergence at the lower 


limit. Integrating by parts, we see that 
*" W(x) (u) "uy, (7) 
[PAP ar Oy (Mae, (14 
/ t u / a 
where 


°u 
¥,(u) V(x)dz. 
fF 


Since ¥,(z)=O(1) by periodicity, and since] x~*dz 


is convergent, we can let u-> @ in (14) to obtain 


a Wr) 9, [" vil) 9, 
/ z * 


” z 
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Thus the existence of 


& 
ad 
J 
established 
We now consider 
. y . >> 5 SIn vs 
J 
dx— | da 
J I ‘ 15 
. Sil 
Sh ls - 
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We 


This inversion may be justified us follows 
he fact that 


"= sin ft cos ft "= cos tdt 
Nit dt - Wy 
; t t J t 


use 


SJ What we have to show is that 
*7 b, sin va ae 
/ » > da - 2 b, sin vr da 
, \ j I \ 
nad 
\ . 
. Sin wa 
cnet d 
"|, : 
both small for N, 7 sufficiently large. Since 
S*b, is absolutely convergent, and since for all », 
Sin ws 


O(T~'), it follows that /,=O(T™'), i.e 


d 
small if 7 is large uniformly for all NV. Fixing T we 
an choose NV so large that J, is small. This estab- 
lishes (15 
We can let ¢t—-0 in (15 


I vht Is uniformly convergent in t 
sm vi l 


| V2) re — 3 b ds 
Jo J ant ‘ a 2 


since the series on the 


This gives 
So, 


and thereby completes the proof of Theorem 1. 
Remark: An example that is not covered by this 


theorem is the case when 6,—(—1)"/n. Here 
y(t) = >°*b, sin nt al t+ . ! 
— { 2 j 
ola 
+, Is not an integer and ¥(f)=0 if t+, is an in- 


teger, and although it is clear that f(x)dz exists, 


the question of the truth or falsity of (1) seems to be 
rather deep. 


3. xf(x) =y¥(x), a Fourier Integral 


Our second theorem asserts the truth of (1) when 
r) is the sine transform of an absolutely integrable 
netion n(t). This, when we take 
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reduces to the case 


Ae fe Zz. 


which has been mentioned by Hardy and Bromwich 


[1] and by Krishnan [9] 


Theorem 2 Let 
ic 
AG n(t) sin rtdz, 
t. 
where 
n(t)\dt 
is finite Then (1) is true 
Proof We have 
. Jf 
hS3 fivh > > sin vit nit dt 
. Sh | 
24 sin vht n(t)dt 
’ v 


This interchange of summation and 
= SG Ta. 
legitimate because >): 


integration 1s 
' sin vAt is convergent and its 
partial sums are uniformly bounded 


It is well known [10, p. 523] that for z positive and 
not an integer, 
>) v-' sin 2rva x{x—|[z]— 3 


1 


where [x] denotes the integral part of x. This gives 


“ee ( At ht | 
S* py! sin vht r— t. 
om (2Qq7r Jr 25 
Substituting, we find 


>> vy sin vht¢ nit dt 


. (At ~ At 1) ; 
|, (2r Espen 
"= At ht 
e t n(t)dt. 
Jo (2x 2n_]}) 


We shall show that the second term on the right tends 
to zero with h To do this, write 


">< ae At } 
t)d 
, \2e Esk , 
“(ht At) Po ¢ A ht) 
5 [5 |i nwaes [ 1D ae cade 


I, (h), 


n(t)dt T 








I, (h) t 


We can choose and fix w so large that 


say. 








I, (h) 


is small for all A 


that 


Hence 


AsS* / vh) 


—_ 


To complete the proof “ 


. 


f(r)da 


Then we 


At 
ar 


is itself small enough to insure that |/,(h 


dt 


can choose h so small 
At At 
Jr Jr 
is small. 


| n(tydt. 


e have to show that 


n | n(tjdt: 


but this follows easily, since 


[, Sorde— Jr J, 


ime 


iJ, 


sin at 
n(t)dt dr 
' 
“? sin rt 
: dr | ” t\dt 
0 } i 


n ‘dt, 


the change of order of integration being permissible, 


since n t)dt 1s, by 


gent 


4. The Cases f(x) 


It is known [5, p. 


DS Junk 
and so, as h--0 

i ty 
Since also !2 p 96) 

Ja 


the relation (1) is 


hypothesis 


Jo(x) and f(x) 


36) that 


absolutely conver- 


h (16) 
(nh l 
idx=1, 


true in this case 


The result (16) can be established directly in the | 


following manner. From 


the asymptotic formula 


using partial summation, it can be shown that t! 


yy 
SeTies 
ri 

( 


J (nh) converges for 0< h mr (and, indeed 


converges uniformly in every closed subinterval! 
We evaluate the sum of this series by determining 





git 
((’,1) sum, using the integral representation 
> *s 
Ja cos (Fr cos tf) dt. 4 
Ji 
and the facts that 
l 
; sin (n+, )2 
(,(4r) S* cos pa t ’ i 
ca 2 l 
2sm og 
l 
sin (x t ) J 
l Z. l 2 
S ( ] , ‘) 
n : J 2 l 
2(n 1) sin I 
» 
From (17), (18), and (19) we obtain 
. 7 sin n | h cost] 
' 2 (*'*}] 2 
s,(h >) do(vh oT | dt 
ue 0 os 
: | 2 sin (5 A cos t) 
I . 
a, (h) S sith 
i _— 


so [erytond 


(n 1) sin {: h cos t)) 





Now sin@*u l/a 
h->O we find 


7/2 810 | 3 ( n-+1)h sin¢ | 


| 
wii Lh? Jo i dt+OV. 


as s->0, and so, for fixed 


sin 


l 
oy R,(n) R,(n) Un~*), 


where FR, and R, are the contributions from th 
ranges (0, #/4) and (#/4, 2/2), respectively. Clearly 


i 
R,(n)- O(n"). 


rnh? sin 


To deal with R, we change the variable ¢ to w 


l . es 
, (n+ 1)A sin t; we find, writing A(n+-1)2 k, that 














J (x) (= Jésin (2 } ; r)+ O(r-*/?), 
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? sin? udu 
rin 2. | | 
rh. \ wll ; 
. ( ( n 1)/ \ 
2 (*/sinwu 
rh ( u )é adi 
2 (* sin u \| Ju ) 
| di 
ad, | u ya ( n l i) | a 
k—>o@ and the first term tends to A-' and 
ndeed A-'4-O0 
Let 2u/( 1)/ then in the range O0<u </ 
have O<97 and 
| 
| 
| | 
| 
, Su 
1 / | h, 
This means that the second term is OG “as 


Gathering our results together, we have 


l 
ath +. fj Ain 


This means that the series $)J,(nh) is summable 
f | to the sum } h when h -0 This com 


letes the sketch of the prool of (16) 
lo deal with the case f(a Ky(r2) we use the for 
ila + & p 301] 


: 7 . | | l 
S* Ky(nh) —r4 r - te x4 
Zh yh in2h? 2ne é 
l | 


. Inh log 4x, 


whi h leads to 


limh S) Ag(nh)=, x 
h—0 yA 
Since |12, p. 388] we have 
. | 
K (r)da . 


is true in this case. We note that here we have a 
logarithmic) singularity at the origin, and have 
sed a summation starting at 1 and not at 0 to avoid 
ouble 


5. The Cases f(x) = J,(x) cos xt and f(x) 
J.(x) sin xt 


It is known [4, p. 59] that for A>0, 0 <t- 





6. The Case f(x) 


1, we have 
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S* J.(nh) cos nht 


—_—_ " 


Multiplying by A and then letting / 
O<t< J 


On the other hand |4, p. 51}, we 


and so (1) is true 


The cise when Wi considel 
f(a Ji (xr) sin at 
is shghtly more complicated It is known 
that for h>0. O<te< 1. we have 
>a nh) sin (nhl 
| 
S ri th } Ql 
Gum | 
S Nanth iD 
where 
I rl th / J 9 th 
ar th h rl 


Now for small A and fixed ¢ 

nth h 
» O( = ) 
Hence as hh 0) 


h>S° J inh) sin (nht 
On the other hand, it is known [4 
i) { l 
J Ax) sin rt da 
Thus (1) holds in this case also 


A, (x) 


bu , 


It is known [2, p. 91] that 


wr (yp 


h r( n+) 


> A (nh ; 





(2’x ")E(u 





and so 


The convergence of this series can be established in a 


nanner 
12, p. 391] on 


similar 


Since also 


the relation (1 
zation of that 
when lu 


discussed 


is true in this case 
ith 
. reduces to the case f(a 


to that used in the case u 0 


IS2 


and Is a venerall- 
This 


sin sod 


section 4 


cause, 


We shall now discuss this case using the Poisson 


summation formula 


yi ia A 0) + SSF (np 


where a8 =22, a >0 
of f defined by 


This formula is true, for 
provided that f(2) is 
Variation in OSs . 


that the Cauchy-Riemann integral 


and Ff 


continuous 
that f(a 4) 


\! . / 
S , 
Va GI an.) (na), 
> 1) 
is the cosine transform 
f(t) cos rt dt 


60 to 64). 
bounded 
and 


instance [7, p 
of 


as d 


and 


j tydt exists 


These conditions are satislied when «> } if 


») 
{(z)=-— | ! <3 
r(ut+s;) 
fir)=0 r>t1 
so that [7, p. 65] 
J a 2 
F (2) . Al 
ge” I'( yu l 
In this case (21) gives 
{1 . } 
+ >) A, (np) 
YT (uy a £S \ 


This formula is true as it stands for w >}; if w= 4 
a term on the right corresponding to n= a 
it must be halved; if —}<yp<} 
true, provided that a@ is 
integer 

We now assume a>! 
Is empty; 


iMac 
the formula remain 
not the reciprocal of ay 
so that the sum on the right 
this means that 


In these circumstances (22) gives 
, 
l sly ‘4 (ng)! (=e) l 
2*T (yu 1)(2 =— ) re | 
2 yw ) 
, 
that is 
~ \ nB \ 7 | fi | | 
ry p l 2 
. I (4 T ) 
> 
which is (20 
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Maximum Likelihood Estimates of Position Derived From 


Measurements Performed by Hyperbolic Instruments’ 


Eugene Lukacs 


Certau electronic surveying tostruments referred to a hyperbolic Instrument 
determine the difference between the distances from two fixed stations to a moving object 
Data obtained from a single observation vield one relation between the three coordinates of 
the movu object and restrict it to a hvyperboloid of revolution lo locate a position it 
space, observatior from at least three pairs of stations are needed If more tha three 
observations are taken, more than three equations for the three coordinates are obtained 
In general these equations will not be compatible; it is then a statistical problem to find ar 
estimate for the Inknown position In this paper maximum likelihood estimate ure 
btained for positions it space derived from such hyperbolic data 
. . 4 > 
l. Introduction coordinates of / \ single observation based on 
the stations S, and Sy gives then a determination of 
sSevera electron surveying instruments mve Tran TS 
the distance PS, — PS, If NV stations are available 


en developed in recent years. These instruments 


effect 


tilize the Doppler and are designed to 
ate the position of a moving object. In- this 
iper we are considering measurements obtained 
y using certain types of these instruments, which 


re often referred to as hyperbolic instruments 
Hyperbolic instruments, such as Raydist (see foot 
the difference the 


Stations to a moving object 


4 determine between 
from two 
obtainable from using 
vield one relation among the three 
the To 
space hyperbolic data nt 
If more than three 
more than three equations 
In veneral 
hese equations will not be compatible; it is then a 
find an estimate for the 
The use of exactly three obser 
vations is justified if it can be assumed that the 
not affected by any 
In general, this assumption is not realistic and should 
not be made without a thorough investigation. It 
appears desirable, therefore, to use more than the 
minimum number of observations and to estimate 
he location of the moving object by statistical 
The purpose of this paper is to obtain 
maximum likelihood estimates for positions in space 
derived from hyperbolic data 


ote 
stances 
Data 


mir of 


one observation one 
stations 
object locate a 


oordinates of moving 


wsition In from least 


observations ure needed 
»bservations are taken 
for the three coordinates are obtained 
tatistical problem lo 
nknown position 


data are observational errors 


methods 


2. Formulation of the Assumptions and the 


Problem 


We denote by N>4 the number of stations 
Sala=1,2, N) and by (d4).@e2,4a;) the known 
coordinates of the station S,. Let P be the unknown 
position of the moving object and (@,, %, 0) be the 


sifled Navord 
Naval 


An unclassified treatment of the mathematical portion of a cla 
t written by the author while a member of the staff of the | 

nance Test Station, Inyokern, Calif 

J. F. MeAllister, Measuring velocity of 


V-2 rockets by Doppler effect, 


Pech. 6, No, 2, p. 56 to 59 and 129 (1947 
Hastings, 
lech, 6, No. 6, p 


Raydist-——a radio navigation and tracking system 


#0) to 33 and 100 to 108 (June 1047 
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it is possible to make at most 


(\) . 


based on different pairs of stations 


llowever, not all these observations will be used in 


observations 


practice Let ns{NCN 1)}/2 be the number of 
observations actually taken We make the follow 
Ing assumptions 

1) The hyperbolic measurements 2,9 are observa 


tions on independent random variables Z,, The 
variates Zs have rectangular distributions with mean 
ua and common range 2r. Further Z,, Zaa und 


Il) The parameters of the n rectangular distribu 
tions may be expressed in terms of the coordinates 
of the unknown point by means of the formula 


k(Z, Me yu! a. Vy a 1) 


(II The distances d., and d, from the position of 
the moving object at time f—0 to the stations S, 
and Sys, respectively, satisfy the inequality 
d, d / 2) 
The nature of the problem as well as its solution 
are to a great extent determined by assumption (1 
To justify this assumption it is necessary to describe 
briefly how hyperbolic data are obtained 

It is desired to determine the instantaneous loca- 
tion of a moving object. A transmitter is carried 
in the object and emits a continuous signal of fre 
queney fo. A transmitter (the reference 
transmitter) emits a signal of a frequency f’o, slightly 
different from f). Each station receives both signals 
simultaneously. A certain beat 


second 


frequency nq 1s 











Ne depends on 


thus associated with each station S,; 
the velocity component of the object away from S, 
and on the frequencies f, and /’ To obtain an 
observation two stations S, and Ss have to be linked 
by wire by radio The frequencies n, and ng 
give then a beat frequency; the quantity recorded ts 
the number of eveles n,» of the beat frequency signal 
occurring during the time of observation It can be 
assumed that the observation 2,5 on Z,»9 can be satis- 
factorily approximated by a linear function of nag 
given by 


In this formula ¢ is the speed of light. The frequeney 
Nag may be recorded continuously to permit reading 
nas to fractions of aecyele. Alternatively it is possible 
to devise a counter that registers full cveles only.* 
In either case it possible to determine a 
range for the reading error of n.» (and therefore also 
of Moreover, the recording devices are of such 
a nature that any error within the possible range is 
equally likely to occur. It therefore, quite 
realistic to assume that the errors have a rectangular 
distribution about zero with a certain range. This, 
however, is equivalent to assumption (1 Equation 
(2) is a condition imposed on the location of the 
stations. As the half range r is rather small (2) 
implies hardly more than d,>d,. This considera- 
tion shows that (2) is not a very serious restriction 
The aim of this paper the determination of a 
maximum likelihood estimate of the unknown posi- 
tion. That the unknown coordinates 6,, 4, 4s 
have to be determined so that the greatest possible 
probability is given to the actual observations. This 
discussion will be based on assumptions (1), (11), 


and (IIT). 


is easily 


seems, 


1s 


1s, 


3. The Maximum Likelihood Estimate 


Let fas(z) be the frequency function of the random 
variable Zs. According to assumption (1) 


Rs 
, fi if pas PS2SKesT!. 
haa(z) <a? 


bo otherwise 
The likelihood function® assumes, therefore, also 
only two values, the value (1/2r)" and the value zero. 
If one wants to determine the unknown parameters 
6, 0, 0 so that the greatest possible probability is 
given to the actual observations, one has only to 
chose the parameters so that the likelihood function 
is not zero. This is the case if and only if each 
observation 2.3 falls within the range of Z,s. 
That is, the likelihood condition 
* Hastings Bulletin R-23, Complete Raydist tracking system for missiles (Sept 
vA description of the maximum likelihood method and a description of the 


likelihood function may be found in, H. Cramér, Mathematical Methods of 
Statistics, p. 40 (Princeton Univ. Press, Princeton, 1946 








Bea? 8=>Mas tl, 4) 
must be satisfied for the n pairs (@,8) used From 
(1) it is seen that the likelihood condition (4) consists 
of a set of inequalities involving 6, 6, @ A point 


A. As, @) of the parameter space 1s said to be a 
maximum likelihood estimate if its coordinates 
satisfy the maximum likelihood condition. In gen 
eral, mequalities (4) will not have a unique solution 
Kither all the pomts of a certain region in the para 
meter space will satisfy (4), or it may happen that 


no maximum likelihood estimate exists. If maxi 
mum likelihood estimates exist, one may consider 
the set of all the maximum likelihood estimates 


This set will be called the maximum likelihood region 
Any point in the interior of the maximum likelihood 
region may be used as an estimate for the unknown 
position. The most important problem is the deter- 
mination of the maximum likelihood region. This 
problem will be solved by expressing the boundaries 
of this region in terms of the actual observations 
This requires certain transformations of the likeli 
hood condition (4). If we substitute the value of 
dos from (1) and transpose we obtain 


32—Tr >> (0 Ges)" S*(0 a . ? 
Ve Ve 
fa) 
We next introduce the following notation 
ho 
n (d. ds r) 
, 
j th) 
ny —~" (d,—ds-+r) 
e 
From (4a) we see ny=ny—2rfp/e< no. 


In the following we distinguish three cases, namely, 


Case (A): Nas >No 
Case (B): Nag Ny 
Case (C’): N; S Nas S No. 


From (3) it is seen that 


243—r >0 and a fortiori z.3+r>0 in case (A) 


Zas+r< 0 and a fortiori z.3—r< 0 in case (B) (5) 


Zap +r OS 2ag—F in case (C). 

We first consider the conditions imposed by a single 
observation 2.3 based on an arbitrary pair of stations. 
It is no loss in generality if we assume a=1 and 8=2 
To simplify the discussion we choose an appropriate 
system of coordinates. Let the position of the mov- 
ing object at time t=0 be the origin O and choose the 
6,-axis so as to be parallel to the line joining S; and 
S, the @,-axis should be in the plane OS,S, and per- 
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») 




















ndicular to S,S, The @-axis is perpendicular to 
6,0-plane; finally some orientation is given to 
axes With this choice of coordinates, we have 


a = il 
a a b 
a a 0 

Wi shall write a a,.a ado: a a. hb 


lo make the formulae more concise, we introduce 
following notation 


; / \é 
2 } 2 } (7a 
R PS \ 0 a 0 b t Oy for i 1,2 7b) 
m a a 2; if 
d a a.)/2: d* a a,)/2 7d 
] Me id ry m ye. ft 
1 m td v7 m 4d" (71 
2d 2d* 
i=- ii=<s (7g) 
2d 2d* 
1, is * (7h) 


Considering that according to assumption (1) 
, we obtain from (7) and (7a) 


* =—~¢@ 


= and 2 (8) 
We see also 
d* d, (Sa) 
and therefore 
Ti+ 7,=T7,+T73=2m, (Sb) 
from these relations we see 
AY A, and A} A\. (St) 
\ simple computation gives RY —R,? 4d(0,—m) 


so that we have with the aid of (7e) and (7f) 
If the origin is located on the line S;S: the three points OS;S; do not determine 


plane. In this case the @-axis may be chosen in an arbitrary plane so as to 
rsect the 6,-axis perpendicularly 








R?— R3—2*= 4d (T,—6,) ) 

R?—R id (T,—¢ 
R?—R?—2?=4d(6,—T* ' 
R?—R?—F*=4d(0,—T?) J 


We next rewrite the maximum likelihood condition 
for each of the three cases (A B) and (C 
Case (AA) 


From (5), (7) and (7a) it is seen that zg r, 
z r. It follows then from (4a) that 0< 2 
R,SR R,. By squaring and transposing we 


obtain the inequalities 


R R R R 2 
R,: = ~ and R,> ~a 
Considering (9), (72), and (7h), we obtain from these 


relations the final form for the maximum likelihood 
conditions 


(A.1) R.s AAT,—96,) and (A.2) R,= A\(7,—0é 


In an analogous manner we obtain the likelihood 
conditions in Case (B) 


(B.1) R,< At(T%—6,) and (B.2) R,> ANT? —0 


It is seen that the inequalities (B.1) and (B.2) may 
be obtained from (A.1) and (A.2) by writing A, in- 
stead of R, and by attaching an asterisk to the A 
and 7’ symbols 

Case (C) Similarly, we obtain the likelihood condi- 
tions 


(C.1) R= A,(7,—94,) and (C.2) R,= ANMTT—¢é 
It is seen that (C.1) A.2) and (C.2) =(B.2) 
4. Discussion of the Likelihood Condition 


Before discussing the likelihood condition, it is 
desirable to derive a simple relation 


From (2) we have d,>d,.; in our coordinate system 


this means 


From this we see easily 


sign d=sign (a,—a sign a, sin d*. (10) 
Here sinh J stands for “the sign of 2’, that ts, 
sign J Lif s >0 and sign z Lif r<0 


We first discuss the case (A) and start with condi- 
tion (A.1). It is seen immediately that (A.1) can 
only be satisfied if A,(7' A)=0. 
Considering (7h) and (10) this necessary condition 
can be written 

T, sign a,= 6 sign a, (11) 
This condition restricts the pom|l P to a half space 
determined by the observations and the location of 


the stations. In the following we assume that (11) is 
satisfied. Wecan, therefore, square (A.1) and obtain 
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Ris AX T,— 6 


it follows that 


(0,.—b R<(A 1)@— 2(A?27.—a.) 0 127 a 
(12) 

We have to distinguish two possibilities (a) If 

1 1 then from (7h) and (7f 7T7.=m—d=a 


The right-hand side of (12) is then zero so that (12) can 


only be satisfied for the points of the line @,= 6, 
6,—0 
We consider next the case where (b) A?#1 From 
(7f) and (7h) it is seen easily that 
AiT.—a m(A l (12a) 


Using this formula we obtain by a simple computa- 
tion from (12 


(12b) 


The equality sign in (12b) determines a quadric 
surface of revolution Q,. The center of this quadrice 
is the point (m,6,0); the axis of revolution is parallel 
to the @,-axis and is therefore the line joining the two 
stations S, and S,. The surface Q, is a hyperboloid 
if A?>1 and an ellipsoid if Aj< 1 Denote by u 
the real (resp. major) semiaxis and by ¢ the imaginary 
Then it is seen from 


(resp. minor) semiaxis of @ 
(12b), (7h), and (7f) that 
u?—=- and v?7== | A | (13) 
} 
Let 2e be the distance between the two focal 


points of Q,, then e?=u* + r’sign(.Aj— 1) or from (13), 


e=d and u=2/2 


e=d and u (13'a) 


The focal points of the quadric Q, are therefore the 
stations S, and S,. The quadric surface Q, divides 
the space into two regions, we call the region con- 
taining the center the interior of Q, and its comple- 
ment the exterior of Q,. We can then say that 
condition (A.1) is satisfied outside Q, if A?>1 and 
inside Q, if A?<c1. If A?=1 then (A.1) is satisfied 
only on points of the line @,—6, @,—0 for which 
A,=sign(7,—9,). 
We next consider condition (A.2). This ine- 
quality is satisfied in two different regions. In the 
first region the right-hand side of (A.2) is positive, 
in the second the right-hand side of (A.2) is negative 
or zero. We assume first that the right-hand side 
of (A.2) is positive. By squaring and transposing 
we obtain then 
2(A?7T, —a,)0, + A?T?—a?. 
(14) 


(0.— b)?+ @= (A?—1)# 





(a) If A?=—1] 
then from (7e) and (7g) 7;=m-d=a,.. The right- 
hand side of (14) is then zero; in this case (14) js 
satisfied for all points; (b) if Aj#1 we see 
from (7e) and (7g) that 


Again we distinguish two possibilities 


easily 


Ai?T,— a,=m(A?—1) (14a 
Using this formula (14) reduces by a simple compu- 


tation to 


A?—1)(0 (0.—b #5: 14b 


my 


The equality sign in (14b) determines again a 
quadric surface of revolution Q.. The center of this 
quadric is the point (m, 6,0); the axis of revolution 
is again the line joining the stations S; and S,. The 


| surface (), is a hyperboloid if A,*>1 and an ellipsoid 


if A’<1. Denote by p the real (major) semiaxis 
and by g the imaginary (minor) semiaxis of (, 
Then it is seen from (14b), (7g) and (7e) that 
, and q? Ai—1 (15 
/ 4 il { 1 4 


Let 2f be the distance between the two focal points 
of @. By an argument similar to the one used 
before, we obtain 


(15a) 


f d and Pp 


so that the quadrics Q, and @, are confocal. We 
summarise these results in the following statement 
If the right-hand side of (A.2) is positive, the 
condition (A.2) is satisfied inside the quadrie (, 
when A,?>1, but outside Q@ when A,’<1. In 
case A,’=1, condition (A.2) is satisfied for all points 
for which A,=sign (7,— 4). 
We investigate next the second region in which 
(A.2) is valid by assuming that the right-hand side 
of (A.2) is nonpositive; that is 


A,(T,—-6,) s 0. (15b 
From (7g) and (10) we see that sign A,;=sign a, so 
that (A.2) is satisfied if 


7, sign a; @, sign ay. (15e 
From (15c) and the necessary condition (11) it is 
seen that (A.2) is satisfied if 6, belongs to the interval 
(7;, T2) in ease (15b) holds. Therefore, (A.2) is 
satisfied either in the above-mentioned region deter- 
mined by (» respectively in a half space or in a region 
bounded by two planes perpendicular to the line 
joining S,; and S). 

In case (A) we have 2=2,.—r>0 so that 2>2. 
Therefore, we see from (7g) and (7h) that A,’ >A; 
and using also (lla) and (15a), up. We are now 
ready to determine the region where (A.1) and (A.2) 
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satisfied simultaneously We have to distinguish 
possibilities 
\ A? >1. It is easily seen that the point 

5.0) [resp. the point (7),6,0] is on the line seg- 
nt joing the vertices of , (resp (). From 

previous results we see that (A.1) and (A.2) are 
tistied inh one-half of the hyperbolic shell bounded 

t/ and (/o. 

. Af >I=A The point (7',6,0) coincides 

th S.. while T,,6,0) is on the line segment jolning 

Conditions (A.1) and (A.2) are 
tisfied in one-half of the outside of the hyper- 
nd @ 

\ 1>A The point (7;,6,0) is then out- 
de ©, and outside Q, while (7),6,0) is inside © 
‘he region consists again of one-half of the outside 

the hyperboloid ( 

1 A l>A; In this case conditions (A.1) and 
\.2) are satisfied only by points of the line S,S,, 
hich obey 11 
5. If 1 > A> A? no region exists where (A.1) and 
\.2) are simultaneously satisfied It is worth while 

remark that cases (2), (3), (4), and (5) are of 
ittle practical interest Krom (72), (7h), and (10 
is seen that A? S1 (respectively A,’ S1) is equiv- 
ent to 


e vertices of / 


Id | « ; respectively (2d) < } 
\ situation where these relations hold is conceiv 
ble but unlikely as long as ris small. This is seen 
f we remember that the difference between two 
sides of a triangle is always less than its third side 

Case (B) can be discussed independently by repeat 
ing the arguments used for case (A) 

We can then summarize the results by distinguish 
¢ the following five cases 

l. A,*>A; 1. In this case (B.1) and (B.2) are 
satisfied in one half of the hyperbolic shell bounded 
by (, and @Q, 

conditions (B.1) and (B.2) are 
2. Ay ol=As atisfied in one half of the out- 
OA 1>A? satisfied in one half of the ou 
side of Qs. 

1, A 1 >A”. Then (B.1) and (B.2) are satis- 

fied only by points of the line S,S, that obey (11) 


5. If 1>A?>A’ no region exists where both 
B.1) and (B.2) are simultaneously satisfied. Cases 
2 3), (4), and (5) offer again little practical 


interest 

We finally consider Case (C) with the conditions 
C.] (A.2) R,2A,(7,—6 and (C.2)=(B.2) 
R,=A,*(7T;*—6,). It is seen that (C.1) is identical 
with (A.2) and (C.2) with (B.2), therefore, the pre- 
vious results can be applied immediately. It is 
advisable to remark that no necessary conditions 

milar to (11)] exist; however, (C.1) as well as (C.2) 
are satisfied if the right-hand sides are nonpositive. 
Considering the previous results, case (C) may be 





summarized by distinguishing the following five cases. 

l. A A’ ] lor A’ A 1] Conditions ((" 1) 
and (C.2) are satisfied in a simply connected region 
bounded by one nappe of Q, and by one nappe of Q, 
so that the convex side of one nappe is turned toward 
S,, while the convex side of the other nappe is 
turned toward So 

2. A} >Aj=1 [or A} >A?—1]. The conditions 
(C.1) and (C.2) are satisfied in the region that is 
bounded by one nappe of one of the quadries and 
that contains the center 

3. AZ >1 >A} lor A} >1>A}]. The conditions are 
again satisfied in the region that is bounded by one 
nappe of one of the quadrics and that contains the 
center 


4. Aj=1 
§. 1>A! 


A? lor A? | 43) { the es sang = 
Ajlor1>4i>4i]| fa 

In case (C) the observation gives only little infor 
mation about the location of the moving object 
However, case (C) can occur only if al<c7 This 
will happen only for a small portion of the path of 
the moving object ifr is small. Moreover, it should 
be possible to locate the stations in such a way that 
case (C) occurs only for one observation at a time so 
that this possibility should not cause any difficulty 
in locating the moving object 


5. The Maximum Likelihood Region 


In section 2, it was assumed that n hyperboli 


observations (based on n pairs of stations) were 
taken to determine the position of the moving object 
In section 3 we derived the likelihood condition 
imposed by a single observation and discussed in 
section 4 the resulting information on the location of 
the moving object It was shown that the moving 
object was confined to a certain region in space with 
known boundaries. We finally consider all the n 
observations. Each observation determines a region; 
if there exists a point set in the parameter space 
common. to all these regions then a maximum likeli- 
hood region FR exists and is identical with this point 
set The boundaries ol R are hyperboloids of revo- 
lution. Any point inside the region PR can be con- 
sidered to be a maximum likelihood estimat« 

With an appropriate choice of stations it is easy to 
devise a computational procedure that leads always to 
an interior point of R. If the maximum likelihood 
region becomes small, then it is practically possible 
to use any point interior to PR as an estimate for the 
unknown position of the moving object 


The maximum likelihood region cannot increase in size when the number 
of observation nere ase On the contrary t to be expected that it will 


decrease 


Wasutneton, August 18, 1950 





Infrared Spectra of Eighteen Halogen-Substituted 
Methanes 


Earle K. Plyler and W. S. Benedict 


In order to study how the vibrational frequencies of methane vary with halogen sub- 


stitutions, the 
studied with a prism spectrometer between 
trifluoromethane, dichlorodifluoromethane, 
trichlorobromomethane, 
chlorofluoromethane, chloroform, 


9 


dichlorodibromomethane, 
dichlorobromomethane, 


infrared spectra of the following halogen-substituted methanes have been 
and 38 microns 
trichlorofluoromethane, 


Carbon tetrafluoride, chloro- 
carbon tetrachloride, 
difluorochloromethane, di 
bromo- 


fluoroform, 
chlorodibromomet hane, 


form, iodoform, methylene fluoride, methylene chloride, and methylene iodide Spectra 
are presented as obtained at room temperature in the vapor and iiquid states, and in a few 
cases in solution The bands are interpreted as fundamentals, overtones, combination, 


and difference bands 


fundamentals of these molecules and other halogen-substituted methanes 
using an approximate potential function with general constants transferred from 


calculations, 


one molecule to another are reported 


1. Introduction 


The infrared absorption spectra of various halogen 
derivatives of methane have been studied by several 
investigators in various regions of the spectrum, and 
under varying degrees of dispersion [1].' The Raman 
spectra of these compounds have also been exten- 
sively studied [2]. The gradual improvement in 
infrared techniques, in particular the development of 
prism spectrographs using thallium bromide-iodide 
(KRS-5) as the dispersing medium, and the increas- 
ing practical importance of the fluorine-containing 
compounds (Freons) have made it seem worth while 
to carry out a comprehensive reinvestigation of these 
compounds over the wide spectral range of 2 to 38 ug, 
using prism spectrographs of various dispersing inter- 
vals, so that good dispersion is obtained throughout 
the entire region. Compounds were studied both in 
the vapor and liquid states, using in many 
comparatively thick cells in order to develop bands 
of low intensity. 

A listing of the compounds studied, together with 
their boiling points and the phases in which they 
were observed, is given in table 1. Representatives 
of all types of halogen-substituted methanes except 
the methyl halides, which have been extensively 
studied under high resolution, were included 

Since the spectra were recorded with pris's spectro- 
graphs, it was not possible to resolve the rotational 
fine-structure. 

It should be possible to account for the observed 
absorption bands in terms of transitions among the 
various vibrational states of the molecules. The 
bands will be of three types: fundamentals, in which 
transition occurs from the ground vibrational level 
to a state in which one mode of vibration is singly 
excited ; overtones and combinations, from the ground 
level to states in which there are two or three quanta 
of vibrational energy in one or more modes; and 
difference bands, in which transition is from one 
singly excited vibrational mode to other higher vibra- 
tional levels. In general the fundamentals will be 


Cases 


' Figures in brackets indicate the literature references at the end of this paper. 


202 


| 
| 
| 
| 
| 
| 


lables and graphs display the regularities existing among the various 


Force constant 


the most intense bands, although in certain of the 
molecules, because of symmetry considerations, the 
intensity in the vapor phase should become zefo 
The combinations and overtones, lying at shorter 
wavelengths than the fundamentals, will be weaker, 
and the difference bands, which, in general, will be 
at longer wavelengths, will be very weak, except in 
the heavier molecules, where the fundamentals will 
be of such low energy that appreciable percentages 
of the molecules will be in excited states at room 
temperature The following sections describe the 
experimental method, present the results, and discuss 
further the vibrational assignments. 


2. Experimental Methods and Results 


The infrared absorption spectra of 18 substituted 
methanes have been measured in the region 2 to 38 
They were measured in the liquid state, at room 
temperature whenever possible, and also several of 
the compounds, as listed in table 1, have been meas- 
ured in both the liquid and vapor state. The eight 


Freons have been measured as gases, except tr- 
chlorofluoromethane, in the 23- to 38-g region. In 


order to observe the bands of low intensity in liquids, 
a cell with a thickness of 1.6 mm was used. For the 
intense bands a film of liquid between windows 
constituted the cell. In such cases the cell thickness 
was of the order of 0.01 mm. For the measurement 
of the gases and vapors the cell thickness varied 
from 5cemtolm. The pressure of the vapors was 
that of the saturated vapor at room temperature 
(approximately 25°C). For the Freons the pressure 
was atmospheric or less for the different cell lengths 
The liquid-cell thickness, gas-cell length, and gas 
pressure are given on graphs of the results. The 
spectra were recorded on a Perkin-Elmer spec 
trometer. 

In order to cover the region from 2 to 38 yg with 
good resolution, prisms of lithium fluoride, sodium 
chloride, potassium bromide, and thallium bromide 
iodide were used. Sharp bands could be measured 
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h an accuracy of 5 cm™~'! or less for the entire 
on. The general method of measurement was 
same as that used in previous studies [3] 

The eight fluormated 
plied by \\ ‘e 


methane derivatives were 
Murray, of Kinetic Chemicals. 
mington, Del. They were all above 99 percent 
purity, except carbon tetrafluoride, whose spec- 
mm showed one band that could not be accounted 
by combination or overtone bands. It appeared in 
same position as an intense band in chlorotrifluoro 
The carbon tetrafluoride gas 


thane at about 9 uw ‘ vr 
tested for purity by I L 


Mohler on the mass 
trometer and found to be 98.2. percent pure 
he Michigan Chemical Co. supplied the chloro- 
bromomethane, dichlorobromomethane. trichloro 
omomethane, and dichlorodibromomethane, and 
ye of the others were from Dow Chemical Co 
he bromoform and chloroform were obtained from 
istman Kodak Co. and were purified by W. Harold 
mith, of the Organic Chemistry Section After the 
st purification 
}4-u region in addition to the intense band at 
After the second purification, the band at 

{ uw disappeared 
lhe results are presented as graphs of percentage 


the chloroform showed a band in 


inclusive. The conditions of meas- 
rement are stated on each graph. In general, the 


ol ‘ 4 ‘ ‘ ‘ . ‘ —" ‘ ‘ ‘ 
2 ’ ‘ s ‘ ’ 7 
WAvVELEnety 


Fiaure |. 


Ihe pressure is atmospheric 


I'he insert was measured at reduced pressures 


Infrared spectra of gaseous CH,F,, 
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spectra are presented in groups of five molecules of 
similar structure, and the spectrum is divided into 
two portions 2 corresponding to measure 
ments with sodium chloride and lithium fluoride 
prisms, and 15 to 38 uw, corresponding to measure 


to 15 pg, 


ments with potassium bromide and thallium bromide 
iodide pristis 

Numerical values of th maximum 
absorption that may be identified as band centers, 
and further data pertaining to their interpretation 


posit 1OTLS of 


are presented in tables 2 to 6, inclusive. These 
include the band centers in wave numbers; the ap 
proximate intensity, characterized as weak, strong 


ete for a more quantitative estimate the yraphs 
should be consulted); the vibrational transition: and 
the svmmetry of the transition In addition to 
bands that appear in the infrared spectrum, there 
are given in parentheses the fundamental frequencies 
that do not appear either because they fall outside 
the range of wavelengths studied. o1 
have 


because they 
intensities. These fre 
quencies are obtained either from the Raman effect 
or from combination and difference The 
values given refer in general to the liqguic state, 
except for compounds 
temperature 


zero: ool very weak 


bands 
that are room 


where there is a pronounced 
difference between the two states will be mentioned 


rauscous al 


( ‘uses 


in the lise ussIoOn 


| 
A 4 ‘ a a ee ‘ ‘ ‘ ‘ | 
. ” " ” ) “ ” 
- wicrons 
CF,, CCl CC }4 and CCLEF from 2 to 14 pa. 





TRanswrss 


PERCEN 


100 


“a Enotes “CRON 
. 
7 
7 , 
Wavelenetn ” wCRONS 


one 
a4 
« . 
—————— | 
j 
. | 
5 4 s 
» ® roo 
- -_ ay Cy 
———4 
~~ 
~3 mm (CMa 
Zc 
‘a 4 
. / 
j 
‘ eae 
Oly 
rum 
_ ——_{ 
4 
4 leeeeeiemenea ieee 
ners 
Ome 
~—~ 4 
\ ; 
\ 
n “ ‘ ‘ SS 
O2e= Mly 
a5 /~ OF cS, 





spectra of ( 
CHBrCl,, 
CBrCh, « 
in the liau 


Fictre 3 
spectra of 


CH,l;, 





Infrared 
H Br, 
C Bré 
ind ( 
d state 


Infrared 
CH,Ch, 
CHC) 


CHBr,, in the q 
uid state and CHI 
in solution from 2 to 


15 w 




















WAVELENGTH 





FiGuRE 4 Infrared spectra of gaseous CUCI,F, CHCIF,, and CHF; fron lol 
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Figure 5 Infrared spectra of gaseous 


CH. F,, CF, CCIF,, CCL F>, and CChl 
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Figure 7. Infrared spectra 
] of CH,Ch, CH,I,, CHC! 
CHBrs, in the liquid state, 
and CHI, in solution fron 
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Ficure 9. Infrared spectra of vapors at saturation pressures at room temperature of ¢ HBr,Cl, CHBrC] 
CH,Ch, CH, Br:, and CBry, from 2 to 14 yu 


I'he insert curves are for reduced vapor pressures. 
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Figure 12 Infrared absorption spectra of satu- 


rated vapor of CCl, at room temperature 


Insert curve is for reduced pressure 
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Comparison of the positions of the fundamental bands of mole- 
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Fiaure 14. Comparison of the positions of the fundamental bands of mole 
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3. Discussion 


3.1. The Fundamental Frequencies 


Each of the related methane derivatives reported 
in this paper has nine fundamental frequencies of 
vibration. The general theory of the vibrations of 
such molecules has been discussed in many places 
and need not be repeated here. It will suffice to 
recall that when there are atoms of different kinds 
joined to the carbon atom, the symmetry of the 
molecule is less than when all four linked atoms are 
alike, and that the particular degree of symmetry 
yoverns the degeneracy of the V ibrations, the general 
form of the molecular motions, and the resulting 
activity in both the infrared and the Raman spec- 
trum. The pertinent information derivable from 
such symmetry considerations is summarized in table 
7. The symmetry class of overtones, combinations, 
and difference bands of the various fundamentals 
follow well-known rules and are given in tables 2 to 6. 

The numerical values of the fundamentals show 
regular variations when going from one molecule to 
another in which a single substituent has been 
altered. This is a qualitative consequence of the 
facts that in all these substituted methanes the 
forces between the central atom and a particular 
peripheral atom are to a first approximation inde- 
pendent of the other peripheral atoms, and that 
these forces vary regularly with the mass of the 
peripheral atom. A more quantitative discussion 
will be given later. The relative positions and 
groupings of the numerical values are also affected 
by the symmetry. The conventional system of 
numerical identification of the fundamentals is based 
on the symmetry. This system fails to display the 
regularities existing among molecules of different 
symmetry. Accordingly, in this paper, one purpose 
of which was to show the consistency of all the 
frequencies of a large number of related compounds 





of different symmetry, use is made of an unconven- 


tional numbering system. The system closely re- 
sembles the conventional one for molecules of the 
tvpe CH,X., where there are no degeneracies. Fo; 
molecules of higher symmetry than C,,, where de- 
generacy occurs, the degenerate frequencies have 
been assigned two or three numbers, so that all nine 
The relation between 


frequencies may be identified 
would 


the numbers adopted here and those that 
follow from application of the conventional rules js 
shown in table 8. This table also collects in one 
place the numerical values of the fundamentals of 
all the molecules, the spectra of which are reported 
in this paper, as well as of related molecules that 
have recently been studied in this laboratory [4] 
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The regularities mentioned in the preceding para- 
graph are displayed in figures 13 to 15, inclusive. 
In these graphs the frequencies of related molecules 
are plotted, the molecules being arranged vertically 
in order of increasing molecular weight. The regular 
decrease in frequency with increasing weight of the 
substituent, and the way in which the degenerate 
frequencies in the more symmetrical molecules break 
up in the less symmetrical ones and then regroup, is 
to be noted. 

Another often-used method of identifying the 
fundamentals, which also possesses a degree of physical 
significance, is to characterize them by the type of 
intramolecular motion involved. In this method, 
v, signifies the stretching of a bond connecting the 
carbon atom with the x substituent, and 6, signifies 
the bending of such a bond. All the molecules con- 
sidered will have four »’s and five 4’s. When there 
are several identical bonds, the suffix s or a may be 
added, according to whether the motion of the 
identical atoms is symmetrically toward the central 
atom or otherwise. In molecules as light as those 
considered here, this method often has limited mean- 
ing, since the motions in each fundamental mode 
cannot be localized as closely as the method would 
imply; in the 6 modes in particular all the atoms 
participate in the motion. Nevertheless, a charac- 
terization of the fundamentals by this method is also 
included in table 8 and will be referred to in later 
discussion. For some such cases no suffix can be 
given; in others, a double suffix identifies the two 
atoms most prominently involved. 
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The relative intensity of the various fundamentals 
is not displayed in figures 13 to 15. The intensity 
depends on a number of factors, principally the 
symmetry of the motion and the distribution of 
electrical charge within the molecule. From the 
spectral graphs and tables it will be noted that there 
is a general parallelism in the intensities of corre 
sponding bands in the various molecules of like 
symmetry. It should be noted that, although for 
molecules of symmetry C,,, the A, frequency « 
should be inactive, it does appear in all the molecules 
considered, if it falls in the observable frequency 
range. It is always quite weak, however, and its 
appearance may readily be explained as due to a 
breakdown of the strict selection rules because of 
intermolecular forces, particularly in the liquid state, 
or to Coriolis interactions with nearby active fre- 
quencies, or to a combination of these causes 


3.2. Force Constant Calculations 


For a quantitative verification of the assignment 
of fundamentals, recourse must be had to a calcula- 
tion in terms of the atomic masses, atomic dimen- 
sions, and a suitable set of force constants of the 
potential energy. In the most general potential 
function there are many more disposable constants 
than there are observed frequencies, so that various 
approximations must be introduced. Certain of the 
constants, especially those involving interactions be 
tween motions of different types, may be set equal 
to zero; other constants may be assumed to have th: 
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Many 
h treatments have been reported for most of the 
lecules considered here 
of the recent 

yecius 15), who was able to calculate 82 frequencies in 


ne numerical value in different molecules 


One of the more success 


{normal and deuterated chloro-bromo-methanes from 
with a mean accuracy of 1.1 percent 
In order to substantiate some of our assignments. 

id to determine whether a similar treatment may 
applied to compounds containing fluorine, the 
ilculations of Decius have been extended to the 
iolecules CF,, CF,H, CF,H,, CFH,, CF,Cl, CF.CL, 
PCL, CCLBr, CCLBr., and CCIBr,, for which all 
ne frequencies were computed; and to the less 
ymmetrical molecules CHF,Cl CHFCL, CHCIBr 
HCIBr., CH.FCI, and CH,CIBr, for which only the 
hree frequencies of symmetry A’’ were computed 
lo evaluate the remaining six frequencies of the 
atter molecules would involve the expansion and 
solution of a sixth-order determinant In these 
aleulations the constants and molecular dimensions 
of Decius were adopted when applicable (one ap 
arent numerical error or misprint, for /} in C—Cl 
New values for the corresponding 
constants involving fluorine were adopted, which 
rave a fair over-all fit of the data considered It is 
probable that additional computations would yield 
somewhat better values of these constants, but this 
does not seem useful, in view of the approximations 
nvolved. The constants used are listed in Decius’ 
notation in table 9. Calculations on “‘mixed”’ mole 
that is 


Ss constants, 


was corrected 


iles those 


containing more than one 
halogen, require additional constants For these 
are used the geometric mean of the corresponding 
PARLE ¥ Force constant for halomethanes (afte Deciu 
A g ta r r : j 
erical u uw ‘ retet ' " l 


r bending and stretching | ) ! ec for 


simpler calculations was that of 


constants involving a single halogen 


For example, 


since the constant fx, corresponding to the 
interaction of motion along two C—X'! bonds, 
C—X'-C— Xi, is 0.96 (in units of 10° dyne/em) for 
X F and 0.332 for X Cl, the constant for inter 
action of motion along C—F and C—Cl would be 
0.96) - (0.332)] 0.565. The constants involving 
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bending of X'—C—X' bonds and their interactions 
are obtained similarly. without introducing further 
hew constants 

The resulting computed frequencies are compared 
with the observed fundamentals in table 10 It is 
seen that, in general, the agreement is good. The 
greatest deviation is percent, and the average 
deviation of 85 frequencies in 16 molecules is 2.22 
percent. This degree of fit was obtained by intro- 
ducing 13 new constants, in addition to the 28 taken 
from Decius Better agreement is not to be exper ted 
in view of the approximations introduced. Tnasmuch 
as the very accurate results of microwave 
troscopy show that the interatomic distances are not 
constant, particularly for the C—F bond, and that 


spec- 


the angles are not tetrahedral, it is clear that the 
force constants must vary to some extent. The 
interaction constants are relatively large: the par- 
ticularly high value of f, for C—F (0.96% 10° dyne 


em is to be noted This is a further reflection of 
the importance of specific forces arising from the 
highly polar C—F linkage, and indicates that for a 
truly accurate calculation, one should take into 
account all the ignored interactions, perhaps basing 
them on some more general force field [6]. The 
present comparison serves merely to confirm the 
choice of fundamentals and to demonstrate the ap- 
proximate similarity of the carbon-halogen bonds in 
a large number of halomethanes 

It should be noted that the C—F constants, par- 
ticularly the major ones involving the stretching, are 
in good agreement with a very recently published 


ealculation based on much the same data and 
methods {7}. It should also be mentioned that pre- 
liminary calculations indicate that the C—F con- 


stants used also give adequate agreement with ob- 
servation in such molecules as C,F,, CF;—CHsg, 
CF,—CCl,, ete. 

There are not enough data on the iodine com- 
pounds to derive constants for this bond, but ap- 
proximate values extrapolated from the 
F—Cl—Br--I may be seen to give adequate agree- 
ment with the few observations. 


series 


3.3. Overtone, Combination, and Difference Bands 


In addition to the fundamentals, the bands that 
appear in the spectrum must be attributed to over- 
tones, combinations, and differences of the funda- 
mentals. They have been so assigned in tables 2 
to 6. The assignments appear eminently reason- 
able, in that the selection rules are obeved, in the 
molecules of higher symmetry, with the exception of 
an occasional weak combination of symmetry A, in 
C,,; that binary combinations are in general much 
stronger than those involving three frequencies; and 
that difference bands arise only from relatively low 
vibrational frequencies. 

Under the dispersion of prism instruments, it is 
not always possible to make an unambiguous assign- 
ment of these higher-order bands. As will be noted 
in tables 2 to 6, there are many instances in which 
two or more combinations can be offered as probable, 
and many other such have not been listed separately. 


| 


The magnitude of the anharmonic corrections, by 
which the frequencies of the overtones and combina- 


tions differ from the sum of the fundamentals. 
ranges from +10 to 50 em and, with nine 
fundamentals, such overlapping is common. Unless 


much higher resolution is available, it is unprofitable 
to attempt to distinguish between the possibilities, 
except to note that the assignments that are set 
forth show a consistent behavior in a large number 
of molecules. Indeed, because of the possibility of 
Fermi resonance among levels of like symmetry, it 
is often meaningless to distinguish between close- 
lving levels. Some of the more striking and definite 
examples of this effect occur in CCL, and CCLF, and 
will be discussed in more detail later. In a very 
heavy molecule of low svmmetry, such as CCLBr,, 
the number of possible difference bands of appreci- 
able intensity arising from the five fundamentals 
below 300 cm@' is particularly large. It will be 
noted in figure 6 that in these molecules there is a 
rather strong background of structureless absorption 
throughout the region 280 to 500 cm which un- 
doubtedly arises from such overlappings. 

It should be mentioned that a number of the over- 
tone and combination bands of the methylene-halide 
compounds CH,F,, CH,FCI, CH,Cl,, CH,CIBr, and 
CH,Br, have been studied under very high resolving 
power in the region 1.6 to 2.5 u. The results will 
be reported separately. They show the expected 
large number of bands, with the rotational structure 
partially resolved, and further confirm the vibra- 
tional assignments reported in this paper. 


3.4. Liquid-Vapor Shifts 


As mentioned above, a number of the compounds 
of low boiling point have been studied both as liquids 
and as vapors, at room temperature. As is com- 
monly observed, the vapor spectra show more in- 
dications of partially resolved rotational structure 
(P, ©, and R branches), but otherwise do not differ 
appreciably from those of the liquids, except for 
slight shifts of the central absorption frequencies. 
Some of these shifts are listed in table 11. It will 
be noted that the frequencies for the vapor are con- 
sistently higher than for the liquid, and that the 
magnitude of the shift is in general greater, especially 
on a percentage basis, for the stretching frequencies 
involving the halogens (»,’s) than for the bending 
frequencies or for the hydrogen-stretching fre- 
quencies. This is undoubtedly due to the fact that 
the shifts result from partial association in the 
liquid, which is predominantly due to forces between 
the dipoles localized in the carbon-halogen bonds. 

The fluorine compounds CH,F, and CHF, were 
studied in this research only in the gaseous state. 
Their Raman spectra have, however, been studied in 
the liquid. Since the C—F bond is considerably 
more polar than C—Cl, one might expect the vy fre- 
quencies to show quite large liquid-vapor shifts, and 
these have been observed. hese two molecules 
have also been listed in table 11, the A now referring 
to the shift between our adopted infrared frequencies 
and the most recent Raman data [8]. 
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3.5. Discussion on Particular Molecules 


It does not seem necessary to discuss all the assign- 
ments for each mole ule, since In most cases the 
assignments are in substantial agreement with the 
work of other investigators who have studied these 
molecules, either by the infrared or Raman spectra 
In a few instances, however, where certain points 
have been discussed in the literature, it seems desir- 
able to justify our choices further It should be 
pointed out that in many cases final decisions can 
hardly be reached until the rotational fine structure 
of the bands in question is studied under very high 
resolution 


3.6. CF, 


The bands may all be assigned as F, combinations 
or overtones of the two infrared active frequencies 
629 em and 1,277 em™', and the Raman-active 
137 em™' (E) and 904 em™ (A,). The fact that 
these pairs are each nearly in the ratio 1:2 would 
make it possible for Fermi resonance to appear in 
many of the combinations; however, as the observa- 
tions give a band of high intensity at each position, 
it appears that the interactions are not strong, and 
that the resonance is only partial. There is a very 
weak peak at ~1,238 cm™', which may correspond 
to the weaker component of the 29, rigs Fy pair; the 
corresponding combinations with », also appear at 
2,123 and 2,160 with a relative intensity of about 
1:5. The unresolved rotational envelopes of the 629 
and 1,277 em~' fundamentals are quite different, the 
former showing well separated P, ©, and # branches, 
with a P—R distance of ~33 em™', whereas in the 
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latter this distance is much less and unresolved 
This is due to the interaction with the rotation of 
angular momentum of vibration, and is in qualitative 
agreement with calculation: the ¢-factor for the 
valence vibration ys. should be about +-0.68 and for 
Yeo about 0.18, so that the rotational spacings for 
the two bands should be in the ratio (0.32/1.18) ~0.8 
The combination bands of these fundamentals with 
vy Show the same behavior; for example, 1,532 is 
found with a P—R spacing of ~38 em 
2.160 is narrow and unresolved 


3.7. CF,;Cl 


The two lowest-frequency fundamentals, w;(E) 
and w,(A,) are very weak in the infrared, but appear 
strongly in the Raman spectrum of the liquid. Al 
though a weak band appears at 450 em~', in the 
infrared spectrum, it seems preferable to adopt 
170 cm"! of the Raman spectra for the vapor-phase 
frequency of w,, as this accounts for the weak ab- 


whereas 


sorption at 635 em~! as a difference band and gives 
satisfactory agreement in a number of combinations. 


3.8. CF,Cl, 


The foree-constant calculations definitely require 
that there be five (deformation) frequencies below 
500 em”! and only four above 600 em Hence the 
assignment of both the strong bands at SS2. and 
922 em™' as fundamentals cannot be correct; it is 
very reasonable to consider this a case of Fermi 
resonance, with the original w, fundamental (B,) at 
906 em! split by interaction with the combination 
w,(A,)455+-,(B,)437. The region from 420 to 480 
em”! must be studied under higher resolution before 
the values for the three overlapping frequencies wy», 
Wy, and @, Can be considered certain The value of 
155 em”! for w, is from the Raman effect 


3.9. CFCI 


As in the other chloro-fluoromethanes, the de 
formation fundamentals are very weak, but Raman 
effect measurements enable a choice to be made 
between the fundamentals and the equally intense 
difference bands at wavelengths greatet than 20u 
The combination band at 932 em™' is more intense 
than the fundamentals on which it is based. This 
may result from a partial Fermi resonance with the 
very intense agg. The vp frequency shows a + 17 em 
shift from the liquid value found in the Raman 
effect [9] 

3.10. CCL, 


This much-studied molecule shows one of — the 
earliest-recognized examples of Fermi resonance in 
the strong pair at 762 to 785 em™', The unsplit 
position of the apg is taken as the exact center, 
773 em since the two bands are nearly of equal 
intensity, in the liquid. In the vapor the resonance 


with 3194-458 is less exact, the intensity ratio now 
being about 4:1, the higher frequency component 
being the stronger. This is in accord with the general 
observation that the vapor-liquid shift is greatest for 
the asymmetric valence vibrations. It is of consider 








able interest that the doubling also appears in two 
difference bands, arising from w.; and wr 


3.11. CCl,Br and CCl,Br, 


In these heavy molecules many fundamentals lie 
beyond the limit of our observation, but are known 
from the Raman effect. There is considerable ab- 
sorption without definite structure beyond 20 yg, 
corresponding to many possible overlapping differ- 
ence and combination bands of the deformation fre- 
quencies, 

3.12. CHF 

This molecule has been studied under very high 
resolution in the photographic infrared [10]. We 
have also examined a number of the combination 
bands under high resolution. The only question in 
the frequency assignment concerns the value of as, 
the vp,, fundamental. Its intensity and polarization 
in the Raman spectrum definitely prove that in the 
liquid it fais at 1,116 em In the infrared spec- 
trum of the gas there is no indication of any weak 
absorption to longer wavelengths of the 1,152 em”! 
peak. It is therefore believed that the vp,, and vy,, 
fundamentals overlap at this point, the difference 
from the Raman effect being due to the vapor-liquid 
shift, whereas the weak band at 1,209 em“ is a 
combination tone. <A high-resolution study of this 
region would be of interest. Our force-constant cal- 
culations, in agreement with those of Rank, Shull, 
and Pace, confirm the assignments of 700 cm™! to A, 
and 507 em™' to E symmetry, although the rotational 
envelopes of the two bands are nearly identical. 


3.13. CHF,Cl 


The assignment of the ws(ve,,) fundamental differs 
from that of Glockler, Edgell, and Leader [11] ob- 
tained from the Raman effect. There are two me- 
dium-intensity bands to choose between, at 833 and 
1,178 em". The higher value fits much better with 
the related molecules. The low frequency, which is 
fairly strong in the Raman effect, and which may 
well be accounted for as 2422 em, undoubtedly 
gains in intensity by partial Fermi resonance with 
w;, at 809 em™. 


3.14. CHFCl,, CHCl,, CHC1,Br, CHCIBr,, CHBr, 


In all these molecules, there is little question as to 
the assignment of fundamentals, when the infrared 
and Raman results are both considered. All show 
absorption in the long-wavelength region. The 262 
em™~' fundamental of CHC, is relatively strong, but 
the long-wavelength limit of the 


3.15. CHI, 


The results for this molecule are quite incomplete. 
Because of its physical properties (solid at room 
temperature and above, subliming with decomposi- 
tion) it could be studied only in solution. Its 
chemical instability prevented the use of high con- 
centration, and rapid surveys of the spectrum were 
required. The Raman frequencies, and rough force- 
constant calculations, show that two of the funda- 
mentals lie at very long wavelengths. The few 


occurs just at 


T1Br-I prism. 
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bands observed find a natural interpretation by 
analogy with the corresponding molecules. 
3.16. CH,F, 

This molecule has been investigated recently under 
high resolution [12]. As mentioned above, high- 
resolution studies in the 1.6 to 2.5 u, region have also 
been made in the radiometry laboratory of the Bureau 
The assignment of fundamentals differs from that of 
the above authors, but is in agreement with the one 
very recently offered by Rank, Shull, and Pace, [8} 
and with the force-constant calculations. Th: 
Raman-effect data leave no doubt that in the liquid 
w is at 1508 em"'. This vy, frequency might well 
be very weak in the infrared and escape detection, 
since it falls in a region where the H,O absorption 
masks weak bands. The band structures found by 
Stewart and Nielsen are equally well accounted for 
by the assignment, since the B, combination at 
2,945 em™'! can be ascribed to w;+,. As mentioned 
above in connection with table 11, a relatively large 
liquid-vapor shift may be invoked to reconcile the 
Raman and infrared measurements in the region of 
1,100 em~'. The depolarized Raman line at 1,054 
em~! can be attributed to ws , rather than to 2a, 
which should be polarized. The w, frequency, whose 
A, symmetry should make it inactive, appears 
weakly at 1,262 cm~', presumably due to Coriolis 
interaction with w, (B.) at 1,176 em™', or to break 
down of the strict selection rules induced by the 
strong polar intermolecular forces. There are two 
weak bands, at 816 and 1,349 em~! which cannot be 
They may be due to impurities. 


3.17. CH,Cl, and CH,I, 

These molecules show no unexpected effects 
There is a close parallelism among the spectra of all 
the methylene halides. These will be discussed at 
greater length in connection with high-resolution 
studies in the overtone region. 


me 
assigned. 


The authors are indebted to N. Acquista and 
Mary A. Lamb for assistance in the observations of 
the spectra. 
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